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ABSTRACT 


An  analytical  and  experimental  study  of  ship  motions 
and  capsizing  in  extreme  seas  has  been  conducted.  The  analysis 
of  linear  and  quasi- linear  one-dimensional  roll  models  has 
revealed  motion  anomalies  not  apparent  from  the  usual  linearized 
ship  motion  theory.  Extensive  tests  have  been  conducted 
using  two  radio-controlled  models  in  the  wind  generated  seas 
of  San  Francisco  Bay.  Directional  spectra  were  computed, 
using  ?,  variety  of  techniques,  from  the  wave  measurements  by 
an  array  of  wave  gages.  Comparisons  of  the  experimentally 
determined  motions  and  a linear  strip  theory  prediction  are 
presented.  A time  domain  numerical  simulation  program  for  motions 
and  capsizing  has  been  used  to  investigate  motions  in  a variety 
of  wave  group  configurations.  The  results  show  good  agreement 
with  observed  capsizing  phenomena  and  have  revealed  a number  of 
important  characteristics  associated  with  large  geometry  changes 
in  waves. 
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I. 


INTRODUCTION 


Stability  against  capsizing  in  heavy  seas  is  one  of  the  most 
fundamental  requirements  considered  by  the  naval  architect  when 
designing  a ship.  At  the  same  time  it  is  one  of  'he  least 
understood  phenomena  in  the  realm  of  ship  motion  theory.  By 
its  very  nature,  stability  is  concerned  with  extreme  motion 
both  of  ship  and  waves.  Therefore,  the  linearized  methods  of 
wave  and  ship  motion  analysis  which  have  received  so  much 
attention  in  recent  years  cannot  be  applied  here  with  any 
expectation  of  obtaining  directly  useable  results,  but  may  be 
used  only  for  general  guidance.  Statistically,  capsizing  falls 
into  the  category  of  events  of  rare  occurence,  again  calling 
for  special  techniques  in  studying  its  probability  distributions. 

The  basic  purpose  of  studying  capsizing  is  to  establish 
prediction  techniques  which  may  be  used  by  naval  architects 
during  the  design  process.  Even  though  traditional  empirical 
methods  have  apparently  been  successful  in  insuring  adequate 
stability  of  conventional  ships  built  in  the  past,  it  is 
expected  that  the  future  will  see  more  rapid  evolution  in  ship 
size,  speed,  geometry,  and  performance  requirements.  Some  of 
the  new,  high  performance  ship  types  may  be  expected  to  experience 
motion,  control,  and  even  stability  problems  when  operating  in 
severe  seas  which  are  not  susceptible  to  prediction  or  remedy 
by  traditional  empirical  methods.  The  more  general  goal  of 
our  study,  therefore,  is  to  establish  a basis  of  understanding 
of  phenomena  associated  with  the  operation  of  ships  in  severe 
seas  and  to  relate  these  ship  behavior  phenomena  to  the 
geometrical  and  operational  ci.°»"acteristics  of  the  ship. 

Model  experiments  conducted  under  conditions  which  are  as 
realistic  as  possible  in  simulating  a severe  storm  sea  environ- 
ment offer  the  best  initial  means  of  gaining  this  understanding . 

At  the  present  time,  there  are  no  towing  tanks  which  have 
the  capability  of  producing  realistic  short-crested  seas  of 
sufficient  severity.  Further,  the  infrequent  nature  of  the 
capsizing  event  indicates  that  experimental  runs  of  considerable 
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length  will  be  required,  especially  in  following  seas  where 
the  frequency  of  wave  encounter  is  low.  The  testing  of  free- 
running  radio-controlled  models  in  an  open  water  area  in  natural 
wind-generated  waves  is  an  alternative  to  tank  testing  that  has 
beer  successfully  used  for  this  purpose  in  the  past. 

As  is  the  case  of  all  experiments  conducted  in  a natural 
environment,  there  is  the  difficulty  of  controlling  the  conditions 
of  the  test.  However,  the  two  principal  requirements  enumerated 
above:  steep  short-crested  waves  and  virtually  unlimited  length 

of  run  are  available  only  in  this  way. 

The  first  experiments  of  this  type  were  conducted  by  a group 
working  under  the  direction  of  Professor  Kurt  Wendel  at  the 
universities  of  Hamburg  and  Hannover  in  Germany.  Their  work 
clearly  pointed  out  the  value  of  open  water  testing  as  a means 
of  studying  capsizing.  Partly  as  a result  of  the  success  of 
Wendel,  the  present  capsizing  study  was  initiated  in  1969  at 
the  University  of  California  under  the  sponsorship  of  the  U.S. 
Coast  Guard.  riMs  research  started  with  a program  of  model 
testing  in  the  open  wateis  of  San  Francisco  Bay.  Concurrently 
with  the  model  experiments,  analytical  work  has  been  carried 
out  in  two  general  areas: 

1.  Work  supporting  the  experiments  themselves,  e.g.,  the 
development  of  procedures  for  data  processing. 

2.  Studies  of  certain  aspects  of  ship  motion  complementary 
to  the  experiments  as  an  aid  in  the  design  of  the  experiments 

or  as  a means  of  explaining  experimental  results.  Especially 
noteworthy  of  work  conducted  in  the  first  category,  a thorough 
evaluation  ht.s  been  made  of  available  techniques  for  the 
directional  analysis  of  wave  spectra  using  measurements  made 
with  sparse  wave  gage  arrays.  This  has  led  to  the  development 
of  a new  method  of  analysis  of  directional  spectra. 

In  the  second  category,  several  reports  (in  the  form  of 
student  theses)  have  investigated  areas  of  nonlinear  ship  rnotic 
analysis  by  either  new  techniques  or  in  greater  detail  than 
be  lore  in  an  attempt  to  obtain  insight  and  understanding  of 


ship  behavior  anomalies  observed  experimentally.  Some  or  these 
results  may  be  found  in  Haddara  (1971),  Perez  (1972)  Oakley  (1973), 
Paulling  and  Wood  (1973),  and  later  sections  of  the  present 
report.  A paper  presented  before  the  Tenth  Naval  Hydrodynamics 
Symposium,  Oakley,  et  al.  (1974)  forms  the  nucleus  for  the 
present  report. 

Experiments  have  now  been  completed  with  two  models,  one 
a conventional  cargo  ship  and  the  other,  a fast  twin  screw 
container  ship.  As  a result  of  these  experiments  and  concurrent 
analytical  work, a great  deal  has  been  learned  about  the  mechanism 
of  capsizing.  This  understanding  has  led  to  the  development  of 
a numerical  simulation  of  the  large  amplitude  motion  of  the  ship 
in  certain  capsizing  situations  . It  is  now  possible  to 
duplicate  on  a computer  certain  important  aspects  of  the  motion 
of  the  model  observed  during  the  experiments.  The  numerical 
computation,  however,  has  one  drawback  in  common  with  other 
time-domain  simulation  techniques  — it  has  an  exorbitant  appetite 
for  computer  time.  It  does  not  appear  economically  feasible 
using  present-generation  machines  to  numerically  simulate  a 
complete  experimental  run,  not  to  mention  a ship  lifetime  of 
experience.  However,  an  important  observation  made  during  the 
experiments  suggests  how  this  difficulty  imy  be  overcome.  It 
was  observed  that  severe  motions,  control  difficulties  and  cap- 
sizings  of  the  model  almost  invariably  occurred  when  the  model 
encountered  a wave  group  consisting  of  several  especially  steep 
and  nearly  regular  waves.  Such  wave  groups  characteristically 
occur  from  time  to  time  in  the  midst  of  a random  seaway.  While 
the  statistics  of  their  occurrence  apparently  have  not  been 
extensively  studied,  there  is  some  indication  that  it  is  possible 
to  derive  the  expressions  relating  the  frequency  of  occurrence 
of  groups  having  specified  characteristics  to  the  energy  spectrum 
of  the  seaway. 

The  experimental  observation  then  suggests  that  the  proba- 
bility of  capsizing  may  be  related  to  the  probability  of  the 
ship  encountering  a wave  group  having  the  characteristics 
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necessary  to  cause  capsizing.  The  numerical  simulation  technique 
noted  above  would  be  used  only  to  study  the  ship  behavior  in 
wave  groups  consisting  of  a small  number  of  individual  waves  in 
order  to  define  the  critical  group  characteristics  for  capsizing. 
The  long  term  statistics  associated  with  the  prediction  of 
capsizing  would  then  be  divorced  from  the  mechanism  of  capsizing 
and  would  be  concerned  only  with  the  prediction  of  the  probability 
of  the  ship  encountering  wave  groups  having  the  critical 
characteristics.  Thus  long  computer  simulation  runs  would  be 
unnecessary. 

Since  the  problem  of  predicting  the  occurrence  of  wave 
groups  has  received  so  little  attention  in  the  past,  it  forms 
a subject  for  study  under  the  present  work.  Evaluations  of 
some  previously  proposed  prediction  techniques  are  given  using 
experimental  wave  data  obtained  in  the  present  study.  Also, 
some  suggestions  are  made  for  the  future  direction  which  this 
work  should  take. 
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II.  LINEAR  AND  QUASI-LINEAR  MOTION 

COMPUTATIONS 

As  noted  in  the  INTRODUCTION  modern  ship  motion  research 
has  been  mostly  concerned  with  small  amplitude  motion  which  is 
amenable  to  a linearized  dynamic-hydrodynamic  analysis.  The 
The  practical  implementation  of  such  procedures  has  been  most 
successfully  obtained  in  the  strip-method  of  computation  as 
exemplified  by  the  work  of  Salvesen,  Tuck,  and  Faltinsen  (1970). 
It  is  safe  to  say  that  the  strip  calculation  now  provides  a 
workable  design  tool  for  predicting  the  average  seakeeping  per- 
formance of  a proposed  ship  early  in  the  design  process. 

Performance  predictions,  however,  are  usually  concerned 
with  the  average  behavior  of  the  ship  over  long  periods  of  time 
encompassing,  e.g.  , the  length  of  a voyage  and  tens  of  thousands 
or  hundreds  of  thousands  of  wave  encounters,  rather  than  the 
short  intervals  which  the  ship  spends  in  extremely  severe  seas. 
Performance  characteristics  which  have  been  successfully  de- 
scribed by  linear  procedures  are  ship  motions,  structural  loads, 
stresses  and  even  occurrence  of  seemingly  nonlinear  large- 
amplitude  phenomena  such  as  the  frequency  of  slamming  and  bow 
immersion.  Prerequisite  to  a study  of  performance  but  not  part 
of  it  is  survivability.  In  contrast  to  the  average  conditions 
with  which  performance  estimates  are  concerned,  a study  of  sur- 
vivability addresses  itself  to  the  extreme  sea  and  motion 
conditions  having  a very  low  probability  of  occurrence.  Thus, 
by  definition,  in  studying  capsizing  of  ships  in  extreme  seas, 
we  are  concerned  with  the  largest,  steepest  waves  which  the 
ship  may  encounter  during  her  lifetime  (sometimes  at  its  very 
end!),  and  with  motions  of  such  large  amplitude,  that  they  are 
well  beyond  the  range  of  a linear  ship  motion  analysis. 
Nevertheless,  it  is  useful  to  review  some  of  the  features  of 
linear  motion  theory  in  hopes  that  its  results  may  provide 
some  guidance  and  ir sight  into  certain  phenomena  of  concern  in 
a study  of  capsizing.  Further  insight  and  new  phenomena  may 
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also  be  revealed  by  certain  nonlinear  theoretical  procedures 
which  have  been  developed  to  augment  linear  ship  motion  theory. 

In  a study  of  capsizing,  what  results  of  interest  may  be 
obtained  from  linear  motion  theory?  The  first  which  suggests 
itself  is  resonant  rolling.  Although  a linear  analysis  of 
roll  does  not  accurately  describe  the  details  of  the  motion 
at  large  amplitudes,  some  of  the  properties  of  this  mode  of 
motion  may  be  revealed.  For  example,  there  is  substantial 
evidence  that  rolling  motion  in  beam  seas  in  the  absence  of 
other  effects  is  not  likely  to  result  in  danger  to  the  ship. 

The  vessel  may  roll  to  very  large  angles  if  the  frequency 
content  of  the  sea  is  in  the  right  relationship  to  the  natural 
frequency  of  the  ship,  but  in  order  for  her  to  be  placed  in 
danger  of  capsizing,  additional  factors  must  be  present. 

These  may  include  strong  wind  heeling  moments,  damage  to  the 
watertight  envelope  of  the  ship  resulting  in  the  flooding  of 
a portion  of  the  interior  and  the  entrapment  of  large  quantities 
of  water  on  deck  by  bulwarks. 

Tests  of  the  range  of  applicability  of  the  linear  ship 
motion  technique  have  been  made  by  Dalzell,  and  some  of  his 
results  for  the  Mariner-type  ship  were  reported  in  a discussion 
to  a paper  by  Ochi  (1964).  Pitch  and  heave  motions  were 
measured  in  regular  waves  of  various  steepness.  For  waves  more 
than  1.75  times  the  model  length  these  motions  show  a linear 
behavior  with  wave  height  for  values  of  wave  height  to  length 
ratio  up  to  1/12.  In  the  shorter  waves,  the  pitch  motion 
appeared  to  drop  somewhat  below  a straight  line  relationship 
for  the  higher  waves,  i.e.  the  graph  of  pitch  versus  wave 
height  was  concave  downward  in  these  waves. 

Figure  II-l  shows  some  results  obtained  at  the  University 
of  California  from  experiments  and  corresponding  calculations  of 
the  relative  rolling  motion  between  the  ship  model  and  wave 
surface.  The  experiments  and  calculations  were  carried  to 
moderately  large  angles  and  show  reasonably  good  correlation. 

The  important  characteristic  of  this  motion  is  its  tendency  to 
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Figure  II-l 


level  off  in  increasingly  higher  sea  states.  Since  the  relative 
rolling  motion  is  more  important  than  the  absolute  roll  in 
relation  to  such  effects  as  water  entrapment  on  deck,  this  suggest 
that  the  effects  of  increasingly  severe  sea  states  tends  to 
level  off. 

These  calculations  were  made  using  a linear  approximation 
to  the  ship  righting  arm  curve  equal  to  displacement  times  the 
initial  metacentric  height.  The  quadratic  roll  damping  was 
approximated  by  an  equivalent  linear  damping  coefficient.  The 
trend  of  this  curve  illustrates  the  degree  of  agreement  that  may 
be  expected  in  comparing  linear  calculations  with  experiments. 

The  maximum  wave  steepness,  however,  is  not  in  the  extreme  or 
survival  range. 

It  is  also  useful,  at  this  point,  to  look  at  some  results 
of  nonlinear  ship  motion  theory.  Two  topics  of  ship  dynamics 
have  been  studied  using  nonlinear  equations  of  motion:  maneu- 
vering in  calm  water,  and  rolling  in  btci.i  seas.  In  both  cases. 
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tfie  procedure  has  been  to  include  terms  in  the  equations  of 
motion  representing  forces  which  depend  in  a nonlinear  manner 
upon  some  of  the  unknown  motion  variables.  Most  commonly,  these 
have  been  damping  forces  assumed  to  vary  as  the  square  of 
velocity  or  more  general  forces  which,  by  reason  of  symmetry  , 
behave  as  odd  functions  and  may  therefore  be  expressed  by 
cubic  terms  in  all  velocities  and  accelerations.  In  principle, 
these  terms  represent  higher  order  terms  in  a Taylor  series 
expansion  of  the  total  force  on  the  ship,  where  the  force 
is  assumed  to  depend  in  some  unknown  way  upon  the  relevant 
instantaneous  motion  variables.  Many  of  the  terms  can,  in 
fact,  be  related  to  specific  physical  effects,  e.g.,  the  above 
noted  viscous  drag  which  varies  in  a nonlinear  manner  with 
velocity. 

The  results  of  including  nonlinear  terms  in  the  equations 
of  motion  are  twofold  . First,  behavior  which  might  be  obtained 
by  a conventional  linear  analysis  is  modified  somewhat  with- 
out changing  its  basic  nature.  Second,  new  phenomena  may 
be  revealed  which  are  not  present  in  the  linear  results.  We 
shall  look  at  examples  of  each  of  these  cases  and  attempt  to 
relate  them  to  the  capsizing  problem. 

In  studying  the  large  amplitude  rolling  motion  of  ships  in 
beam  seas,  the  method  of  equivalent  linearization  has  been 
utilized  by  several  investigators,  e.g.  Vassilopoulos  (1971)  and 
St.  Denis  (1967) . Equivalent  linearization  is  a method  which 
is  generally  suitable  for  describing  a dynamic  system  in  which 
large  deviations  from  linear  behavior  are  not  anticipated, 

(weak  nonlinearity) . A reasonable  approximation  to  the  exact 
behavior  of  the  real  system,  therefore,  is  given  by  an  equiv- 
alent linear  system  having  linear  coefficients  appropriately 
selected.  The  principal  shortcoming  to  the  method,  in  common 
with  most  approximation  methods,  is  that  the  limits  of  appli- 
cability or  degree  of  approximation  are  difficult  to  estimate. 

In  applying  equivalent  linearization,  one  writes  a set  of 
linear  differential  equations  to  describe  the  behavior  of  the 
real  nonlinear  system.  The  coefficients  of  the  linear  equations 


are  then  chosmin  such  a way  that  the  behavior  of  the  solution 
is,  in  some  way,  equivalent  to  the  behavior  of  the  real  system. 
To  illustrate  this  technique  let  us  consider 


I *+DU)+R(*)=K At)  , (ii-l) 

X V 

where  D($)  and  R($)  are  the  danpiig  and  restoring  moments  respectively. 
Let  us  further  restrict  the  problem  by  assuming  sinusoidal  ex- 
citation and  an  exact  damping  moment  which  varies  quadratically 
with  roll  velocity 

(II-2 ) 

The  equivalent  linear  damping  coefficient,  for  this  case 
has  been  derived  by  Blagoveshchensky  and  is 

. (II-3» 

where  is  the  amplitude,  and  w the  frequency  of  the  rolling 
motion.  This  expression  is  obtained  by  requiring  that  the  linear 
damping  results  in  the  same  dissipation  of  energy  per  cycle  of 
motion  as  the  quadratic  resisting  moment. 

In  the  case  of  random  seas,  the  concept  of  energy  interchange 
"per  cycle"  is  not  meaningful  since  the  motion  is  no  longer 
periodic.  Rather,  one  must  deal  with  time  averages,  but  the 
concept  of  equivalent  energy  interchange  still  applies.  Both 
St.  Denis  and  Vassilopoulos  estimated  equivalent  linear  coef- 
ficients for  roll  motion  in  random  beam  seas  by  the  following 
procedure.  Assume  that  the  actual  force  (moment)  which  is 

function  of  motion  variable  X.,  is  given  by  F.(X.).  This  is 

J «7 

to  be  approximated  by  an  equivalent  force  consisting  of  a linear 
term  in  the  appropriate  motion  variable  plus  an  error  term,  E 

F . (X  .)*C  .X  .+E.  (II-4) 

t j t j 

We  shall  require  that  C.  be  selected  such  that  the  time  average 

z< 
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of  the  mean  square  of  the  error  term,  E,  will  be  a minimum. 

(II-5) 


3<£2  > 

3C  . 
x 


= 0 


For  the  case  of  quadratic  damping  this  leads  to  the 
following  expression 


(II-6 ) 


Thus  it  is  seen  that  the  evaluation  of  ( I I — 5 ) leads  to  an 
equivalent  linear  force  coefficient  for  random  motion  which, 
as  in  the  case  of  sinusoidal  motion,  is  a function  of  the 
motion  amplitude  itself,  in  this  case,  a-  , the  RMS  value  of 
the  velocity . 

Applying  the  same  reasoning  to  a nonlinear  restoring  moment, 
let  us  represent  the  righting  moment  curve  by  a fifth-order,  odd 
polynomial  in  the  roll  angle  <J> , 

R($)=C  $+C  <t>3+C  ( II  — 7 ) 

1 3 5 


Designating  the  equivalent  linear  restoring  moment  by  an  error 
term  plus  cr<J>,  and  again  assuming  the  roll  motion  to  be  a zero 

t 

mean,  Gaussian  process,  having  an  RMS  amplitude  given  by  , 
minimization  of  the  error  term,  equation  (II-5)  leads  to 

C=C  +ZC  0?  +15C  a*  { 1 1 — 8 ) 

El  3 $ 5 


By  introducing  equivalent  linear  forces  and  moments  for  the 
"exact"  nonlinear  terms,  we  see  that  we  have  obtained  a linear 
system  in  which  these  forces  depend  on  the  mean  amplitude  of  the 
motion.  The  solution,  for  a given  amplitude  of  the  excitation, 
must  therefore  be  obtained  by  iteration.  That  is,  one  begins 
with  an  estimate  of  o^,  and  compu-r  the  moment  coefficients 
by  (II-6)  and  ( 1 1— 8 ) . Those  are  then  inserted  into  the  roll 
equation  of  motion.  . The  solution  for  roll  amplitude  is  then 


11 


obtained  and  its  RMS  value  compared  with  the  assumed  o^.  Con- 
vergence is  usually  quite  rapid. 

Now  consider  the  nature  of  the  coefficients  appearing  in 
(II-8)  for  a typical  ship.  Referring  to  FigureA-3  it  is  noted 
that  C ^ will  normally  be  positive  and  is,  in  fact,  the  product 
of  the  metacentric  height  and  displacement,  C will  be  positive 
corresponding  to  the  concave  upward  nature  of  the  initial  part 
of  the  curve,  and  C is  negative  corresponding  to  the  reversal 
of  the  curvature  beyond  approximately  30°.  Thus,  it  is  seen 
that  if  the  RMS  roll  is  in  the  medium  range,  the  equivalent 
linear  coefficient  will  be  greater  than  the  small  angle  linear 
value  (stiff  restoring  spring)  and  the  ship  will  have  maximum 
roll  response  to  wave  component  frequencies  higher  than  the 
small  angle  natural  frequency.  For  large  average  roll  ampli- 
tude however,  C predominates  resulting  in  a soft  equivalent 
effective  spring  constant,  and  the  maximum  response  is  to  wave 
frequencies  lower  than  the  small  angle  natural  frequency. 

This,  qualitatively,  describes  a situation  in  which  the 
consideration  of  certain  nonlinear  forces  in  the  equations  of 
motion  has  the  effect  of  modifying  the  results  that  might  have 
been  obtained  by  a linear  analysis,  but  without  introducing 
any  radical  changes  in  the  behavior. 

Let  us  now  consider  another  kind  of  nonlinear  term  which 
may  be  included  in  the  roll  equation  of  motion.  Recall  that  the 
higher  order  damping  and  restoring  terms  in  (II-2)  and  (II-7) 
may  be  considered  as  higher  order  derivatives  in  a Taylor  Series 
expansion  of  the  perturbation  force  about  a mean  state  of 
motion  of  the  ship.  In  such  an  expansion,  there  will  be  higher 
order  mixed  derivatives  as  well,  exemplified  by  a roll  re- 
storing moment  term  of  the  following  form 

R(t,y)=K^y.  (II-9) 

This  may  be  thought  of  as  expressing  the  change  in  the  linear 
restoring  moment  coefficient,  , with  changes  in  the  heave 
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position  of  the  ship  designated  by  y.  The  method  of  calculating 
, as  well  as  details  of  the  analysis  whose  results  are  given 
here,  are  found  in  a paper  by  Paulling  and  Rosenberg  (1959). 
Including  just  this  nonlinear  term,  the  roll  equation  of  motion 
( 1 1— 1 ) becomes 

I $4Cri+C  <t>y=K  (t).  (11-10) 

x L i < $y  e 


Let  us  now  assume  that  the  ship  is  proceeding  in  regular  head  or 
following  seas  so  that,  by  symmetry,  the  roll  excitation,  Kg(t) 
is  zero.  In  this  case,  the  heave  motion  is  sinusoidal  and  is 
given  by  y0  cos  ut.  Equation  (11-10)  then  becomes 

I $+C . f + $ ( C + X u sos  ut)=0  (II —11) 

x i*  l $:/*'  o 


For  the  time  being,  ignore  the  damping,  ?■  r . The  resulting  equation 
is  known  as  Mathieu's  equation  and  has  solutions  expressible  in  the 
form  of  special  functions.  It  should  be  noted  that  this  is  a linear 
differential  equation  in  $.  The  unique  feature  of  the  equation  is 
the  presence  of  a time-dependent  coefficient  of  the  roll  motion 
variable  $. 


The  solutions  to  Mathieu's  equation  have  a property  of  con- 
siderable importance  in  ship  rolling  problems  in  that  for  certain 
values  of  the  frequency,  u>,  the  solution  is  unstable.  Physically, 
this  implies  that  if  the  roll  motion  described  by  (-II-11)  is  taking 
place  in  an  unstable  region,  the  amplitude  will  grow  until  limited 
by  a physical  constraint  not  included  in  this  equation.  The  unstable 
frequency  may  be  found  as  follows. 

Introducing  the  change  of  variable  t =u>t , equation  (11-11) 

(with  CT  neglected)  may  be  rewritten 


where 


(II-12) 
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£ = kJjL, 

I U)2 
X 

(d  */C  /I  * natural  frequency  in  roll. 
rt  i x 

Figure  II-2  is  a graph  of  e versus  6 for  the  Mathieu  equation  in 
which  regions  corresponding  to  stable  solutions  have  been  shaded. 
It  is  seen  that  for  small  values  of  c,  corresponding  to  small 
heave  motion,  yot  we  have  unstable  solutions  of  Mathieu' s equation 
and  therefore,  unstable  rolling  motion  occurring  for 

6=1/4,  4/4,  9/4,.  . .=(n/2)2;n=l,2,3.  . . 
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FIG. 1-2  STABLE  AND  UNSTABLE  REGIONS 
FOR  THE  MATHIEU  EQUATION 
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Since  6 the  square  of  the  ratio  of  the  nat  ral  frequency  to 
exciting  frequency#  this  means  that  unstable  roll  occurs  when 
this  frequency  ratio  takes  on  half-integer  values.  The  lowest 
and  widest  unstable  region  occurs  at  an  exciting  frequency  of 
twice  the  natural  frequency.  Thus#  unstable  roll  will  be  ex- 
cited if  the  ship  encounters  head  or  following  seas  with  a 
frequency  of  encounter  equal  to  twice  the  natural  frequency  of 
roll.  For  usual  ship-wave  proportions  in  which  capsizing  is 
likely,  this  occurs  only  in  following  seas.  Note  also  that 
the  instability  may  occur  at  an  encounter  frequency  equal  to 
the  natural  frequency  of  roll#  2/3  times  the  natural  frequency 
and  so  on.  The  width  of  the  unstable  region  for  these  nigher 
frequencies  is  less, however , than  that  at  the  lowest  frequency 
with  the  practical  resu)  t that  for  unstable  motion  to  be  ex- 
cited at  the  higher  frequency,  the  amplitude  of  the  variable 
term#  e,  and  therefore  yof  must  be  greater.  For  practical 
purposes,  only  the  two  lower  unstable  frequencies  have  any 
significance. 

The  existance  of  this  parametrically  induced  unstable 
motion  has  been  investigated  theoretically  and  experimentally 
by  Grim  (1952),  Kerwin  (1955),  and  Paulling  and  Rosenberg  (1959). 
The  former  two  utilized  a ship  model  equipped  with  a pair  of 
weights  which  rotated  at  constant  angular  velocity  in  the  ver- 
tical plane.  In  this  way,  the  vertical  position  of  the  center 
of  gravity,  thus  the  metacentric  height  was  forced  to  vary 
sinusoidally.  The  latter  performed  a constrained  model  ex- 
periment in  which  the  ship  model  was  mechanical) y forced  to 
oscillate  vertically  in  heave  but  with  mechanical  attachment  to 
the  oscillating  mechanism  which  allowed  freedom  in  roll.  Thus 
the  condition  of  sinusoidal  heave  was  induced  artificially  in 
calm  water.  In  both  types  of  experiments,  an  unstable  rolling 
motion  was  found  to  be  set  up  at  frequencies  predicted  by  the 
Mathieu  instability  chart.  Grim  showed  that  the  motion  would 
ultimately  reach  a steady  state  value  determined  by  nonlinear 
damping,  Kerwin,  in  addition  to  his  experiments,  carried  out 
a numerical  integration  of  the  equation  of  roll  motion.  He 
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concluded  from  this  solution  that  the  transient  behavior  in  the 
region  of  instability  was  quite  sluggish  and  that  a large  number 
of  oscillations  would  be  experienced  before  the  amplitude  of 
roll,  starting  from  a small  arbitrary  disturbance,  would  grow  to 
appreciable  magnitude.  This  conclusion  is  not  in  accord  with 
observations  of  the  model  behavior  during  some  of  the  experi- 
ments which  form  the  subject  of  later  parts  of  the  present 
paper  and  we  shall  return  to  this  point. 

The  above  discussions  has  described  the  phenomenon  of 
parametrically  induced  unstable  roll  motion.  The  excitation 
consists  of  the  periodic  variations  in  static  transverse 
stability  which  result  from  the  variation  in  geometry  of  the 
immersed  hull  during  periodic  heave  oscillation  of  the  ship. 

When  the  ship  moves  through  head  or  following  seas,  the  immersed 
geometry  variations  result  from  the  changing  water  profile 
along  the  ship  or  wave  surface,  as  well  as  from  the  ship  motion. 
The  relationship  between  transverse  stability  and  the  wave 
proportions  in  following  seas  has  been  investigated  by  Paulling 
(1961)  . It  is  shown  that  significant  variations  in  the  righting 
arm  curve  are  caused  by  the  changes  in  the  immersed  geometry 
as  the  ship  is  overtaken  by  successive  crests  and  troughs  in 
following  or  head  seas.  In  general,  for  waves  of  length  nearly 
equal  to  the  ship  length  the  ordinates  of  the  righting  arm 
curve  are  increased  if  a wave  trough  is  near  amidships  and 
decreased  when  a crest  is  near  amidships.  The  effect  is  most 
pronounced  in  steep  waves  and  for  ships  of  low  freeboard. 

Typical  curves  of  righting  arms  for  these  two  conditions,  crest 
amidships,  and  trough  amidships  are  shown  in  Figure  A-4. 

If  the  entire  righting  arm  curve  shown  in  Figure  A-4  rather 
than  just  the  initial  slope  is  considered  in  a Mathieu  type 
analysis,  one  obtains  an  equation  of  the  form 

1 xb+cI}+RM($,t)=° . (II-13) 

Here,  HM(<$>,t)  represents  the  time-dependent  variation  between 
the  extreme  righting  arm  curves,  trough  amidships  to  crest 
amidships,  as  waves  move  past  the  ship.  Neither  a solution  nor 
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a motion  stability  chart  is  available  for  (11-13) . Some  char- 
acteristics of  the  large-amplitude  motion  under  the  influence 
of  typical  ship  righting  moment  curves  have  been  obtained  by 
numerically  integrating  the  following  differential  equation 

I 'i+CTi+  C (<b)+c  (<t>)co8  u t=0.  (11-14) 

x l o i e 

In  this  equation  is  the  mean  of  the  trough  amidship  and 

crest  amidships  righting  moment  curves,  and  C^($)  is  one-half 
their  difference.  Two  important  features  of  the  resulting 
motion  were  obtained  in  these  studies: 

(1)  Motion  instabilities  were  obtained,  just  as  in  the 
case  of  the  linear  Mathieu  equation,  resulting  in  a growth  of 
motion  from  an  initial  small  disturbance.  As  the  motion  grew 
however,  the  "effective  stiffness"  of  the  restoring  moment 
(see  Equation  II-8)  changed  from  its  small-angle  value.  The 
system  therefore  became  "detuned"  from  the  exciting  frequency 
and  the  motion  would  tend  to  limit  itself.  With  a different 
value  for  the  initial  disturbance,  the  instability  would  be 
found  at  a somewhat  different  frequency  corresponding  to  the 
altered  effective  stiffness  associated  with  this  amplitude. 

(2)  In  Figure  A-4  it  is  observed  that  the  variations  in 
stability  with  wave  position  are  not  symmetrically  arranged 
about  the  still  water  curve.  The  crest  amidships  curve  nearly 
coincides  with  the  still  water  curve  up  to  about  twenty  degrees, 
while  the  trough  amidships  curve  is  significantly  higher  than 
the  other  two.  This  asymmetry  is  less  pronounced  in  ships  with 
low  freeboard  amidships  and  relatively  little  flare  in  the  ends. 

The  effect  of  the  asyiu.ietrical  nature  of  the  righting  arm 
curves  is  to  make  the  ship  respond  to  parametrically  excited 
motion  as  though  the  mean  restoring  moment  initial  slope  were 
greater  than  the  still  water  value.  Thus,  for  the  case  illustrated 
in  Figure  A-4#  the  first  Mathieu  unstable  frequency  will  be 
somewhat  higher  than  that  predicted  by  using. the  still  water 
metacentric  height.  Further,  since  the  stability  variation 
between  trough  and  crest  varies  with  wave  height,  this  unstable 
frequency  will  also  be  different  for  different  wave  heights. 


17 


This  departure  of  the  Mathieu  instability  frequency  is 
thought  to  explain  certain  important  differences  in  the  ob- 
served motion  response  of  the  two  different  ship  models  which 
were  tested.  It  will,  therefore,  be  referred  to  again  in  a 
later  chapter  in  which  the  experimental  results  are  discussed. 

If  the  ship  is  moving  through  a following  sea,  it  is 

apparent  that  the  ship  will  experience  periodic 
variations  in  its  transverse  stability,  and  these  variations  will 
affect  the  roll  motion  of  the  ship.  Although  the  preceeding 
discussion  was  concerned  with  regular  waves,  a somewhat  similiar 
phenomenon  but  perhaps  of  a transient  nature  might  be  expected 
if  the  ship  encounters  a wave  group  having  sufficient  regularity 
and  steepness.  Such  behavior  has,  in  fact,  been  observed 
during  some  experiments  to  be  described  in  the  next  section  and 
the  nature  of  this  observed  motion  suggests  that  the  wave  in- 
duced stability  changes  and  resulting  parametrically  induced 
roll  play  an  important  role  in  extreme  motions  and  capsizing. 

The  examples  just  described  illustrate  the  point  that 
contemporary  linear  and  nonlinear  ship  motion  theory  can  reveal 
many  features  of  ship  behavior  of  interest  to  a study  of 
capsizing,  but  it  falls  short  of  describing  all  aspects  of  the 
problem.  This  suggests  that  experiments,  conducted  in  conditons 
which  are  as  realistic  as  possible,  offer  the  best  possibility 
for  an  initial  study  of  capsizing.  The  primary  objective  will 
be  to  gain  insight  and  understanding  of  the  basic  phenomena.  A 
second  purpose  will  be  to  gather  some  data  upon  which  to  base 
conclusions  and  to  be  used  as  calibration  for  theoretical 
studies  of  the  motion. 

As  noted  in  the  INTRODUCTION,  experiments  utilizing  radio- 
controlled  free  running  ship  models  have  been  conducted  several 
times  in  the  past  to  study  aspects  of  ship  behavior  and  control 
for  which  adequate  theory  is  lacking  or  for  which  conventional 
towing  tanks  are  inadequate.  In  order  to  conduct  free  running 
experiments,  one  requires,  in  addition  to  a suitable  test  site, 
certain  hardware  items. 
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These  include  the  test  model  itself,  complete  with  radio 
controls,  scaled  autopilot  steering  systems,  self-contained 
propulsion  and  monitoring  instrumentation  to  record  the  per- 
tinent motion  variables,  apparatus  for  recording  the  wave 
conditions  in  the  test  area,  and  some  kind  of  floating  plat- 
form from  which  to  conduct  the  experiment.  Such  experiments 
performed  in  San  Francisco  Bay  during  the  past  three  years 
are  described  in  the  following  section. 


! 
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EXPERIMENTS 


A . Sea  Spectra  Measure ment  & Analysis 

If  the  model  tests  are  to  reflect  realistic  extreme  sea 
situations,  it  is  obviously  important  that  the  waves  at  model 
scale  adequately  represent  full  scale  conditions.  An  important 
part  of  the  program  has  therefore  been  the  recording  and  analysis 
of  the  seaway.  The  following  sections  provide  a brief  analysis 
of  the  measured  wave  data  and  a qualitative  comparison  with  full 
scale  storm  conditions.  Appendicies  D and  E contain  detailed 
tabulations  of  wave  statistics  and  plots  of  both  point  and  direc 
tional  sea  spectra.  Chapter  VI  contains  a discussion  of  wave 
encounter  and  wave  group  phenomena. 

1.  Test  Site,  Array,  and  Procedures 

The  test  area,  known  as  the  "Berkeley  Flats" » is  a region  of 
gradually  shoalling  water  with  an  average  depth  of  from  eight  to 
twelve  feet  and  a mud  bottom.  Tidal  currents  are  generally 
under  one  knot  and  the  test  area  is  outside  of  the  regular 
commercial  traffic  lanes. 

During  the  summer  months,  the  overall  wind  pattern  is 
extremely  stable  and  repeats  itself  day  after  day.  The  hot  air 
rising  in  the  central  valley  draws  the  cooler  ocean  air  through 
the  Golden  Gate.  The  test  site  is  located  down  wind  and  slightly 
off  to  one  side.  The  fetch  is  between  three  and  six  nautical 
miles  in  the  form  of  a diverging  funnel,  two  miles  wide  at  the 
Gate  and  over  four  just  before  the  test  area.  The  overall 
jeometry  of  the  area  is  shown  in  Figure  III-l.  On  a typical  test 
day,  the  winds  would  be  light  in  the  early  morning  and  build 
gradually  to  about  fifteen  knots,  depending  on  location,  by  mid 
afternoon.  The  waves  are  a typical  "wind  sea."  Ocean  swell  is 
rarely  noticeable  in  the  Bay  and  is  of  no  consequence  in  the 
test  area. 

The  waves  were  measured  using  an  array  of  step  type  wave 
gages  (Unddara,  1971  and  Chou  et  al.,  1972  ) . Three  of  the  gaqes 
were  mounted  on  ants  extending  from  a wave  buoy  and  the  fourth 
gage  was  located  at  the  center.  The  wave  buoy  is  of  the  tension 
ion  variety,  each  interconnected  column  being  held  down  with 
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its  own  anchor  line.  The  buoy  is  very  stable  and  has  virtually 
no  motion  in  heave  and  probably  less  than  ten  degrees  in  yaw 
in  extreme  conditions.  Since  the  support  vessel  supplied 
electrical  power  and  carried  the  tape  recorder,  it  was  necessary 
to  anchor  and  connect  with  the  buoy  via  a long  electrical  cable 
on  the  bottom.  The  wave  measurements  were  monitored  with  a 
strip  chart  recorder  and  the  observed  wind  and  sea  conditions 
recorded.  The  analog  records  were  later  digitized  at  0.2 
second  intervals. 

2.  Point  Spectra  and  Statistics 

The  waves  were  recorded  for  fifteen  minutes  approximately 
every  hour  to  hour  and  one-haif  or  when  there  was  a significant 
change  in  conditions.  Observations  of  wind,  tide,  currents, 
wave  height  and  direction,  and  other  potentially  useful  data  were 
noted  on  the  data  sheets.  Statistics  were  computed  from  the  digitized  time 
series  and  derived  from  the  spectrum  computations.  Appendicies 
D and  E contain  the  summaries  of  the  rur.s  and  statistical 
information  for  both  the  1972  and  1973  test  seasons.  The  1973 


Figure  III-l  Test  Site  and  San  Francisco  Bay  Area 


overview  tabulation.  Table  E-  1,  contains  a summary  of  the  notes 
taken  during  the  tests.  Figure  III-2  shows  the  distribution  of 
wave  recordings  with  respect  to  significant  wave  height  (4*RMS) 
for  both  1972  and  1973. 

The  spectral  analysis  was  conducted  using  a fast  Fourier 
-nalysis  technique  with  2048  points  and  a 0.4  second  time  interval. 
Smooti.ing  was  done  by  block  averaging  over  eight  or  sixteen 
frequencies . 

An  outline  of  the  method  is  given  in  Appendix  H of  Haddara,  et  a l 
(1971).  The  point  spectra  of  one  probe  recorded  on  each  day  are 
shown  in  Figure  E-2.  The  growth  of  the  spectra  appears  to  follow 
the  general  description  of  a fetch  limited  sea  as  reported  by 
JONSWAP,  The  Joint  North  Sea  Wave  Project  (Hasselman,  et  al,  1973)  . 
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Figure  1 1 1-2  Distribution  of  Significant 
Wave  Height,  1972  and  1973. 
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A parameterized  one-dimensional  spectrum  was  used  to  fit  that 
data.  The  spectrum  was  given  as  the  product  of  a Pierson- 
Moskowitz  (1964)  type  spectrum  with  parameters  a and  / : 

SF,y(j')=OLc;2(2T\)~uj~bexp  l- 1 . 25  ( f^/f)  * ] , (III-l) 

and  a "peak  enhancement"  factor: 


Y exp 


(III-2) 


where 


c = 


C o for  f < f 
la  J J r 


c;.  for  > fm 


width  of  the  spectral  peak 


= frequency  of  the  peak 


rq c p • f 

Y = ;/s‘  (fj 


Typical  values  of  the  constants  can  be  estimated  from  the  non- 
dimensional  fetch 


x = xg/U 
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where  x is  fetch  and  £/,_  is  the  wind  velocity  at  ten  meters, 

x U 


using 


f = 9 


n 


3.5  x 


0*3  3 


, a = 0.076  x 


0*22 


(III-3) 
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the  (JONSWAP)  estimates  of  y,  and  showed  much  scatter 

but  were  assumed  to  be  approximately  3.3,  0.07,  and  0.05, 
i.espectively . The  fetch  derived  from  Equation  III-3,  given  the 
joectral  peak  f and  the  measured  wind  velocity  U - U /l. 11, 

w 5 10 

i jnqed  from  approximately  one  nautical  mile  at  low  wind  speeds 
to  about  seven  nautical  miles  at  the  higher  wind  speeds.  The 
scatter  in  the  present  data  is  due  in  part  to  inaccuracies  in 
wind  speed  measurements  since  these  were  not  taken  over  the  full 
fifteen  minute  wave  recording  period.  However,  the  derived  fetch 
appears  tc  correlate  with  the  direction  of  the  current  in  the 
generating  area. 
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Figure  III-3  Wind  Speed  vs  Significant  Wave  Height 

The  significant  wave  height  versus  wind  speed  is  shown  in 
Figure  HI-3.  As  expected,  an  ebb  current  produces  larger  waves 
since  it  opposes  the  mean  • ind  direction.  The  anomalous  low 
wind  speed  and  high  significant  wave  heights  are  usually 
associated  with  decaying  spectra  and  cross  sea  situations. 

The  peak  ratio  y,  unlike  JONSWAP,  varied  smoothly  from 
approximately  3.3  for  low  spectra  to  0.98  for  the  highest.  The 
fit  of  the  parameterized  model  with  the  experimental  spectrum 

is  shown  in  Figure  III- 4.  The  parameterized  model  was  computed 

OBS 

using  the  observed  peak  frequency  and  amplitude  S ( fn) 

and  appears  to  agree  reasonably  well.  It  may  be  possible  to  use 
such  models  as  aids  in  interpolating  for  spectra  at  times  other 
them  those  measured. 

On  one  occasion,  nine  wave  recordings  were  made  during  a 
four  hour  period  in  order  to  investigate  wave  stationarity . 
Consecutive  records  yield  similar  spectra  with  amplitude  and 
directional  properties  that  correlate  very  well  with  the  observed 
values.  Futhermore,  the  changes  appeared  to  be  smooth  functions 
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SEA  WAVE  SPECTRA  OF  SAN  FRANCISCO  BAY 


of  time.  If  there  were  no  unusual  changes  in  the  observed  wind 
and  wave  characteristics,  it  is  therefore  probably  reasonable 
to  interpolate  between  consecutive  measurements  for  an  estimate 
of  intermediate  sea  conditions. 

All  of  the  spectra  peaked  at  f = 0.28  Hz  or  above  with  one 
exception  which  peaks  at  approximately  0.22  Hz.  The  large 
spectrum  on  that  day  is  the  only  one  with  any  energy  at  f - 0.2  Hz, 
which  corresponds  to  a deep  water  wave-length  of  56.9  feet.  Even 
at  this  frequency  there  is  only  a nine  percent  decrease  in  the 
wavelength  due  to  shallow  water.  Since  this  is  still  3.1  times 
longer  than  the  models  being  tested,  the  effect  of  the  shallow- 
water  is  not  considered  to  be  significant. 

Some  spectra  exhibited  a small  spike  below  f = 0.02  Hz 
peaking  at  zero  frequency.  These  correlate  wen  with  rapid  changes 
in  the  tide  at  the  time  of  recording. 
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The  significant  wave  height  computed  from  the  point 
spectrum  differs  considerably  from  the  analysis  of  the  time 
series.  The  spectral  statistics  agree  with  those  of  time  series, 
however,  if  any  zero  frequency  spike  is  removed  from  the 
spectral  area  and  the  significant  wave  height  is  corrected 
using  the  broadness  facter  (Loukakis,  1970  ). 

H i/  - Hi.  *sfl-e2/z 

/ 3 ~ / 3 

time  series  spectrum 

where 

C2  * l-m2/(m  m ) 

2 0 4 

00 

m - f wnS(u)  ddi. 


As  noted  earlier,  the  measured  spectra  do  not  appear  to 
be  significantly  distorted  due  to  the  limited  fetch  or  water 
depth.  The  amplitude  and  frequency  content  also  seem  to  scale 
reasonably  well.  Figure  III-5  is  a plot  of  a number  of  scaled 
San  Francisco  Bay  point  spectra  assuming  a scale  ratio  of 
\=30.  Figure  III-6  is  a collection  of  full  scale  storm  spectra 
recorded  at  a number  of  locations.  These  represent  a variety 
of  shapes  and,  it  might  be  noted  are  not  the  largest  available. 
The  model  and  full  scale  spectra  are  seen  to  agree  quite  well. 
These  are  admittedly  rathor  severe  conditions  that  occur  only 
rarely,  but  they  probably  represent  the  very  situations  where 
capsizing  becomes  an  important  danger.  Figure  III-7  gives 
the  equivalent  Pierson-Moskowitz  spectra  for  comparison  purposes. 
Our  knowledge  of  the  shape  of  storm  sea  spectra,  let  alone  the 
detailed  time  history  of  the  wave  elevation,  is  woefully 
lacking.  The  detailed  experimental  information  provided  in  the 
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Figure  III-5  Typical  San  Francisco  Bay  Spectra  Scaled  using  A=30 


FULL  SCALE  STORM  SPECTRA 


[1]  Wave  Spectra  estimated  from 
355  North  Pacific  Wave  Record 
Station  Papa 

[2J  Ocean  Weather  Station  A, 

62*N  32.2*W,  12  Jan.  1973 

[3]  Data  taken  during  SkyLab-FLT 
730606,  Hurricane  Ava. 
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Figure  II1-6  Full  Scale  Storm  Spectra 
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PIERSON-MOSKOWITZ  SPECTRA 


*SJ 


Figure  III-7  Pierson  Moskowitz  Spectra 

appendix  cannot  be  adequately  judged  since  comparable  full  scale 
data  is  just  not  available.  On  the  other  hand,  this  may  not  be 
the  case  in  the  near  future.  Long  range  monitoring  of  the 
sea  surface  using  radar  reflections  off  of  the  ionosphere  and 
other  sophisticated  measuring  techniques  are  currently  under 
development.  These  show  promise  of  providing  sufficient  data 
to  evaluate  quantitatively  the  adequacy  of  these  experimental 
spectra  and  the  often  employed  "Standard"  sea  spectra. 

3.  Directional  Spectra 

The  major  reason  for  experimenting  in  the  open  waters  of 
San  Francisco  Bay  was  the  desire  to  investigate  ship  capsizing 
in  more  or  less  "reasonable  extreme"  sea  conditions.  This  of 
course  implies  the  possibility  that  the  sea  spectrum  ould  be 
multi-directional.  A measure  of  the  directionality  was  therefore 
considered  vital  for  the  interpretation  of  the  resulting  motions 
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A number  of  analysis  techniques  were  evaluated.  The  first 
method  involved  computing  the  first  few  terms  of  a Fourier  series 
expansion  of  the  directional  spreading  function  at  each  frequency 
(Oakley,  1973) . This  method  has  not  proven  to  be  well  suited 
to  the  analysis  of  a small  number  of  gages.  Spurious  peaks  and 
negative  spectral  estimates  are  frequent  problems.  The  peak 
direction  can  be  accurately  computed  but  not  the  amount  of 
spreading.  Knowing  th^  peak  direction,  0O  „a  ooen(Q-Qa)  spreading 
function  could  be  assumed.  The  current  literature  on  the  subject 
of  directional  spreading  advocates  a cosine  squared  spreading 
function  within  the  generating  area.  A narrower  spreading  relation 
is  more  appropriate  the  farther  one  moves  from  the  region  of  wave 
generation.  The  test  area  appears  to  be  within  the  wind  wave 
generating  region  so  a moderately  broad  spreading  function  should  be 
applicable.  However,  this  does  not  appear  to  be  the  case  from 
observations.  The  waves  appear  to  be  quite  unidirectional  and 
relatively  narrow  banded  in  character.  Another  analysis  technique 
investigated  was  the  maximum  likelihood  filter  developed  by 
Capon  (1969).  This  procedure  requires  only  the  inverse  of  the 
cross-spectral  matrix  and  is  therefore  rather  simple  to  use.  The 
wave  number- frequency  spectrum  estimate  is  given  by 

K — j 

S(k,f)  - E I Q.1(f)exp(-ik*X.1)) 
j,l^l  JL  - 

where 

A 

[Q.^(f)]  = inverse  of  the  cross  spectral  matrix  with 

3 elements  corresponding  to  the  gage  pair  3-1. 

k * wave  number  (directional)  for  water  of  depth  h. 

X..  = distance  between  gage  pair  j-l  in  an  array  of 

~3  K gages. 

This  method  offers  significantly  improved  resolution  over 
the  Fourier  series  method  and  yields  positive  definite  spectral 
estimates.  Cross  plots  of  the  directional  spectrum  S(Q,f)  for 
September  17,  1973,  run  number  3,  at  selected  frequencies  are 
shown  in  Figure  III-3.  The  directions  are  shown  relative  to 
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the  buoy.  The  peak  direction  and  relative  nanowne.';  agi  >c  ’i  . • 
well  with  observation.  The  narrowest  peak  corresponds  to  n 
cosine  spreading  function  raised  to  the  6.3  power.  Resol uti 
appears  to  be  acceptable  for  wavelengths  up  t<  an  order  01 
magnitude  longer  than  the  average  distance  between  gage  pai 
Note  however  that  with  only  four  gages,  resolution  is  not  w 
The  filter  is  only  able  to  achieve  complete  garntum  of  a 
of  sine  waves  if  they  are  separated  by  at  15  oi  more 

degrees  in  azimuth.  Since  the  spectral  estimate  is  given  L *. 
inverse  of  a trigonometric  series,  theii  is  a tendency  to  si  rr  u 
energy  into  all  directions.  The  amount,  of  spreading  inci>  a., 
with  a decrease  in  frequency  or  the  number  of  gages.  Sprea 
also  increases  with  the  amplitude  of  the  peak  of  the  distrii  1 
with  such  a small  array.  Care  must  therefore  bo  taken  in 
applying  the  results.  Small  amounts  of  spun,  us  energy,  < 
to  the  peak,  going  in  the  wrong  direction  nay  shew  up  in  th< 
prediction  of  ship  accelerations,  for  example 

Inproved  resolution  is  possible  using  a oneji  plug:  air 
approach  (Oakley,  1973).  This  technique  attempts  to  minimii 
the  largest  deviation  between  the  measured  cross  spectrum  anu 
a parameterized  estimate.  The  method  capable  of  achievi  ig 
a very  high  degree  of  resolution,  even  fci  a limited  numbei  of 
gages,  by  looking  only  at  those  directions  containing  energy. 

If  pushed  too  far,  however,  the  directional  distr  Lbut  ion 
estimate  begins  to  resemble  a distribution  ot  spikes  even  ± £ th 
actual  distribution  is  smooth.  Some  care  in  interpreting  t,he 
results  is  therefore  in  order.  If  the  spectrum  is  known  to  se 
"smooth",  a smoother  parameterization  than  the  tent  functions 
used  in  Oakley  (1973)  should  be  consideied.  The  linear  program- 
ming technique  was  used  to  check  the  lesults  usmy  Capon's 
analysis.  The  results  support  the  previous  observation  cha 
the  directional  distributions  tend  to  be  very  narrow,  especially 
for  the  higher  seastates  and  lower  frequencies. 

The  point  spectra  and  directional  spent. ra  for  the  lg7J  test 
season  are  plotted  in  Appendix  E.  The  directional  spectra 
are  presented  as  contour  plots  relative  to  magnetic  North. 

These  are  somewhat  easier  to  interpret  than  tin  cluttered  cross- 
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plots  of  spectral  density  versus  direction.  Note  however  that 
the  amount  of  directional  spreading  is  deceptive.  The  contour 
plots  obscure  the  relative  narrowness  of  the  peaks  of  the 
larger  spectra.  Note  also  that  the  contouring  levels  differ 
from  run  to  run. 

Hun  4 on  August  24,  1973  offers  a particularly  dramatic 
example  of  the  power  of  the  directional  analysis.  The  winds  were 
rather  light  and  erratic  during  the  morning  and  early  afternoon 
at  the  initial  test  location.  Potentially  better  westerly 
wind  and  sea  conditions  were  apparently  developing  in  the  test 
area  closer  to  the  Berkeley  pier.  Upon  arrival  at  this  new 
location  a West  or  Northwesterly  breeze  had  also  developed  to 
the  north  of  tne  previous  test  site.  The  resulting  seaway  was  a 
classic  example  of  two  wave  trains  crossing  at  right  angles, 
both  of  approximately  the  same  amplitude  and  wave  length.  The 
two  spectral  peaks  are  clearly  shown  in  the  contour  plot.  The 
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connecting  lines.  Figure  F-3.27  are  spurious  directional  spreading 
at  constant  frequency. 

The  wave  buoy  direction  was  not  recorded  on  three  occasions. 
Since  the  wind  direction  and  spectral  peak  were  always  in  good 
agreement,  the  buoy  orientation  was  assumed  to  be  the  value 
required  to  bring  them  into  line.  This  fact  is  noted  on  the 
contour  plots.  The  wave  array  and  model  headings  are  probably 
reliable  only  to  within  + 10  degrees  due  to  reading  errors  and 
variable  amounts  of  deviation. 

An  additional  seaway  topic  of  importance  is  that  of  wave 
groups.  A discussion  of  the  theoretical  and  experimental 
statistics  is  presented  in  Chapter  VI. 

B.  Model  Experiments 

Extensive  experiments  have  been  conducted  in  the  open 
waters  of  San  Francisco  Bay  for  over  three  seasons.  Two 
seventeen  foot  models  were  tested.  The  first  was  a medium  speed, 
full  hull  form  (\=20,  American  Challenger)  of  wood  construction 
(Haddara,  et  al,  1971).  The  second  model  represented  -\  high 
speed,  narrow  hull  form  {>.=££,  Sealand-7  Container  Clare), 

(Chou,  et^  al^,  19  73)  and  was  the  more  fully  instrumented.  Body 
plans  and  principal  dimensions  of  the  two  models  are  given 
in  Appendix  A.  Comparisons  given  in  this  report  are  based 
upon  the  same  assumed  scale  ratio,  \=Z0  for  both. 

The  second  model  is  a fiberglass  reinforced  epoxy  sandwich 
construction  with  a balsa  core.  Large  hatches  with  O-ring  seals 
were  provided  for  easy  access  and  proved  to  be  completely  \.  • t*»r- 
tight.  Table  A-I  lists  the  particulars  of  the  twe  models. 

' 

These  models  were  self-propelled,  using  a 24  volt  DC  PM  field 
motor,  and  were  radio  controlled.  The  rechargeable  bettory 
circuitry  is  divided  to  provide  power  to  the  instruments  end 
propulsion  separately.  Four  batteries  are  also  used  as  movable 
ballast  for  GM  changes  in  the  second  model. 

An  analog  tape  recorder  in  the  mode1  was  used  to  record 
four,  and  later  seven,  channels  of  data  selected  from  the 
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tol lowing: 

a. )  pitch,  roll,  yaw,  and  yaw  rate  angles  (gyros) 

b,  ) heave,  sway  and  surge  accelerations  for  a location  near 

the  CG  of  the  model. 

e.)  motor  RPM,  relative  water  speed,  and  rudder  angle. 

Tlu  autopilot  system  was  designed  to  simulate  a typical  full 
;c.ile  system  with  heading  error  and  yaw  rate  feedback,  a dead  band 

i:  * lag.  A detailed  description  is  presented  in  Haddara,  et  al. 
( nd  an  analysis  of  the  actual  system  characteristics  is 

in  Appendix  C of  this  report. 

The  righting  arm  curves  for  the  two  models  in  calm  water 
no  in  waves  are  given  in  Figures  A-3  and  A-4. 

After  recording  the  waves,  the  wave  buoy  cable  was  dis- 
connected and  left  with  the  anchored  dinghy.  The  model  was  then 
launched  such  that  the  wave  buoy  would  be  passed  close  aboard 
■v  approximately  half  way  through  the  run.  The  compass  heading 
the  model  was  recorded  but  the  heading  characterizations, 
j .■  - following,  quartering,  or  beam,  were  selected  strictly  by 
eye.  During  the  course  of  the  run  the  model  would  occasionally 
deviate  from  the  desired  heading  due  to  initial  heading  errors 
tiu'iny  uncaqina,  drift  of  the  yaw  gyro,  or  a change  in  the 
bs«?vod  principal  direction  of  the  seas.  Course  corrections 
r<  radioed  to  fhe  mode’  in  the  form  of  two  degree  steps  in  the 
• w gyro  heading  angle.  If  the  model  was  able  to  maintain 
> use,  the  relative  heading  was  held  substantially  constant, 
inee  the  sea  conditions  were  often  extreme,  the  model  would 
casionally  experience  a broach  or  a series  of  large  yaw  forces 
nd  would  be  thrown  off  course  into  a large  quartering  or  even 
b*  am  sea  heading.  At  high  speeds  there  was  sufficient  directional 
tto]  for  the  rudder  (generally  hard  over)  to  bring  the  model 
1 ewly  back  on  course  as  soon  as  the  large  group  of  waves  had 
•jsed.  At  the  lowest  Froude  number  and  largest  seas,  the  model 

number  of  occasions  was  unable  to  recover  for  long  periods 
wallowed  in  beam  seas  until  a long  lull  arrived.  The  spectral 
Jysis  of  these  runs  generally  show  very  large  values  at  low 
' r jenci.es  in  the  ya’> 


and  rudder  angles. 
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If  a capsize  occurred  before  the  end  of  a fifteen  minute 
run,  the  model  was  recovered,  the  yaw  gyro  re- caged,  and  the 
procedure  repeated  until  there  was  a noticeable  change  in  sea 
conditions  or  approximately  one  hour  had  elasped  since  the  last 
wave  measurements.  Launching  and  recovery  in  high  seas  was  a 
tricky  operation  requiring  both  skill  and  time.  Two  long  runs 
or  three  to  four  short  runs  were  generally  all  that  could  be 
done  between  wave  measurements.  Nevertheless,  a substantial 
number  of  experiments  were  conducted  in  a variety  of  sea  conditions, 
speeds,  GM's,  and  headings. 

Tables  D-V  and  E-III  provide  an  overview  and  summary  of  the 
model  runs  during  the  1972  and  1973  test  season.  The  statistics 
were  derived  from  an  analysis  of  the  time  series  rather  than 
from  the  spectrum  of  motions. 

At  the  completion  of  the  day's  testing  the  onboard  model 
tape  recorder  was  removed  and  the  analog  signal  inspected  on  a 
strip  chart  recorder.  The  data  were  digitized  each  week  and  the 
results  preserved  in  both  analog  and  digital  form. 

Many  runs  were  too  short  to  make  either  a statistical  or 
spectral  analysis  meaningful.  Only  runs  over  seven  minutes  have 
been  selected  for  statistical  analysis,  the  majority  being  over 
ten  minutes.  Many  of  these  involve  rather  violent  motions, 
especially  in  roll,  yaw,  and  rudder  deflections.  These  are 
necessary  conditions  for  a capsize  to  be  likely,  but  this  makes 
a comparison  with  a linear  theory  analysis  rather  difficult. 

The  mean  relative  speed  between  the  model  and  water  was 
measured  on  numerous  runs  using  a paddle  wheel  type  speedometer. 
Figure  III-9  is  a plot  of  the  mean  speed  versus  significant  wave 
height  (interpolated  to  the  time  of  the  run)  for  the  various 
headings.  The  measured  speed  is  seen  to  decrease  with  increasing 
significant  wave  height,  even  in  following  seas.  If  there  were 
significant  surge  motions,  the  speedometer  would  be  expected  to 
read  a lower  value  due  to  the  decreased  relative  velocity  between 
the  model  and  water  velocity  on  the  face  of  the  wave.  Computations 
were  made  using  the  linear  motions  program,- using  both  the  calm 
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water  speeds  and  the  lower  speeds  reported  by  the  speedometer. 

Due  to  the  nature  of  the  folding  of  the  motion  spectra  for  con- 
version to  the  encounter  domain,  the  frequency  distribution  is 
quite  sensitive  to  speed  changes.  Better  frequency  agreement 
was  found  for  speeds  close  to  the  speedometer  readings,  rather 
than  the  calm  water  speeds.  This  would  suggest  that  the  mean 
velocity  did  indeed  decrease  and  that  in  the  mean  surging  was 
not  too  significant.  This  agrees  with  the  observations  made 
in  the  Bay.  An  estimate  of  surging  is  possible  by  using  the 
surge  accelerometer  data  taken  on  about  a third  of  the  runs. 

The  RMS  surge  motion  and  RMS  surge  velocity  were  obtained  by 
integrating  inverse  moments  of  the  acceleration  spectrum.  Typical 
motion  amplitudes  were  less  than  0.5  feet  and  speed  deviation  less 
than  eight  percent  about  the  mean  in  the  larger  seastates. 

The  principal  cause  of  the  speed  decrease  is  presum- 
abJy  the  large  rudder  deflections  required  to  keep  the 
model  on  course  in  the  following  and  quartering  seas.  A plot 
of  the  RMS  rudder  deflections  versus  significant  wave  height, 
speed,  and  headings  showed  a large  amount  of  scatter.  This 
points  to  the  need  for  some  further  experiments  expressly  to 
evaluate  the  speed  decrement  due  to  rudder  motions. 

Motor  RPM  was  not  regularly  recorded  due  to  channel  limitations. 
It  was  checked,  however,  and  found  to  he  sufficiently  constant. 

The  heave  and  sway  accelerometers  ( + 2g  ) were  adjusted 
to  have  ranges  of  approximately  two-thirds  of  a g.  Subsequent 
analysis  suggt  ts  that  there  may  have  been  somewhat  more 
structural  damping  involved  in  the  heave  acceleration 
measurements  at  high  frequencies  than  anticipated.  The  surge 
accelerometer  ( + lg)  had  a higher  gain  factor  and  is  believed 
to  have  worked  well. 

Pitch  and  roll  angle  measurements  appear  to  be  accurate. 

The  full  scale  yaw  range  needs  to  be  large,  and  it  should  be 
remembered  that  very  small  yaw  motions  are  only  a small  per- 
centage of  the  recorded  range. 


In  order  to  compare  the  experimental  results  with  cal- 
culations using  an  inviscid  linear  theory,  the  actual- roll 
damping  of  the  model  is  required.  The  damping  coefficient 
was  computed  from  two  sets  of  experiments,  one  at  zero  speed 
and  one  at  forward  speed. The  analytical  procedure  and  results 
are  reported  in  Appendix  B. 
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Figure  III-9  Model  Speed  vs  Significant  Wave  Height 


C.  Observations 

A pronounced  difference  was  observed  between  the  two  models 
at  the  relative  frequency  of  occurrence  of  low  cycle  resonance 
at  the  first  Mathieu  frequency.  This  first  unstable  frequency 
(of  encounter)  it  will  be  recalled  is  equal  to  twice  the  natural 
frequency  of  roll.  During  the  experiments * it  was  observed 
that  the  Challenger  model  experienced  low  cycle  resonance  at 
this  frequency  relatively  frequently*  in  perhaps  one-third  to 
one  half  of  all  capsizes.  The  SL-7  model  on  the  other  hand 
experienced  virtually  no  capsizes  at  the  first  Mathieu  unstable 
frequency. 

The  explanation  is  obtained  by  reference  to  Figures  III-10 
(a)  and  (b)  where  righting  arm  curves  in  still  water  and  waves 
of  various  heights  are  plotted  for  the  two  models.  A pronounced 
difference  is  seen  in  the  nature  of  the  curves  for  the  two 
different  hull  forms.  In  the  case  of  the  Challenger , the 
trough  amidship  and  crest  amidship  curves  are  relatively 
symmetrically  disposed  about  the  still  water  curve.  That  is* 
the  increase  in  righting  arm  for  the  trough  amidships  case  is 
about  equal  to  the  decrease  in  righting  arm  for  the  crest 
amidships  position  of  the  wave  of  the  same  height.  In  contrast* 
for  the  SL-7  the  entire  family  of  crest  amidships  curves  nearly 
coincide  with  the  still  water  curve  up  to  an  angle  of  heel  of 
about  18  degrees*  while  all  of  the  trough  amidships  curves  are 
much  steeper  and  higher  over  the  entire  range  . It  appears 
from  the  nature  of  these  curves  that  the  SL-7  operating  in 
regular  following  seas  would  experience  a mean  righting  arm 
curve  lying  midway  between  the  trough  amidships  and  crest  amid- 
ships curves,  which  would  have  higher  values  than  the  still  water 
curve.  In  other  words,  for  the  SL-7  model,  the  apparent  mean 
ordinates  of  the  righting  arm  curve  in  following  seas  are  greater 
than  the  calm  water  values. 

This,  in  turn*  influences  the  frequency  of  wave  encounter 
at  which  Mathieu  type  resonance  is  experienced.  While  the  Challenger 
displayed  such  resonant  behavior  in  the  roll  motion  at  an  encounter 
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frequency  of  almost  exactly  twice  the  natural  roll  frequency  in 
caln  water.  The  SL-7  model,  on  the  other  hand,  experienced 
almost  no  cases  of  roll  motion  at  half  the  encounter  frequency, 
but  did  show  large  motion  and  frequent  capsizes  at  a roll 
frequency  equal  to  the  encounter  frequency.  The  explanation 
is  given  by  the  difference  between  the  two  vessels  righting 
arm  curves  in  following  seas.  The  Challenger's  natural  frequency 
in  waves  is  nearly  equal  to  that  in  calm  water  because  of  the 
symmetry  in  the  righting  arm  curves  for  the  crest  and  trough 
amidships  positions.  Mathieu  resonance  at  the  first  unstable 
frequency  was  observed  to  occur  quite  frequently  during  the 
experiments  because  the  combination  of  model  speed  and  average 
wave  lengths  resulted  in  the  occurrence  of  encounter  frequencies 
in  this  range.  For  the  SL-7  on  the  other  hand,  the  mean 
effective  righting  arms  in  waves  was  generally  substantially 
greater  than  in  still  water.  This  resulted  in  a relatively  high 
"apparent"  natural  frequency  of  roll  in  waves.  The  combinations 
of  model  speeds  and  wave  lengths  in  which  tests  were  conducted 
resulted  in  frequencies  of  encounter  which  were  more  nearly 
equal  to,  rather  than  twice , the  natural  frequency  of  roll. 

This  observation  was  borne  out  by  the  results  obtained 
by  numerically  integrating  the  approximate  equation  of  motion 
(11-14)  which  was  described  in  an  earlier  section.  Mathieu 
instability  in  the  results  of  the  numerical  solution  was  ob- 
served to  occur  at  a frequency  of  excitation  equal  to  twice  the 
(small  angle)  natural  frequency  which  would  be  computed  by 
using  the  mean  initial  slope  of  the  trough  amidships  and  crest 
amidships  righting  arm  curves  rather  than  the  still  water  GM. 

All  capsizes  observed  in  the  course  of  our  experiments 
occurred  in  following  or  quartering  seas,  and  these,  as  was 
pointed  out  in  an  earlier  section,  are  the  headings  of  the  ship 
relative  to  waves  in  which  the  stability  is  most  strongly 
affected  by  the  waves.  Trom  observing  the  capsizes  and  the 
motion  picture  records  of  them,  it  was  clear  that  the  attenuation 
of  stability  by  the  v/aves  played  a very  important  role  in  nearly 
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all  capsizes.  Further,  it  was  possible  to  distinguish  three 
distinct  modes  of  capsizing  which  may  be  described  as  follows. 

Mode  1:  Low  Cycle  Resonance.  This  refers  to  an  oscillatory 

rolling  motion  which  builds  up  rapidly,  i.e.,  in  two  to  five 
cycles,  to  a very  large  amplitude,  culminating  in  a capsize. 

The  phenomenon  appears  to  occur  in  approximately  the 
following  sequence.  The  model,  while  operating  in  following 
or  quartering  seas,  encounters  a group  of  especially  steep  and 
regular  waves.  When  the  crest  of  a wave  is  about  amidships, 
the  stability  of  the  model  is  greatly  reduced  and  it  takes 
a large  roll.  This  wave  moves  on  past  the  model  and  a 
trough  comes  into  the  amidships  position  while  the  model  is 
heeled  over,  resulting  in  sharply  increased  stability.  This 
causes  the  model  to  'snap"  back  upright,  acquiring  a high  roll 
angular  velocity  bv  the  time  it  reaches  the  upright  position. 
Another  wave  crest,  meanwhile,  is  moving  into  the  amidships 
position,  resulting  in  diminished  stability  once  again  as  the 
ship  starts  rolling  past  upright  and  to  the  other  side.  The 
ship  then  rolls  far  over  to  the  side  against  a diminished 
restoring  moment.  If  now  another  trough  moves  into  the  amid- 
ships position  with  the  correct  timing,  the  roll  will  be 
stopped  and  the  model  snaps  upright  again.  This  process  con- 
tinues until  either  the  model  capsizes  or  it  mc  js  out  of  the 
wave  group  and  the  motion  dies  down.  This  mwae  of  capsize  is 
seen  to  be  related  to  the  Mcthieu  motion  instability.  It 
results  directly  from  the  periodic  stability  variations  ex- 
perienced by  the  ship  moving  through  waves  and  in  its  most 
pronounced  form  takes  place  at  one-half  the  encounter  frequency, 
thus  at  the  first  Mathieu  unstable  frequency. 

Mode  2:  Pure  Loss  of  Stability. 

This  usually  occurs  in  a following  sea  at  high  speed.  The 
model  is  observed  to  encounter  one  or  more  very  steep  and  high 
waves  and,  with  little  or  no  preliminary  rolling  motion,  simply 
loses  all  stability  when  a crest  moves  into  the  amidships 
position  and  "flops"  over.  The  essential  prerequisite  for  this 
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to  occur  is  a model  speed  nearly  equal  to  the  wave  phase  velocity 
so  that  the  model  remains  almost  stationary  relative  to  the  crest 
for  a sufficient  length  of  time  to  capsize.  The  necessary  wave 
would  be  of  about  the  same  length  as  the  model  and  the  height 
would  be  sufficient  to  immerse  the  deck  in  the  crest  with  the 
model  upright.  This,  of  course,  implies  a high  model  speed  since 
a Froude  number  of  0.4  is  required  for  the  model  speed  to  be 
exactly  equal  to  wave  speed  in  waves  of  length  equal  to  model 
length.  From  motion  picture  records  of  several  capsizes  of  this 
nature,  it  appeared  that  a model  speed  lying  between  the  group 
velocity  (one-half  the  phase  velocity)  and  the  phase  velocity 
could  result  in  this  mode  of  capsize. 

Mode  3:  Broaching.  This  is  the  most  dynamic  mode  in 

appearance,  and  has  received  the  most  attention  in  the  previous 
literature.  In  this  mode  of  capsize,  the  model  is  struck  from 
astern  by  three  or  four  steep  breaking  seas  in  succession.  As 
each  wave  strikes  it,  the  model  is  forced  to  yaw  off  course  to 
such  an  extent  that  the  steering  system  is  unable  to  correct  the 
heading  in  the  time  interval  between  waves.  The  breaking  seas 
striking  the  model,  combined  with  the  dynamic  heeling  moment 
resulting  from  the  turn,  combine-  to  cause  capsizing,  again  with 
the  crest  of  a wave  amidships.  The  essential  features  of 
broaching  are  the  breaking  waves  striking  the  model  in  series, 
and  the  large  heading  deviation  and  associated  angular  velo- 
city. 

On  several  occasions,  broaching  was  observed  to  occur  with- 
out capsizing  but  with  such  total  loss  of  directional  control  that 
the  model  swung  through  ninety  degrees  from  a following  sea  course 
to  beam  seas.  This  was  observed  to  occur  most  frequently  in  the 
light  displacement  condition  where  the  rudder  was  less  deeply 
immersed  and  therefore  less  effective.  On  several  occasions,  at 
low  speed,  the  model  was  totally  unable  to  regain  it  original  course, 
but  remained  in  the  beam  sea  condition  even  with  the  rudder  hard 
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IV.  LINEAR  ANALYSIS  AND  EXPERIMENTAL  COMPARISONS 

This  chapter  is  an  attempt  to  investigate  cl  e applicability 
of  linear  ship  motions  theory  to  the  case  of  extreme  following 
seas.  There  are  a number  of  obvious  reasons  for  pursuing  a 
linear  analysis.  No  other  three-dimensional  theory  is  so  well 
developed.  Strip  theory  has  yielded  surprisingly  good  comparisons 
with  measurements  despite  a flagrant  disregard  of  the  restrictions 
on  the  input  wave  frequency,  wavelength,  and  amplitude.  Capsizing 
ia  a rare  event  and  should  be  viewed  in  a probablistic  sense.  A 
spectral  description  of  the  motions  is  far  more  amenable  to  long 
term  extreme  motion  prediction  them  a time  series  derived  from 
a nonlinear  model. 

On  the  other  hand,  there  are  a number  of  reasons  why  one 
should  not  expect  particularly  good  agreement.  Strip  theory 
requires  the  input  wavelength  to  be  of  the  order  of  the  beam. 

For  a ship  operating  under  conditions  varying  from  beam  to  follow- 
ing seas,  the  significant  responses  stem  largely  from  the  excit- 
ation of  waves  ranging  in  wavelength  from  the  order  of  the  beam 
to  the  length  of  the  vessel.  The  solution  is  singular  for  the 
unrestored  motions  in  stern  seas  when  ui e=0 . Significant  vari- 
ations of  the  roll  hydrostatic  restoring  coefficient  are  not 
included  in  the  computations  since  linearized  hydrostatics  are 
determined  from  the  calm  water  properties  of  the  waterplane 
section.  A rudder-autopilot  system,  having  important  nonlinear 
properties,  may  be  of  some  consequence  in  large  following  seas. 
Obtaining  a good  linearized  estimate  of  the  roll  damping 
(including  the  viscous  contribution)  and  a measure  of  the  waves 
actually  encountered  also  introduce  complications. 

These  are  all  potential  problems.  The  significance  of  each 
one  must  be  evaluated  under  actual  test  conditions.  The  aim 
here  is  to  see  if  there  is  sufficient  agreement  between  theory 
and  experiment  to  warrant  further  refinements  in  the  linear 
prediction  procedure  and  to  determine  if  an  estimate  may  be  made 
as  to  the  bounds  to  which  this  theory  may  be  applied. 
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Linear  motion  predictions  were  made  using  the  Salvescn,  Tuck, 
Faltinsen  theory  (1970) . The  six  degree  of  freedom  ship  Motions  pro- 
gram, SEAWAY , developed  at  U.  C.  Berkeley , uses  a modified  ver  si  on  of 
the  Frank  close-fit  method  for  calculating  the  two-dimensional 
hydrodynamic  etJe'f ficients . first  and  second  "irregular 

frequencies"  have  been  removed.  Coefficients  at  the  desired 
encounter  frequencies  are  obtained  by  interpolation  from  a 
computed  tab^e.  The  program  will  generate  a Pierson-Moskowitz 
sea  specfrum  or  accept  either  a point  or  a directional  empirical  spectra 
as  input.  The  point  spectrum  can  be  spread  using  a cosine  to 
the  power  r distribution  if  desired.  An  additional  feature  is 
tne  Froude-Krylov  computation  of  surge  motion.  Surge  is  assumed 
to  be  a second  order  quantity  in  most  linear,  head  sea  ship 
motion  predictions.  This  is  not  the  case,  however,  for  long 
waves  in  following  and  quartering  seas  where  the  surge  exciting 
force  can  be  significant. 

Before  proceeding  directly  into  the  evaluation  of  the  ex- 
perimental results  and  theoretical  predictions,  it  is  of  value  to 
comment  upon  a few  of  the  previously  mentioned  potential  theoretical 
problems  and  to  discuss  examples  of  the  calculated  linear  transfer 
functions  with  regard  to  sensitivity  and  possible  singularities. 

Experimental  research  and  analysis  carried  on  within  the 
last  20  years  has  clearly  shown  that  linear  strip  theory  may  be 
used  to  predict  pitch  and  heave  motions  in  head  and  bow  seas  with 
reasonable  accuracy  even  at  high  Froude  numbers  and  large  sea 
states.  The  success  of  the  theory  is  believed  to  stem  largely 
from  the  predominance  of  the  reasonable  linear  pitch  and  heave 
hydrostatic  restoring  coefficients  in  the  solution  of  these 
coupled  equations  of  motion.  Variations  in  the  geometry  of  the 
waterplane  due  to  waves  or  a heel  angle  induce  dramatic  changes 
in  the  roll  hydrostatic  restoring  moment  so  that  as  indicated  in 
Figure  A- 3,  the  righting  arm  curve  is  not  linearly  related  to 
the  roll  angle.  In  a following  seaway  the  restoring  moment  will 
actually  vary  sinusoidally  between  the  two  extremes.  It  is  necessary 
then  to  choose  a suitable  value  for  the  GM  so  that  GMb  will  give 
a reasonable  approximation  of  the  righting  arm  up  to  15  degrees 
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or  some  appropriate  figure  of  the  RMS  roll  amplitude. 

Viscous  effects  become  apparent  in  roll  in  the  neighborhood 
of  resonance.  To  account  for  this,  a constant  is  added  to  the 
computed  wave  damping  coefficient  which  is  equal  to  the  dirference 
between  the  total  measured  damping  and  the  computed  damping  at 
resonance.  Experiments  conducted  elsewhere  have  demonstrated 
that  roll  damping  is  increased  at  forward  speed.  To  verify  this 
change  and  to  determine  the  incremental  variation  with  speed, 
values  for  the  total  damping  were  determined  from  roll  decay 
experiments  conducted  with  the  SL-7  model  in  calm  water  at  zero 
speed  and  at  the  three  speed  settings.  The  linear  damping  term 
at  forward  speed  was  approximately  three  times  greater  than  the 
zero  speed  ertimate.  Large  irregular  variations  in  the  total 
damping  occurred  between  the  three  r’roude  numbers.  Values  for 
the  increased  damping  calculated  from  a relationship  which  was 
linear  in  forward  speed  approximated  the  general  increase  which 
was  measured.  As  can  be  expected  changes  due  to  GM  were 
significant  due  to  the  predominance  of  the  roll  period  in  the 
calculations.  A more  complete  discussion  of  tne  roll  decay 
experiments  and  of  the  calculations  of  the  total  roll  damping  term 
is  presented  in  Appendix  B. 

During  experimental  runs  in  San  Francisco  Bay,  the  SL-7 
model  was  observed  to  experience  large  heel  angles  due  to  the 
effect  of  the  rudder,  particular ly  at  high  speed.  If  the  model 
broached,  as  it  was  prone  to  do  in  a large  following  or  quartering 
seaway,  and  if  its  speed  was  not  greatly  reduced,  then  the  model 
would  heel  significantly  as  the  rudder-autopilot  system  tried  to 
correct  cne  course.  This  induced  roll  moment  coupled  with  a 
succeeding  large  wave  was  often  sufficient  to  capsize  the  model. 
Clearly,  the  effect  of  a rudder-autopilot  system  is  important  ii. 
ship  motion  calculations  in  large  seas,  at  least  for  a fine 
domain  solution.  But  how  would  such  a system  affect  linear  ship 
motion  predictions?  Intuitively,  it  was  felt  that  the  magnitude 
of  the  related  motions  would  be  reduced.  To  verify  this  assump- 
tion, a linear  rudder-autopilot  system  was  added  to  SEAWAY . 

Normally,  any  autopilot  system  is  designed  to  operate  in  a 
linear  manner.  The  rudder  force  or  lift  is  proportional  to  the 
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rudder  angle,  6,  which  in  turn  is  a function  of  the  yaw  gain  and 
the  yaw  rate  gain. 

6 = k ip  + k \p 
1 2 

The  nonlinearities  arise  due  to  any  deadband  in  the  system, 
due  to  time  lags  inherent  in  the  mechanical  system  and  due  to  a 
maximum  limitation  of  the  rudder  angle. 

It  is  possible,  though,  to  use  a least  squares  fit  method 

to  determine  values  for  k and  k which  give  a better  approximation 

12 

of  the  actual  operation  of  the  rudder-autopilot  system.  See 
Appendix  C.  With  acceptable  values  for  k ^ and  k the  rudder 
coefficient  of  lift,  and  the  center  of  pressure  of  the  rudder, 
the  appropriate  coefficients  may  be  determined  and  added  to  the 
coupled  equations  of  motion  in  roll,  yaw,  and  sway. 

The  response  amplitude  operators  (RAO)  are  the  motion  components 
computed  for  unit  wave  amplitudes.  These  were  calculated  for  one 
condition  of  the  SL-7  model  with  and  without  an  equivalent  linear 
rudder-autopilot  system.  Results  are  presented  in  Figures  IV-1.1 
through  IV-2.3.  As  anticipated,  the  yaw  and  sway  motions  were  decreased. 
Surprisingly,  the  roll  RAO's  showed  a small  but  definite  increase. 

Note  that  the  results  are  presented  for  wave  frequencies  and  not  for 
encounter  frequencies.  The  latter  would  require  a folding  of  the 
spectra  or  a concentration  of  the  energy  within  a narrower  band. 

Calculation  of  the  RAO's  using  SEAWAY  was  limited  by  a 
maximum  of  24  wave  frequencies.  In  order  to  avoid  missing 
resonant  peaks  in  the  spectra  and  to  assure  that  linear  inter- 
polation between  spectral  ordinates  would  be  reasonable,  a maximum 
wave  frequency  of  8.0  radians  per  second  was  chosen.  Selection 
of  this  upper  limit  was  due  to  the  fact  that  most  of  the  measured 
wave  point  spectra  contained  very  little  energy  beyond  this 
frequency.  Even  so,  some  of  the  RAO's  do  show  a somewhat 
squarish"  peak  which  would  indicate  that  the  peak  was  rather 
narrow  and  that  the  resonant  frequency  occurred  midway  within  the 
interval. 
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The  heave  and  pitch  RAO's  shown  in  Figures  IV-3.1  and  IV-3.2 
display  a gradual  change  in  shape  and  magnitude  between  different 
headings  with  respect  to  the  wave  direction.  The  largest  vari- 
ations occur  near  beam  seas.  For  x = 90°  the  response  is  very 
broad  and  therefore  susceptible  to  excitation  from  higher 
frequencies,  even  after  folding,  since  for  this  condition 

* w.  This  point  is  significant  for  theoretical  predictions 
for  a quartering  seas  course  when  the  directional  srectrum  is 
reasonably  broad.  Resulting  motion  spectra  may  be  heavily 
influenced  by  high  frequency  components  from  the  beam  direction. 

In  order  to  check  SEAWAY'S  output,  comparisons  were  made 
with  published  experimental  data.  Good  agreement  was  obtained 
for  pitch  and  heave  RAO's  for  all  directions  and  a range  of 
speeds. 

Unfortunately,  there  is  a lack  of  experimental  results  for 
sway,  yaw  and  roll  RAO's  which  does  not  permit  a reasonable 
comparison  to  be  made  with  SEAWAY's  output.  Sway  and  yaw  RAOs 
(Figures  IV-1.1  and  IV-1.3)  are  generally  characterized  by  sharp 
peaks  at  the  resonant  frequency  for  each  heading  or  direction. 

The  roll  RAOs  are  also  closely  related  to  the  natural  frequency. 
However,  the  spectra  shown  in  Figure  IV-1.2  exhibit  much  broader 
qualities  and  often  a second  peak  due  to  the  coupling  effects 
of  yaw  and  sway. 

One  aspect  of  linear  strip  theory  predictions  which  has  not 
been  touched  upon  is  the  seaway  or  directional  wave  spectra. 

This  is,  of  course,  along  with  speed  and  direction,  one  of  the 
major  variables  in  linear  theory  predictions.  The  directional 
spectra  used  as  input  for  SEAWAY  were  computed  by  Capon's  maximum 
likelihood  filter  from  the  point  spectra  measured  with  the  four 
probe  wave  buoy  between  various  experimental  runs.  As  discussed 
in  Chapter  III  the  measurement  and  computation  of  the  directional 
spectra  has  been  very  successful.  Possible  inaccuracies  may  arise 
due  to  the  step-measurement  of  the  wave  height  or  due  to  the 
spreading  or  smoothing  which  is  done  in  calculating  a 360  degree 
spectra  from  a number  of  point  spectra.  However,  these  points 
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die  of  relatively  minor  importance  compared  with  two  other 
considerations.  First,  point  spectra  were  measured  at  one 
location  while  the  model  during  a very  long  run  may  have  covered 
a distance  of  2/3  of  a mile  or  more.  Second,  the  seaway  was  not 
measured  simultaneously  with  each  run  but  rather  at  infrequent 
intervals  between  runs.  All  model  runs  were  bracketed  by  wave 
spectral  measurements.  Interpolation  between  wave  spectra  for 
the  tine  of  the  model  run  should  be  possible  for  many  of  these 
experiments.  However,  little  is  known  about  the  changes  in 
spectral  shape  with  time.  Comparisons  of  linear  theory 
predictions  with  experimental  results  are  made  using  the  direc- 
tional spectra  before  and  after.  Good  agreement  would  require 
that  the  computed  model  motions  bracket  the  experimental  motions 
in  a manner  consistent  with  the  time  inter  -Is  between  measured 
events. 

The  magnitude  of  the  variation  in  the  seaway  during  a long 
run  is  somewhat  difficult  to  assess.  Courses  were  chosen  so 
that  the  model  would  pass  near  the  wave  buoy  approximately  half 
way  through  a 900  second  long  run.  Disregarding  the  necessary 
interpolation  in  time,  the  wave  spectra  may  be  assumed  to  represent 
the  mean  wave  conditions  of  a long  run.  Judging  from  the  time 
hi  stories  of  the  measured  motions  and  from  observations  of  the 
waves  during  experiments  this  assumption  is  not  altogether  unrea- 
son ab  ie  . 

The  significance  of  using  the  seaway  measured  before  and 
after  a run  instead  of  an  interpolated  spectrum  between  the  two 
will  be  discussed  later  on  in  light  of  the  actual  comparisons 
between  theory  and  experiment. 

One  of  the  aims  of  studying  and  comparing  linear  theory 
predictions  with  measured  experimental  motions  was  to  attempt  to 
determine  an  upper  bound  or  limit  to  which  the  theory  may  be 
applied.  Certainly  the  comparisons  will  not  represent  a complete 
range  of  the  input  variables,  particularly  the  sea  state.  Vir- 
tually all  of  the  experiments  were  conducted  in  rather  large 
following  and  quartering  seas.  Comparisons  with  linear 
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computations  in  many  cases  are  not  particularly  fair  since  large 
waves  and  motions  present  a totally  new  geometry  and  substan- 
tially altered  stability  characteristics.  The  model's  speed  and 
heading  were  subject  to  variation  such  as  a real  ship  might 
experience  in  heavy  seas.  Courses  and  speeds  used  for  linear 
calculations  represent  the  best  possible  choice  for  the  mean 
conditions  of  a run.  In  spite  of  the  seeming  inability  to  deter- 
mine the  input  parameters  exactly  due  to  their  variation  and  of 
the  general  disregard  for  some  of  the  basic  restrictions  of  linear 
theory,  it  was  hoped  that  the  predictions  would  still  give 
reasonable  estimates  of  the  measured  motions  or  values  which 
showed  a definite  comparable  trend  that  might  possibly  be  ex- 
plained and  therefore  corrected.  Comparisons  of  the  computations 
for  all  six  degrees  of  freedom  with  measured  values  have  shown 
a large  amount  of  scatter.  In  fact,  the  agreement  in  general 
has  been  poor.  In  order  to  indicate  the  degree  of  agreement,  to 
point  out  some  possible  sources  of  error  and  to  assess  qualita- 
tively the  value  or  usefulness  of  the  theory  for  predictions  of 
ship  motions  in  extreme  seas,  a few  experimental  runs  will  be 
discussed. 

Rather  than  pick  a number  of  experimental  runs  at  random  for 
comparison,  a matrix  of  the  runs  was  formulated, based  upon  a mini- 
mum length  of  500  seconds  and  the  "reasonableness"  of  the  data 
recorded,  to  facilitate  an  analysis  over  a range  of  soastates. 

The  seastate  used  to  characterize  a run  and  determine  its  place 
in  the  matrix  was  found  by  a linear  interpolation  between  the 
seastate  measured  before  and  after  the  run.  Most  of  the  long 
runs  with  the  OL-7  model  were  made  for  speeds  2 and  3 at  uV 
conditions  4 and  6.  At  lower  speeds  the  model  was  less  likely 
to  capsize  and  therefore  many  of  the  long  runs  for  large  seaways 
occur  at  low  speeds.  However,  for  such  conditions  it  was  observed 
that  the  model  often  experienced  difficulty  in  maintaining  its 
prescribed  course  and  once  knocked  about  to  beam  seas  could  not 
return  to  a quartering  or  following  seas  heading.  Linear  motions 
were  calculated  for  speed  2 (7AM  and  for  speed  3 (!'■'>'  of  the 
matrix.  Generally  the  results  from  the  two  conditions  w^re 
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were  similar.  For  the  sake  of  brevity  only  the  surge,  pitch, 
heave,  roll  and  yaw  motions  for  a few  runs  for  the  former 
condition  will  be  discussed.  The  results  are  presented  in 
Figures  IV-4.1  through  IV-4.22. 

Translational  surge  and  heave  motions  were  calculated  as 
displacements.  In  order  to  compare  them  with  the  measured 
accelerations  their  spectra  were  multiplied  by  u4,  thus  emphasizing 
any  motions  at  higher  frequencies. 

Surge  accelerations  were  measured  for  only  three  of  the  six 
runs  and  are  presented  in  Figures  IV-4.5,  4.13,  and  4.18.  The 
computed  surge  motions,  which  were  determined  solely  from  the 
Froude-Krylov  forces,  display  two  very  distinct  peaks.  The  first 
peak  occurs  at  fairly  low  encounter  frequencies  and  is  significantly 
larger  than  the  second.  Note  that  no  similar  experimental  motions 
were  measured;  thus  pointing  out  that  neglect  of  the  damping  due 
to  the  forced  motion  of  the  model  is  a serious  omission.  The 
second  peak  shows  a slight  agreement  both  in  magnitude  and  in 
frequency  range  with  the  measured  peak. 

Both  measured  and  calculated  heave  accelerations  showed  an 
upward  trend  toward  the  higher  frequencies.  See  Figures  IV-4.1,  4.9, 
4.14  and  4.19.  This  characteristic  is  readily  attributable  to  beam 
sea  components  cf  the  directional  spectra.  Elimination  of  the 
ramps  from  computed  results  could  be  effected  by  truncating  the 
directional  spectra  to  remove  all  input  from  these  directions. 

In  general  comparison  between  measured  and  calculated  heave 
results  were  disappointingly  poor.  Unfortunately,  no  reasonable 
explanation  may  be  offered  for  these  anomalies,  particularly  in 
lieu  of  the  fact  that  good  agreement  between  measured  and  cal- 
culated heave  motions  was  obtained  for  head  sea  tests.  As  mentioned 
earlier,  comparison  with  other  published  heave  measurements  has 
also  been  good  even  for  a wide  range  of  speeds  and  headings. 

Agreement  between  measured  and  calculated  pitch  motions  was 
reasonable  for  most  runs.  The  computed  spectra  usually  exhibited 
a much  narrower  spectral  peak. 

The  encounter  frequency  at  which  the  computed  pitch  spectrum 
peaks  is  a useful  measure  of  the  speed  of  the  model.  This  value 
together  with  the  results  obtained  from  the  relative  velocity 
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transducer  permitted  an  estimation  to  be  made  of  the  decrease 
in  speed  due  the  waves.  In  essence  the  speed  input  to  SEAWAY 
was  varied  for  a number  of  runs  until  the  peaks  of  the  measured 
and  calculated  pitch  spectra  coincided.  The  influence  of  the 
magnitude  of  the  seaway  upon  the  speed  of  the  model  is  presented 
as  a graph  of  relative  water  speed  measured  versus  significant 
wave  height  in  Figure  III-9.  Due  to  the  scatter  obtained,  a 
constant  speed  loss  was  assumed  for  each  speed  setting.  For  the 
SL-7  model  in  the  light  condition,  the  speeds  assumed  for  linear 
computations  were  4.40,  5.75,  and  7.10  fps.,  corresponding  to 
reductions  of  0.25,  0.25  and  0.35  ips.  for  speed  settings  1,  2, 
and  3 respectively. 

The  calculated  pitch  spectra  for  run  11.05  (Figure  IV-4.4) 
differ  significantly  due  to  a very  large  increase  in  the  seaway 
between  successive  measurements.  The  development  of  the  seaway 
was  probably  not  gradual.  As  such,  interpolation  between  the 
directional  spectra  would  probably  not  have  been  valid.  Figure 
E-1.10  gives  a time  sequence  of  the  runs  and  wave  records  for  the 
day  and  some  insight  as  to  the  gio.rth  of  the  seaway. (See  Figure  E-2.11) 

Like  pitch,  the  measured  roll  spectra,  shown  in  Figures  IV-4.2, 
4.6,  4.10,  4.15  and  4.20  exhibited  broader  peaks  than  those 
predicted  by  linear  theory.  The  total  energy  or  area  under  the 
curve  was  often  approximately  double.  Peak  values  for  the 
examples  were  reasonable  though  consistently  low.  This  is  be- 
lieved attributable  to  the  use  of  a damping  coefficient  which  was 
calculated  using  the  natural  period  of  small  roll  rather  than  the 
period  of  roll  at  large  angles.  The  latter  is  approximately 
10  percent  smaller  and  would  have  produced  a similar  10  percent 
reduction  in  the  damping  coefficient  which  in  turn  would  have 
increased  the  calculated  roll  peaks  slightly.  See  Appendix  B 
for  a more  detailed  discussion  of  the  roll  damping. 

Both  measured  and  computed  yaw  spectra  displayed  peaks  at 
the  low  frequency  end  of  the  spectrum.  However,  the  difference 
in  magnitude  was  quite  large  so  that  any  further  comparisons  are 


50 


seemingly  impractical.  In  fact,  only  run  18.03  shown  in  Figure 
IV-4 . lb  displays  some  correlation  in  shape  between  the 
measured  and  computed  spectral  curves. 

Based  upon  the  foregoing  results,  linear  theory  does  not 
appear  to  be  a useful  tool  for  predicting  motions  for  all  six 
degrees  of  freedom  in  extreme  following  seas.  Unfortunately  the 
aim  of  assessing  the  bounds  or  limitations  of  the  theory  has 
not  been  realized  except  to  conclude  that  it  is  inadequate  for 
the  aforementioned  conditions.  To  accomplish  the  latter  task 
would  require  a controlled  series  of  experiments  in  which  the 
effect  of  each  variable  might  be  determined  separately. 

Although  it  may  not  be  feasible  to  use  a computed  spectrum 
for  a statistical  analysis  of  roll  in  extreme  following  seas 
due  to  the  discrepancy  in  spectral  area,  it  would  appear  that 
a useful  value  for  the  maximum  roll  angle  may  be  determined 
from  linear  theory  provided  that  a reasonable  "adjusted"  value 
for  the  W.  is  chosen  and  that  the  damping  coefficient  include 
the  effect  of  forward  speed  and  is  calculated  for  large  roll 
angles  as  proposed  in  Appendix  B. 
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Figure  IV--3.1 
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Figure  IV-4.9 
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V.  NUMERICAL  SIMULATION  OF  THE  MOTION  IN  ASTERN  SEAS 

A complete  analytical  solution  for  the  motion  of  a ship 
in  waves  requires  first  that  the  hydrodynamic  forces  acting 
on  the  ship  be  found.  Solutions  of  the  hydrodynamic  problems 
have  heretofore  been  obtained  only  under  the  assumption  of 
small  motion  amplitudes,  in  which  case  the  forces  acting  on  the 
ship  are  computed  as  though  the  instantaneous  position  of  the 
ship  differs  but  little  from  its  mean  position.  Such  an 
assumption  cannot  be  used  in  the  present  case  where  large 
deviations  in  position  from  the  mean  are  an  essential  feature 
of  chc  phenomenon.  Instead,  we  observe  that  at  high  speed  in 
following  and  quartering  seas  the  frequency  of  wave  encounter 
will  be  low  and  the  ship  motion  will  be  determined  largely  by 
the  hydrostatic  forces.  This  enables  us  to  retreat  from  the 
necessity  of  determining  the  hydrodynamic  forces  with  great 
accuracy  but  to  concentrate  instead  on  the  hydrostatic  forces 
which  nay  be  computed  for  the  exact  position  of  ship  and  waves. 
These  forces,  plus  additional  external  forces  representing, 
e.g.,  the  steering  and  controls,  plus  a simplified  approximation 
to  the  relatively  unimportant  hydrodynamic  terms  then  are  used 
as  tho  right  hand  side  of  the  rigid  body  equations  of  motion. 

A standard  numerical  procedure  is  employed  to  integrate  the 
equations  of  motion  leading  to  a step-by-step  approximation 
of  the  vessel's  motion. 

Formulation  of  the  Problem 

The  ship  is  assumed  to  behave  as  a rigid  body  having  six 
degrees  of  freedom.  Newton's  second  law  may  be  written  for 
the  body  in  the  form 

;vy  rn  V = f (V-l) 


ind 


7 j 


I u = g 


(V-2 ) 
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whore 

t = time, 

m = mass  of  the  body  (ship) , 
v = velocity  vector  of  the  mass  center, 
f = force  vector, 

I = inertia  matrix, 

w = angular  velocity  vector,  and 

g * moment  of  the  force  about  the  mass  center. 

The  force  and  moment  result  from  the  gravitational  force  acting 
at  the  mass  center  and  the  interaction  between  the  ship  and  the 
sea.  The  force  and  moment,  in  general,  depend  on  the  time  history 
of  the  position  of  the  ship  in  the  sea.  Under  appropriate  cir- 
cumstances, however,  this  history  may  be  characterized  by  the 
instantaneous  position,  velocity,  and  acceleration  of  the  ship. 
The  general  problem  is  nonlinear  in  the  motion  variables  in 
that  the  force  and  moment  are  nonlinear  functions  of  the  motion 
of  the  ship,  and  the  rate  of  change  of  angular  momentum  in 
equation  (V-2)  contains  nonlinear  terms.  As  noted  previously, 
we  shall  focus  on  an  exact  computation  only  of  the  hydrostatic 
part  of  the  force. 

B.  Coordinate  Systems 

Since  large  amplitude  motions  are  to  be  computed,  it  is 
necessary  to  clearly  define  the  relati^rships  between  several 
coordinate  systems  to  be  used  in  describing  the  ship  and  water 
motion.  The  coordinate  systems  described  below  are  right  hand 
Cartesian  systems. 

A Newtonian  reference  frame  is  formed  by  the  dxyz  system 
which  is  fixed  in  space  and  so  oriented  that  the  xz-plane  is 
the equilibrium  sea  surface,  and  the  y- axis  is  directed  upwards. 

A body  coordinate  system  Gxyz  is  fixed  in  the  ship  such 
that  the  origin,  G,  coincides  with  the  center  of  gravity  of 
the  ship.  In  a ship  of  usual  form,  the  x-axis  is  parallel  to 
the  baseline  and  directed  forward,  the  xy-plane  is  parallel  to 


r 
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b 


the  centerplane  of  the  ship,  the  y-axis  is  directed  upward  and 
the  s-axis  to  starboard. 


The  position  of  the  ship  mass  center,  G,  may  be 
specified  in  the  fixed  coordinate  system  by 

* " 
y * yc 


*G  * 


This  may  be  represented  by  the  vector 

. (•■ 

* - < y, 


(V-3) 


The  velocity  of  the  mass  center  is  represented  by  the  vector 


XG 

d , 

4 It  *G 

d , 

JFt  *G 


(V-4) 


Any  rotation  of  the  ship  coordinate  system  is  uniquely 
defined  by  the  modified  set  of  Eulerian  angles  described 
below.  These  angles  are  similar  to  the  ones  given  by 
Blagoveshchensky  (1962),  but  differ  from  the  ones  used  by 
Euler.  The  angles  are  defined  as  follows. 

Consider  the  ship  coordinates  in  a position  before  rotation 

with  the  x,  y , and  x-axes  parallel  to  the  fixed  x , yt  i-axes. 

This  is  the  orientation  of  Gx  u m in  Figure  V-l.  Rotate  the 

ii» 

triad  about  the  y^-axis  to  the  yaw  angle  $.  This  postions  the 
frame  as  Gx  y a ii  the  figure.  Next,  rotate  about  the  m -axis 

III  2 

to  the  pitch  angle  ip.  The  figure  shows  the  yawed  and  pitched 
orientation  as  Gxy  x . The  final  rotation  is  aJ>out  the  x>axis 

3 2 

to  the  roil  angle  6.  The  orientation  of  the  ship  coordinates, 
Gxyat  is  indicated  in  the  figure. 


1 
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The  Eulerian  angles  <|/  and  6 are  represented  by  the 
vector 


a 


(V-5) 


The  angular  volocities  about  the  ship  coordinate  axes  are  de- 
noted by  p,  qt  and  r corresponding  to  components  of  the 
angular  velocity  vector  along  the  x , y , and  s-axes.  These 
angular  velocities  may  be  expressed  in  terms  of  the  Eulerian 
angles  and  their  derivatives: 


r m ft  + H 5in  * 

q m cos  0 cos  ip  + sin  0 (V— 6) 

r * if  cos  0 • sin  ®c^8  ♦ 


The  notation  is  simplified  by  representing  the  angular  velocities 
by  the  vector 


u * 


If  we  define  the  matrix 


B 


1 sini|i 

0 cos  8 cos 

0 -sinOcost 


0 

sin6  , 
cos  6 


(V-7) 


(V-8) 
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and  note  that 


It 


da 

•*** 

It 


f 


(V-9) 


then  equations  (V-6)  are  represented  by 

da 

? 37  * s • 


(V-10) 


The  moments  and  products  of  inertia  in  the  angular 
momentum  equation  (V-2)  are  represented  by  the  matrix 


I 
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- 1 

XX 

xy 

xz 

-I 
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-I 

xy 

yy 

y* 

-i 

-i 

i 

xz 

yz 

zz 

(V-ll) 


The  moments  of  inertia  are  defined  as 


Ixx  = 1 (y2  + zZ)3 

I = ^ m'  (z  2 x2 ) j 

yy 


and 


1 zz  = 1 m'  (x  + y )t 


(V-12) 


where  the  summations  are  taken  over  all  particles  of  mass  m 
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comprising  the  ship.  The  products  of  inertia  are 


ZXy  * 1 «' 

I » Z mf  xz. 
xz  ' 

and 

I * l m'  uz. 
yz  * 


(V-13) 


These  moments  and  products  of  inertia  are  constants  in  the 
movina  ship  coordinate  system,  Gxyz. 

In  the  Gxyz  coordinate  system,  the  rate  of  change 
of  angular  momentum  is  given  by 


d 

It 


I u 


I ££  u)  + uj  x I to , 


(V-14) 


C.  Time  Domain  Integration 

The  equations  of  motion  are  solved  by  numerical  integration 
in  the  time  domain.  In  order  to  perform  the  integration  using 
standard  algorithms,  the  equations  of  motion  are  rewritten  as 
first  order  ordinary  differential  equations. 

The  position  of  the  ship's  center  of  gravity  is  determined 
by  the  linear  momentum  equation  (V-l) . This  equation  and 
equation  (V-4)  are  rewritten  as 

d - 1_ 

It  : * m t (V-15) 


and 


v 


(V-16) 
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In  equations  (V-15)  and  (V-16)  the  vectors  are  referred  to  the 
Newtonian  reference  frame,  Qxyz. 

The  rotations  of  the  ship  are  governed  by  the  angular 
momentum  equation  (V-2) . Combining  equations  (V-2)  and  (V-14) 
and  rewriting  equation  (V-10)  give 

u - I 1 [g  - u)  * Iw]  ( V— 17) 


and 

a = B_1U)  . (V-18) 

In  equation  (V-17)  the  moment  vector  is  referred  to  the  ship 
coordinate  system,  Gxyz. 

The  vector  equations  (V-15) , (V-16) , (V-17) , and  (V-18) 
form  a system  of  twelve  simultaneous  first  order  ordinary 
differential  equations  which  may  be  integrated  by  standard 
numerical  procedures.  The  current  version  of  the  program 
uses  a fifth-order  Adams  type  predictor-corrector  algorithm 
developed  by  Glauz  (1960) . This  features  a variable  time 
step  to  control  integration  error  and  an  interpolation  pro- 
cedure to  avoid  the  calculation  of  derivatives  at  the  times 
chosen  for  the  output  of  results.  At  each  time  step  in  the 
integration,  the  algorithm  predicts  values  for  the  instan- 
taneous position  of  the  ship  and  its  velocity.  The  remainder 
of  the  program  is  devoted  to  the  computation  of  the  force  and 
moment  used  to  evaluate  (V-15)  and  (V-17) . The  integration 
routine  corrects  the  position  and  velocity  based  on  the  values 
of  the  derivatives  (V-15)  through  (V-18) , but  if  the  error  is 
small  these  derivatives  are  not  recomputed. 

For  computation  purposes  equations  (V-15)  and  (V-17)  are 
combined.  The  generalized  force  vector  includes  both  force 


73 


and  moment: 


r > 
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f 2 

f — 

ft  i 
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0* 

(V-19) 


The  changing  angular  momentum  in  the  rotating  ship  coordinates  is 
included  in 


H 


where  the  moment  components  are 


(V-20) 


g 


< 


g - U x X u> 

% % v 


In  (V-19)  and  (V-20)  the  first  three  elements  are  components 
of  force  in  ship  coordinates,  and  the  last  three  are  com- 
ponents of  the  moment.  The  generalized  acceleration  vector  is 


fu  1 


a 


(V-21) 
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The  inertia  of  the  ship  is  represented  by  the  matrix 
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(V-22 ) 


Using  this  notation  the  momentum  equations  (V-15)  and  (V-17) 
are 


A*  a = f 


(V-2  3) 


After  determining  the  force  this  may  be  solved  for  the 
acceleration,  a . For  the  integration,  the  first  three 
components  of  a , the  linear  accelerations,  are  transformed 
into  the  fixed  coordinate  system  to  provide  velocity  de- 
rivatives (V-15)  , and  the  derivatives  of  the  angular 
velocity  (V-17)  are  the  last  three  elements  of  a . 

D.  Computation  of  Force  and  Moment 


The  present  version  of  the  computer  program  for  the 
time-domain  simulation  of  large  amplitude  ship  motions  assumes 
that  the  force  and  moment  acting  on  the  ship  may  be  modeled 
using  an  accurate  computation  of  the  hydrostatic  or  Froude- 
Krylov  forces  plus  approximations  to  the  hydrodynamic  forces. 

Since  large  amplitude  motions  and  finite  amplitude  waves 
are  assumed,  the  hydrostatic  restoring  and  coupling  coefficients 
computed  for  the  equilibrium  position  cannot  be  used.  It  has 
been  shown  by  Paulling  (1961)  and  others  that  there  can  be 
significant  variations  in  the  roll  restoring  moment  as  a wave 


I 
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progresses  ale  he  ship's  length  as  well  as  the  change  in  this 
moment  caused  large  amplitude  roll  angles.  The  Froude-Krylov 
force  that  is  computed  by  the  numerical  simulator  includes  both 
the  motion  exciting  force  and  the  restoring  force  and  moment 
that  result  from  the  situation  of  the  ship  in  the  system  of  waves 
at  any  time  step  during  the  simulation. 

The  sea  surface  elevation  is  given  by  the  sum  of  sinusoidal 
waves  in  the  fixed,  dxyz,  coordinate  system.  The  water  surface 
is  given  by 

N 

r\(ttx,z)  - l T\.(l,xtz) 

. , 1 (V-24) 

t = l 

where 

n = the  y coordinate  of  the  surface, 

and 

N = the  number  of  wave  components  (in  the  present  version 
of  the  program  0 < N <_  20)  . 

The  component  wave  amplitude  is: 

n.=  A.  cos  (xk . cos  6.  - zk . sin  6.  + <b.  - a.  t) 

it  i tv  i 

where 

= the  amplitude  of  the  i-th  wave, 

* the  circular  frequency, 

<J>.  = initial  phase  angle, 
k.  = c.2/g  - wave  number, 
g = the  gravitational  acceleration, 

and 

6^  = the  direction  of  the  wave  propagation. 
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The  wave  pressure  is 

p( 


N 

-p gy  + l Pi(ttx,yt  z) 

t=l 


k .y 
P ge  x* 


( V—  25) 


where 

P g = the  specific  weight  of  the  water. 

The  Froude-Kyrlov  force  and  moment  may  be  obtained  by 
integrating  the  pressure  over  the  entire  wetted  surface  of 
the  ship.  By  applying  Gauss'  Theorem  the  force  and  moment 
are  given  by  integrals  of  the  pressure  gradient  over  the 
submerged  volume  of  the  ship.  The  components  of  the  force 
and  moment  in  the  ship  coordinate  system,  Gxyz  , are 

- - !!S  H iv 

t,  - - Hr  dv 

t,  - - "t  If iv 

* nt  '*&  - 

f<  - 1"  - 

which  may  be  represented  by  the  generalized  force  vector,  f. 

In  equation  (V-26)  the  volume  element  is  dV;  f , f and  f 

12  1 

are  forces  in  the  x,  y,  and  3-directions;  and  f , f and 

4 5 

/ are  moments  about  these  x,  y3  and  z - axes.  The  integrals 

are  taken  over  all  volume  up  to  the  instantaneous  sea  surface 
within  the  envelope  of  the  ship. 


(V-26) 


V&)  dv 
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x^-)  dV 


77 


The  ship  hull  is  approximated  by  a number  of  polygons 
representing  the  stations  of  the  ship.  Each  polygon  is  in 
a plane  defined  by  a constant  value  of  x in  the  ship 
coordinate  system.  A maximum  of  24  line  segments  are  used 
for  each  closed  polygon  station  in  the  ship  and  a maximum 
of  25  stations  may  be  used.  The  stations  may  be  unsymmet- 
rical  and  unequal  station  spacings  are  permitted.  The  position 
of  the  center  of  gravity  may  be  in  any  fixed  position  relative 
to  the  ship. 


The  integrals  of  the  pressure  gradients,  velocities  and 
accelerations  over  each  station  made  up  with  straight  line 
segments  are  evaluated  exactly,  but  with  two  restrictions  on 
the  angle  of  pitch.  First,  the  pitch  angle  must  not  become  so 
large  as  to  cause  the  intersection  of  a station  plane  and  the 
instantaneous  sea  surface  to  define  multiple  legions  or  a 
closed  contour  in  the  station  plane.  Second,  the  magnitude 
of  a quantity  like  the  product  of  the  pitch  angle  and  the  slope 
of  all  component  waves  must  be  "small".  These  two  restrictions 
are  satisfied  for  vessels  of  usual  proportions  in  waves  with 
realistic  slopes. 


The  two-dimensional  forces  and  moment  at  each  station  are 
evaluated  as  functions  of  the  form 


fix) 


N 

V'  + l [C . cos  (k  . ' x ) 

*•  t 

i = l 


+ S . sin  (fc  . ’x)  ] 
i i 

(V-27) 


whe  re 

V ' is  obtained  from  the  static  (P gy)  part  of  the  pressure, 

C . and  S . result  from  the  sinusoidal  pressure  fluctuation 
for  the  i-th  wave  component,  and 

k . ' is  a projection  of  the  wave  number  onto  the  x-axis  of 
the  ship. 
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The  relative  magnitude  of  C.  and  S.  depends  on  the  phase 

of  the  waves  relative  to  the  center  of  ship  coordinates  at  each 
instant  of  time.  The  integrals  and  moments  of  the  functions 
like  ( V—  2 7)  along  the  length  of  the  ship  are  evaluated  with 
the  assumption  that  V' ,C^,  and  S ^ vary  linearly  in  x be- 
tween adjacent  stations  of  the  ship. 

E.  Wave  Diffraction,  Added  Mass  and  Damping 

Approximations  are  used  for  the  hydrodynamic  forces 
resulting  from  the  diffraction  of  the  incident  waves  and  motion 
of  the  ship.  These  forces  are  computed  using  constant  two- 
dimensional  added  mass  and  linear  damping  coefficients  for  each 
station  combined  with  averages  of  the  water  acceleration  and 
velocity  relative  to  the  stations.  The  hydrodynamic  approxi- 
mations are  not  expected  to  lead  to  serious  errors  if  the 
Froude-Krylov  force  is  dominant.  This,  as  noted,  is  expected 
to  be  the  case  in  the  most  severe  capsizing  situations  in 
following  or  quartering  seas. 

The  hydrodynamic  force  resulting  from  the  diffraction  of 
the  waves  is  approximated  in  the  following  manner.  Two  di- 
mensional added  mass  and  damping  coefficients  for  heave,  sway, 
roll  and  roll-sway  coupling  are  entered  into  the  program  as  con- 
stants for  each  station  of  the  ship.  Each  time  the  sectional 
Froude-Krylov  forces  are  computed  by  integrating  the  pressure 
gradient  over  a station  of  the  ship,  average  values  of  vertical, 
hortizontal  and  "roll"  water  velocities  and  accelerations  are 
also  evaluated  for  the  station.  The  "roll"  velocity  and  accel- 
eration components  are  the  first  and  second  time  derivatives  of 
the  slopes  of  constant  pressure  lines  in  the  plane  of  the  station. 
The  two-dimensional  coefficients  are  scaled  by  the  instantaneous 
submerged  area  of  the  station,  and  the  products  of  the  average 
water  velocities  and  accelerations  with  these  coefficients  yield 
two-dimensional  diffraction  forces  which  are  added  to  the  two- 
dimensional  Froude-Krylov  forces  before  the  longitudinal  integra- 
tions are  performed. 
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The  force  resulting  from  ship  velocities  in  heave,  sway, 
roll,  yaw  and  pitch  are  computed  in  the  above  manner  using  the 
scaled  damping  coefficients  and  the  components  of  ship  velocity 
at  each  station.  Since  the  linear  and  angular  accelerations  of 
the  ship  are  unknown  when  the  forces  are  being  computed,  the 
force  resulting  from  the  ship  acceleration  cannot  be  omputed 
with  the  same  procedure.  Instead,  a matrix  of  three-dimensional 
added  mass  coefficients  is  computed  using  the  scaled  two-  dimensional 
added  mass  coefficients.  The  product  of  this  added  mass  matrix 
and  the  vector  of  accelerations,  gives  the  required  force  vector. 

The  generalized  acceleration  vector,  a , is  defined  by  equation 
(V-21) . The  hydrodynamic  force  and  moment  resulting  from  this 
acceleration  is 


h = -A"  a ( V-2 8) 

where  A"  is  the  added  mass  matrix.  The  use  of  longitudinal 
moments  of  the  two-dimensional  added  mass  coefficients  leaves 
the  elements  ,A'.  which  are  related  to  surge  (i  or  j equal 
to  1)  undefined.  They  are  taken  as  zero. 

In  order  to  use  this  hydrodynamic  force,  the  momentum 
equation  (V-23)  is  rewritten  as 

A'  a = f*  + h (V-29) 

where  the  right  hand  side  is  the  sum  of  all  forces  on  the  ship. 
The  inertia  matrix  for  the  ship  is  A' . 

The  substitution  of  equation  (V-28)  into  equation  (V-29) 
gives  the  momentum  equation  including  added  mass: 

A a = f'  ( V- 30 ) 


wher'1 


A = A"  + A'  . 


Equation  (V-30)  is  solved  to  obtain  the  accelerations  required 
for  integration.  • 
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F.  Steering  System 

The  steering  system  of  a typical  ship  consists  of  three 
components — an  autopilot,  steering  machinery  and  the  rudder. 

The  autopilot  computes  a rudder  angle  which  should  correct  or 
prevent  errors  in  the  ships  heading.  The  machinery  attempts 
to  rotate  the  rudder  to  the  angle  specified  by  the  autopilot. 
This  rotation  is  mechanically  limited  to  some  maximum  rate  of 
rotation  and  some  maximum  rudder  deflection  angles.  The  rudder 
acts  as  a lifting  surface  in  the  water  which  generates  forces 
and  moments  on  the  ship  which  are  used  to  maintain  the  desired 
course.  The  steering  system  for  the  numerical  simulation  incor- 
porates these  three  components . 

The  autopilot  model  computes  a required  rudder  angle  which 
is  a linear  combination  of  the  yaw  rate,  the  yaw  angle  (heading 
error),  and  the  time  integral  of  the  yaw  angle.  The  propor- 
tionality factor  or  "gain  parameter"  for  any  of  these  heading 
functions  may  be  set  to  zero  resulting  in  a simpler  autopilot. 
For  example,  the  yaw  integral  gain  parameter  is  zeroed  to 
simulate  the  autopilots  used  for  the  Challenger  and  SL-7 
models  that  were  run  on  San  Francisco  Bay.  The  yaw  angle  used 
for  the  heading  error  and  yaw  integral  is  the  angle  <f>  which 
is  always  measured  about  a fixed  vertical  axis.  The  yaw  rate 
is  measured  in  the  ship  coordinates  rather  than  about  a fixed 
vertical  axis.  This  is  the  same  as  used  in  the  ship  model 
autopilots  but  it  may  differ  fxom  that  used  in  some  full  sized 
autopilots.  Two  "^ead  band"  parameters  are  provided  in  the 
simulated  autopilot.  If  the  magnitude  of  the  heading  error  is 
less  than  the  fits  dead  band  parameter,  the  autopilot  will 
require  a zero  rude  r angle  rather  than  the  value  computed  using 
the  gain  parameters.  This  type  of  dead  band  is  typical  of  the 
"weather"  adjustment  c ship  autopilots.  The  othar  dead  band 
available  in  the  simulate  as  incorporated  in  the  ship  model 
autopilots.  With  this  form  of  dead  band  the  magnitude  of  the 
required  rudder  angle  is  reduced  by  the  value  of  this  dead  band 
parameter. 
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The  steering  machinery  is  modeled  with  constant  values 
for  the  rate  of  rotation  of  the-  rudder  and  limits  to  the 
magnitude  of  rudder  deflection.  If  the  autopilot  requires 
a rudder  deflection  rate  that  is  less  than  the  machinery  rate 
and  an  angle  that  is  less  than  the  mechanical  stops,  the 
machinery  simulation  sets  the  rudder  to  the  autopilot  angle. 
Otherwise,  the  rudder  angle  lags  that  computed  by  the  autopilot. 

The  rudder  is  simulated  by  a vertical  line  (in  ship 
coordinates)  through  an  effective  center  of  the  rudder.  The 
average  of  the  water  velocity  relative  to  points  on  this  line  is 
computed  at  each  time  st.-p.  The  water  velocity  across  this 
line  is  the  superposition  of  the  motion  of  the  ship,  the 
motion  of  the  water  particles  in  the  waves,  and  a constant 
wake  velocity.  The  lift  and  drag  forces  on  the  rudder  are 
assumed  to  be  proportional  to  the  instantaneous  submerged 
rudder  area,  the  square  of  the  average  relative  water  velocity, 
and  to  the  angle  of  incidence  between  the  rudder  and  the 
average  water  velocity.  Rudder  lift  is  limited  by  a stall  angle. 
The  lift  and  drag  forces  are  resolved  into  force  and  moment 
components  in  the  ship  coordinate  system. 

G.  Other  Forces 

An  additional  damping  force  may  be  computed  for  any  of  the 
six  motion  components.  Let  f.  be  one  of  the  components  of  the 
force  or  moment  in  the  ship  coordinate  system,  and  let  v.  be 

X’ 

the  corresponding  component  of  linear  or  angular  velocity.  The 
force  or  moment  is  computed  from 

fi  - fi'  - \Li  - I Of  Iff  (V-31) 

where 

f * 

J i * force  or  moment  due  to  the  waves  and  steering  system, 

L.  = coefficient  of  linear  damping  for  the  particular 
motion  component, 
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and 

Q.  * coefficient  of  quadratic  damping  for  the  particular 

"Lf 

motion  component. 

There  are  three  force  and  three  moment  equations  using  these 
three-dimensional  linear  and  quadratic  damping  constants,  but 
in  practice  the  only  non-zero  coefficients  are  for  surge  and 
roll. 

Surge  damping  may  be  controlled  by  equation  (V-31)  or  by 
an  optional  table  of  resistance  versus  speed  data.  When  this 
option  is  used,  a surge  damping  (or  accelerating)  force  is 
obtained  by  interpolation  in  the  resistance  table.  This  surge 
force  is  equal  to  the  resistance  at  the  desired  speed  minus 
the  resistance  for  the  velocity  predicted  at  each  time  step. 

H.  Examples  of  Numerical  Simulations 

Several  examples  of  motions  simulated  by  the  numerical 
time  domain  integration  are  presented  here.  These  computations 
were  performed  for  the  American  Challenger  class  of  cargo  ship 
for  which  a model  was  previously  tested  on  San  Francisco  Bay. 
Details  of  the  experiments  with  this  vessel  are  given  in  the 
reports  by  Haddara,  et  al.  (1971)  and  Paulling,  et  al.  (1972). 
Characteristics  of  the  vessel  and  its  displacement  condition 
for  the  simulations  are  given  in  Table  V-I.  Table  A-I  in 
Appendix  A contains  a detailed  compilation  of  the  vessel's 
characteristics.  A bodyplan  is  also  included  in  that  appendix 
as  Figure  A-l.  The  offsets  for  the  simulations  were  limited 
to  13  stations.  Motion  calculations  were  carried  out  at  ship 
scale  rather  than  at  the  scale  of  the  model  used  for  the  ex- 
periments. The  ratio  of  ship  length  to  model  length  is  30.189, 
and  the  time  scale  for  Froude  number  similarity  is  about  5.5. 
Metacentric  heights  corresponding  to  model  ballast  weight 
positions  2 and  3 were  used  for  the  simulations.  The  righting 
arm  curves  for  the  Challenger  are  shown  in  Figure  A- 3 in 
Appendix  A.  In  addition  to  the  curve  for  still  water,  the  righting 
arms  for  the  ship  in  a longitudinal  wave  with  crest  amidship  and 
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trough  amidship  are  plotted. 

Three  simulations  with  a single  sinusoidal  wave  are  shown 
in  Figures  V-2  through  V-7.  The  wave  is  overtaking  the  ship 
from  dead  astern.  The  circular  wave  frequency  is  0.575  radians 
per  second,  corresponding  to  a wave  length  of  610  feet. 

The  first  example  is  for  a wave  amplitude  of  6.5  feet. 

The  initial  roll  angle  is  15  degrees  to  starboard.  To  de- 
crease the  effect  of  starting  transients  caused  by  the  choice 
of  initial  conditions,  all  forces  applied  to  the  ship  are 
multiplied  by  a ramp  function  which  increases  linearly  with 
time  from  zero  to  one  for  the  first  10  seconds  of  the  integration. 
The  roll  history  is  plotted  as  Figure  V-2,  and  that  for  pitch 
is  Figure  V-3.  The  pitch  record  may  be  used  to  judge  the  number 
of  wave  encounters.  In  this  example  the  rolling  occurs  as  one- 
half  of  the  wave  encounter  frequency,  and  the  roll  amplitude 
builds  up  rapidly.  The  first  roll  to  port  has  an  amplitude  of 
about  19  degrees,  and  the  last  starboard  roll  is  to  an  amplitude 
of  28  degrees.  The  ship  capsizes  to  port. 

The  next  example  is  for  the  same  wave  condition,  but  the 
initial  roll  angle  is  5 degrees  to  starboard.  The  roll  record 
is  Figure  V-4,  and  pitch  is  Figure  V-5.  Rolling  again  occurs  at 
one  half  the  pitch  and  wave  encounter  frequency.  The  roll  motion 
peaks  at  eleven  degrees  after  130  seconds,  and  it  appears  to 
stabilize  thereafter. 

The  third  example  repeats  the  initial  conditions  of  the 
previous  one  with  the  five  degree  initial  roll,  but  the  wave 
amplitude  is  increased  to  7.5  feet.  The  roll  and  pitch  responses 
are  Figures  V-6  and  V-7.  The  interesting  thing  here  is  that  the 
resulting  roll  amplitude  is  smaller  for  this  case  with  a larger 
wave  amplitude.  The  roll  amplitude  for  this  run  is  about  eight 
degrees . 

The  remaining  examples  are  for  a wave  system  consisting  of 
two  sinusoidal  waves  traveling  in  the  same  direction.  The 
circular  frequencies  of  the  waves  are  0.433  and  0.601  radians 
per  second,  which  correspond  to  wave  lengths  of  1076  and  558 
feet.  The  amplitudes  of  both  wave  components  were  equal  ? n 
magnitude.  Three  simulation  runs  for  following  seas  and  this 
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two- component  sea  are  shown  in  Figures  V-8  through  V-10,  and 
V-12  through  V-17.  The  metacentric  height  for  these  runs  is 
0.558  feet  (GM  position  3). 

To  initiate  rolling,  the  waves  approach  the  ship  from  the 
direction  of  ten  degrees  to  starboard  of  dead  astern.  The 
amplitude  of  each  component  wave  is  five  feet  making  a twenty 
foot  maximum  wave  height  in  the  group.  The  first  simulation 
begins  with  the  two  wave  components  in  phase  at  the  position 
of  the  ships  center  of  gravity.  The  starting  ramp  is  applied 
to  the  first  60  seconds  of  the  run.  Figure  V-8  shows  the  wave 
amplitude  at  the  position  of  the  center  of  gravity  of  the  ship. 
The  roll  record  is  displayed  in  Figure  V-9.  After  about  120 
seconds  the  rolling  begins  to  take  place  at  one  half  of  the 
frequency  of  wave  encounter.  As  the  wave  amplitude  builds, 
the  rolling  becomes  erratic:  rolling  alternates  with  the 

encounter  frequency  and  one  half  of  that  frequency.  Figure 
V-10  shows  the  pitch  record  for  this  run.  The  ramp  function 
used  for  starting  the  simulation  results  in  the  attenuation 
of  the  pitch  at  the  beginning  of  the  record.  Note  that  the 
pitch  angle  gives  a good  indication  of  the  wave  slope  for  these 
waves  which  are  longer  than  ship  length.  The  pitch  and  roll 
experiment  records  for  the  SL-7  model  shown  in  Figure  V-ll 
exhibit  this  same  erratic  rolling  with  the  roll  motion  occurring 
at  the  wave  encounter  frequency  and  then  at  one  half  of  this 
frequency. 

Figures  V-12  through  V-14  show  the  wave,  roll,  and  pitch 
records  for  a simulation  with  the  same  initial  conditions 
as  before,  except  that  the  phase  angles  of  the  two  waves  are 
initially  one  half  cycle  apart  causing  the  amplitudes  to  cancel. 
The  starting  ramp  is  reduced  to  30  seconds.  In  this  run  the 
rolling  is  at  the  same  frequency  as  the  wave  encounter,  and  the 
amplitude  is  less  than  half  that  of  the  previous  example. 

The  last  example,  for  the  two-component  following  seas, 
shows  the  wave,  roll  and  pitch  records  in  Figures  V-15  through 
V-17.  The  initial  conditions  are  the  same  as  for  the  previous 
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example  except  that  an  initial  roll  angle  of  fifteen  degrees 
to  starboard  has  been  added.  Rolling  begins  at  one  half  the 
encounter  frequency  which  is  close  to  the  natural  roll  fre- 
quency for  calm  water  and  this  roll  amplitude.  As  the  wave 
amplitude  and  roll  exciting  moment  build  up  the  ship  is  forced 
to  roll  at  the  encounter  frequency. 

Figure  V-18  shows  a Challenger  model  experiment  in  following 
seas  for  the  corresponding  speed  and  GM  condition.  Although 
the  experiment  ended  in  a capsize  with  low  cycle  resonance 
and  these  simulations  did  not,  both  the  experiment  and  the 
first  and  last  following  sea  calculations  show  a tendency  for  the 
roll  frequency  to  switch  back  and  forth  between  the  encounter 
frequency  and  one  half  the  encounter  frequency. 

Two  examples  of  capsize  simulations  in  quartering  seas 
are  shown  in  Figures  V-19  through  V-25.  The  two-component 
wave  system  that  was  used  for  the  following  sea  runs  is  used 
here,  but  the  direction  is  such  that  the  waves  approach  the 
ship  from  35  degrees  to  starboard  of  dead  astern.  At  the 
beginning  of  the  runs  the  two  wave  components  are  in  phase. 

The  first  quartering  example  is  for  a metacentric  height 
of  0.257  ft  (GM  position  2).  The  wave  record  for  this  run 
is  shown  in  Figure  V-19.  The  starting  ramp  was  60  seconds  in 
length,  and  ended  near  the  end  of  the  first  wave  group.  Roll 
and  pitch  records  are  in  Figures  V-20  and  V-21.  At  about  forty 
seconds  from  the  start  a wave  crest  comes  amidship  and  vessel 
rolls  j 4 degrees  to  port.  The  next  wave  crest  comes  amidships 
at  about  sixty  seconds  and  a starboard  roll  into  the  wave  of 
4 degrees  is  reached.  The  next  roll  is  22  degrees  to  port  just 
after  the  next  crest  passes.  The  ship  rolls  to  starboard  as 
the  waves  begin  to  build  in  the  second  wave  group.  The  maximum 
starbf  ard  roll  of  23  degrees  is  reached  as  the  next  crest  moves 
away.  The  roll  momentum  imparted  the  ship  by  the  increased 
righting  arm  of  the  next  trough  coming  amidship  together  with 
the  reduced  stability  in  the  crest  that  follows  causes  the  ship 
to  capsize  to  port — clearly  an  example  of  low  cycle  resonance. 
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Figure  V-22  shows  the  last  two  and  a half  minutes  of  a model 
run  corresponding  to  the  same  speed,  heading  relative  to  the 
waves  and  GM  condition.  The  experiment  also  ended  with  a 
capsizing  to  port  caused  by  low  cycle  resonance. 

An  example  of  pure  loss  of  stability  is  shown  in  Figures 
V-23  through  V-25.  The  metacentric  height  was  0.557  ft. 

(GM  position  3) . The  wave  amplitude  at  the  center  of  gravity 
of  the  ship  is  shown  in  Figure  V-23.  Roll  and  pitch  are  in 
Figures  V-24  and  V-25,  Rolling  is  at  one  half  the  frequency 
of  wave  encounter  between  40  and  110  seconds  into  the  run. 
After  that  the  wave  amplitude  and  the  roll  moment  are  so 
large  that  the  ship  is  forced  to  roll  at  the  encounter  fre- 
quency. The  vessel  finally  capsizes  with  the  combination 
of  a large  port  roll  and  a wave  crest  amidships. 
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TABLE  V - I 

CHARACTERISTICS  OF  THE  AMERICAN  CHALLENGER 
USED  FOR  THE  TIME  DOMAIN  SIMULATION 


Length  Between  Perpendiculars 

Characteristic  length 

Beam 

Depth3 

Mean  Draft 

Displacement 

Speed 

Froude  number 

Meuacentric  heights: 

GM  position  2 

GM  position  3 


SHIP  MODEL 


529.0  ft 

17.523 

ft 

521.0  ft 

17.258 

ft 

75.0  ft 

2.484 

ft 

46.5  ft 

18.484 

ft 

29.75  ft 

11.827 

in 

19652.5  tons 

1600  lbs 

27.7  ft/sec 

2.99 

ft/sec 

0.215 

0.215 

0.257  ft 

0.102 

in. 

0.558  ft 

0.222 

in. 

aThis  depth  is  greater  than  the  depth  of  the  actual  ship 
because  the  model  was  constructed  up  to  the  level  of  the  top 
of  the  bulwark  rather  than  to  the  main  deck  level  of  the  ship. 
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Figure  V- 3 Pitch  Simulation 
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Figure  V-4  Roll  Simulation 


Figure  V-5  Pitch  Simulation 


Tine  ( SECONDS 


Figure  V-7  Pitch  Simulation 


Figure  V-8  Wave  Simulation 


Tint  C SCCONOS 


I 


Figure  v-10  Pitch  Simulation 
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Figure  V-ll  SL-?  Model  Experiment  Records 
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Figure  V-13  Roll  Simulation 
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Figure  V-15  Wave  Simulation 


Figure  V-16  Roll  Simulation 


Figure  V-17  Pitch  Simulation 
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Figure  V-18  American  Challenger  Model  Experiment 


Figure  V-19  Wave  Simulation 
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Figure  V-20  Roll  Simulation 
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Figure  V-23  Wave  Simulation 


Lyure  V-24  Roll  Simulation 


I 


Figure  V-25  Pitch  Simulation 
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VI  ENCOUNTER  PHENOMENA,  WAVE  GROUPS,  AND  PREDICTION 

It  is  obvious  that  intact  ship  capsizing  is  due,  in  large 
part,  to  the  configuration  of  the  waves  encountered  by  the  ship. 
There  is  a great  deal  of  difference,  however,  between  the  wave 
shape  observed  in  a fixed  coordinate  system  and  the  resulting 
excitation  experienced  by  the  moving  vessel  in  following  seas. 
The  following  discussion  is  intended  as  a review  of  some  simple 
wave  encounter  relations  and  an  exploration  of  the  potential 
usefulness  of  wave  group  statistics. 


A.  Wave  Encounter  Phenomena 


A number  of  striking  phenomena  are  observed  in  the  wave 
pattern  encountered  by  a vessel  moving  in  following  seas.  None 
of  these  observations  are  new,  yet  some  discussions  of  waves  in 
the  encounter  domain  is  in  order  in  light  of  the  capsize  modes 
discussed  in  previous  chapters. 

The  most  dramatic  event  in  the  encounter  domain  occurs  at 
zero  frequency.  A vessel  moving  at  or  near  the  velocity  of  a 
wave  crest  may  experience  a significant  non-oscillating  force 
and  a critical  loss  of  stability.  Capsizing  due  to  broaching  or 
to  a pure  loss  of  stability  can  result.  The  zero  occurs  at  an 
absolute  wave  frequency  of  g/  ( UcosX)  where  U cob  X is  the 
ship  velocity  in  the  direction  of  wave  travel.  A singularity 
appears  at  this  frequency  in  the  prediction  of  the  unrestored 
motions  using  linear  theory.  Normally  few  displacement  vessels 
move  at  such  large  Froude  numbers.  Exceptions  might  be  forth- 
coming as  the  speed  of  containerships  increases  and  fast  vessels 
are  driven  at  their  limit.  Small  craft,  such  as  fishing  boats 
and  tugs,  can  experience  problems  when  they  operate  in  shallow 
water.  The  phase  speed  of  wave  travel  can  be  reduced  to  suffi- 
ciently small  values  that  such  craft  have  no  difficulty  in 
running  with  a crest  amidships.  For  example,  a modern  tuna 


seiner  144  feet  on  the  waterline  with  a design  speed  of  approx- 
imately 15  knots  travels  at  the  phase  speed  of  a 125-foot  wave. 

In  thirty  feet  of  water,  a deep  water  wave  of  160  feet  will  have 
a length  of  approximately  140  feet,  it  will  be  steeper,  and  will 
run  at  this  same  ship  speed. 

The  folding  of  the  wave  spectrum  in  the  encounter  domain  is 
also  familiar.  Figure  VI-1  shows  a typical  experimental  wave 
spectrum  in  the  absolute  and  in  the  frequency  of  encounter  domain 
for  two  speeds.  The  collapse  of  virtually  all  of  the  energy  to 
a narrow  band  of  rather  low  frequencies  should  be  noted  and  is 
potentially  significant  for  purposes  of  mathematical  simulation. 

I 

; The  frequency  collapse  with  respect  to  wave  length  is  dramatically 

shown  in  Figure  VI-2  presented  by  Kastner  (1973) . It  is  evident 
that,  for  a given  ship  speed  V , a quartering  situation  may  involve 
the  most  serious  concentration  of  energy.  The  situation  could  be 
serious  if  this  frequency  is  near  the  natural  frequency  of  rolling 
or  in  the  neighborhood  of  the  frequency  for  low  cycle  resonance. 

It  has  often  been  stated  that  the  seaway  appears  to  be  a confused 
collection  of  waves  with  random  phasing  and  the  only  thing  regular 
about  the  sea  is  its  irregularity.  Intuition  would  therefore 
, suggest  that  extended  excitation  at  a single  frequency  or  a very 

narrow  frequency  band  is  unlikely,  especially  when  the  spectrum 
tends  to  be  rather  broad  banded.  It  is  apparent,  however,  that  in 
following  seas  the  picture  is  radically  altered. 

Excitation  by  a series  of  large  waves  at  or  near  a single 
frequency  appears  to  be  a requirement  for  capsizing  due  to  low 
cycle  resonance.  These  conditions  may  be  a far  more  common  event 
than  might  first  be  assumed.  For  purposes  of  discussion,  it  will 
be  convenient  to  consider  a wave  group  or  packet  formed  by  the 
superposition  of  two  regular  wave  components: 

i.  n (x,t)  = (k.x  - u.t) 

+ cos  ( k2x  - w t) } . 
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Figure  VI-2  Encounter  Frequency  vs  Heading 
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The  number  of  waves  in  a group  at  any  instant  of  time  is  given  by 


k + k 

N[t=o]= — 1 *- 

2(kx  - k ) 


and  is  the  number  that  would  be  observed  from  a photograph  or  film. 
This  is  not  the  same  group  size  obtained  from  a record  made  at  a 
fixed  point,  howevar.  A wave  probe  located  at  x=0  would  observe 


GO  + to 

N [x=0]  1 2 

2 (to  j - u>2  ) 


waves  in  each  group.  Figure  VI-3  (a)  is  a plot  of  the  sum  of  two 
waves  of  frequency  0.433  and  0.601  rad/sec.  The  number  of  waves  in 
each  group  therefore  represents  an  increase  by  a factor  of  1.9 
over  the  number  that  would  be  observed  at  a fixed  instant.  This 
discrepency  is  often  overlooked  when  comparing  wave  records  with 
motion  records. 

A dramatic  change  occurs  with  forward  speed  in  following  seas. 
The  number  of  waves  encountered  by  a vessel  moving  at  a speed  U cos  X 
is 

(k  +k  ) U CO  3 Y - (to  "f  to  ) 

12  12 

N [x=Utcosx]=  * 

2{(k  )UcosX-  (to  -to 2 ) } 


Here  the  observed  wave  frequency  is 


f e = {(k  x+k  2)Uco8X-  (Wj+oOj  ) }/(4n) 


which  is  zero  when 


I'cosx 


Gu  t GO 
1 2 


(W 


As  noted  earlier,  this  is  often  a very  high  speed  for  more  dangerous 
waves  of  sufficient  length  and  amplitude.  The  ratio  of  the  number 
of  waves  encountered  at  forward  speed  to  the  number  observed  at  any 
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Figure  VI- 3b  Simulated  Wave  Records 
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one  instant  is  therefore 

' at  + W ) 

1 ? 

s>. 


I'-'O?  x-  

k3) 

Figure  VI-3  (b)  shows  the  sane  set  of  waves  pictured  above  but 
viewed  in  the  encounter  domain  where  l'^cay=211.?  fps.  There  is 
a gain  by  a factor  of  4.5  in  the  number  of  waves  in  each  group 
over  an  observation  at  fixed  time  and  a gain  of  3.1  over  an  ob- 
servation at  a fixed  point.  These  facts  are  well  known  but  need 
to  be  repeated  to  be  appreciated.  If  the  ship  moves  at  or  near 
the  group  velocity  it  will  encounter  a very  large  number  of  waves 
even  when  the  observed  wave  aroup  or  packet  appears  to  be  rather 
small.  This  also  implies  that  the  frequency  of  encounter  could 
be  very  pure  during  such  a wave  group.  If  the  wave  group  con- 
tains one  or  more  large  waves,  their  repeated  encounter  in  the 
moving  domain  could  become  a favorable  situation  for  low  cycle 
resonance  or  a combination  of  resonance  and  large  stability  loss. 

B.  Wave  Group  Statistics 

The  experiments  described  in  Chapter  III  provided  ample 
evidence  of  the  importance  of  "wave  groups"  in  the  process  of 
broaching  and  capsizing.  Further  insight  was  gained  by  the  nu- 
merical simulation  and  analysis  of  Chapter  V.  It  was  also  ob- 
vious that  numerical  simulations  of  long  runs,  though  faster  and 
less  costly  than  open  water  experiments,  are  also  very  expensive. 
Yet  some  sort  of  time  domain  simulation  is  required  due  to  the 
nonlinear  nature  of  the  phenomena.  The  statistical  results  from 
many  long  runs  would  probably  provide  the  most  accurate  data  base 
for  extrapolation.  The  success  of  such  numerical  simulation  de- 
pends upon  the  knowledge  of  the  wave  system  actually  encountered  by 
the  vessel.  Running  a finite  number  of  capsize  simulations  for 
rather  specific  wave  configurations  appear  to  be  the  only  practical 
approach.  This  of  course  does  not  permit  a complete  coverage  of 
all  relevant  wave  configurations  which  means  that  a statistical 
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analysis  based  solely  on  the  runs  executed  is  not  possible 
( Paulling,  e_t  al^  1972).  However,  the  complexity  of  the 
wave  system  and  the  apparent  irregularity  suggest  that  only  a 
probabilistic  description  of  the  seaway  and  approach  to  ex- 
treme motion  prediction  will  be  successful.  Attempting  to 
assign  probabilities  to  any  but  the  most  simple  wave  events 
is  presently  out  of  the  question.  Until  seastate  monitoring 
techniques  can  provide  Su . f iciently  detailed  and  accurate 
information  of  extreme  values,  the  possibility  of  using  a 
simplified  description  of  wave  events  leadinq  to  extreme  motions 
or  capsizing  should  not  be  overlooked.  This  is  where  the 
concept  of  wave  groups  may  permit  engineering  estimates  of  the 
probability  of  capsizing  and  a measure  of  ship  safety. 

Ideally  it  would  be  possible  to  separate  the  study  of  cap- 
sizing from  the  occurrence  of  relevant  wave  configurations.  For 
example,  it  may  be  possible  to  exercise  the  simulator  to  deter- 
mine those  wave  group  parameters  that  result  in  a capsize  or 
nearly  so.  This  might  be  referred  to  as  the  "capsize  surface". 
The  probability  of  capsiz  g would  then  be  expressible  as  the 
product  of  the  probability  of  capsizing  in  a given  wave  group 
(zero  or  one)  and  the  probability  of  encountering  such  a wave 
group.  The  required  contributions  of  the  naval  architect  and 
the  oceanographer  thereby  become  somewhat  more  clearly  defined. 

Two  questions  must  be  answered  if  such  an  approach  is 
worth  pursuing.  The  first  is  whether  a relatively  smooth  cap- 
size surface  exists  and  is  determinable.  The  second  question 
is  whether  the  probability  of  occurence  of  any  given  wave  group 
configuration  can  be  estimated. 

Chapter  V presented  a short  analysis  of  the  sensitivity  of 
capsizing  to  changes  in  wave  configuration.  Though  by  no  means 
exhaustive,  the  first  indications  of  the  study  were  that  moderate 
changes  in  the  input  wave  group  resulted  in  smooth  chances  in 
the  motions. 

The  second  question  involved  the  determination  of  the 
probability  of  encountering  a given  wave  group.  Theoretical  and 
experimental  studies  of  wave  group  statistics  to  date  are  not 
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sufficiently  complete  for  practical  use.  However,  the  following 
review  and  analysis  suggests  that  practical  results  may  be 
possible  and  forthcoming. 

The  exact  definition  of  a wave  group  has  yet  to  be  given. 

The  definition  is  in  fact  somewhat  arbitrary  and  depends  on  its 
uses.  Most  investigators  have  employed  the  most  convenient  form 
for  their  own  purposes.  Each  definition  aims  to  characterize 
the  beat-like  appearance  of  the  wave  amplitude  which  is  often 
visible  in  an  irregular  wave  train.  Low  frequency  components 
and  asymmetry  are  usually  removed  prior  to  statistical  analysis 
by  using  the  wave  height  or  the  envelope  of  the  process  as 
the  basic  parameter.  A group  is  usually  defined  by  employing 
one  or  more  of  the  following  parameters: 

Hg  the  largest  wave  height  in  a sequence 

Hl  the  level  which  a sequence  of  N waves 

exceed 

N the  number  of  waves  in  the  group 

Tg  the  average  period  of  the  component  waves 

T.  the  period  of  the  envelope  exceeding  the 
level  Hl 

T2  the  period  between  up-crossings  of  level 
Hl  by  the  envelope 

H the  height  measured  between  a crest  and 
the  average  of  neighboring  troughs 

A group  may  be  defined  as  a series  of  waves  for  which  the  wave 
height  envelope  exceeds  a specified  level.  An  alternative 
definition  might  involve  the  period  of  each  individual  wave 
component.  The  group  defined  by  Thompson  (1974)  was  based  on 
the  height  of  the  largest  wave  in  a sequence  where  the  approxi- 
mate period  of  each  wave  component  was  within  a given  range  of 
the  average  period  for  the  group.  The  largest  height,  H^,  had 
to  be  greater  than  or  equal  to  one  third  the  significant  wave 
height. 
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Figure  VI-4  Wave  Group  Definition  Sketch 

Other  definitions  can  be  easily  constructed.  The  largest  wave 
in  a group  is  likely  to  be  the  most  important  parameter  for  pre- 
diction of  capsizing  due  to  pure  loss  of  stability.  On  the  other 
hand,  the  average  period,  overall  level  HL,  and  length  of  run  T, 
should  be  more  relevant  for  studies  of  resonant  buildup  of  motions. 
Additional  experience  with  numerical  simulation  should  help  to 
clarify  the  important  parameters. 

Once  the  parameters  are  chosen,  it  is  necessary  to  construct 
the  joint  probability  distribution,  e.g.  P (H,  Hg,  Tg,  N)  . This 
can  be  done  by  fitting  curves  to  the  actual  data  records.  Some 
parameters  have  been  studied  in  this  way  (Thompson,  1974). 

Assistance  from  theory  is  also  available,  but  there  is  not  enough, 
as  yet,  to  construct  the  entire  joint  probability  for  typical 
wave  spectra.  A detailed  study  of  wave  groups  is  necessary.  As 
a first  step,  a brief  review  of  some  of  the  theoretical  highlights 
and  the  correlation  with  data  recorded  in  San  Francisco  Bay 
will  be  given. 
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Several  of  the  references  listed  at  the  end  of  this  report 
contain  discussions  of  certain  aspects  of  "wave  groups".  Among 
them,  two  different  approaches  were  used  for  the  analysis.  For 
a given  level,  the  first  approach  considers  the  number  of  waves 
in  a group  while  the  second  approach  considers  the  elapsed  tire 
of  the  envelope  above  the  level. 

Thte  first  approach  is  represented  by  Goda  (1970)  and  Ewing 
(1972  and  1973) . Ewing  dealt  with  the  mean  number  of  waves  in 
or  between  groups  of  waves  in  a record  for  a given  level.  Mean 
value  statistics  however,  are  not  likely  to  be  sufficient  for 
the  above  purposes. 

C.  The  Probability  of  Having  a Given  Number  of  Waves  in  a Group 

In  Goda's  (1970)  analysis  of  "wave  groups",  he  defined  the 
length,  j , of  a run  of  height,  h , as  the  wave  heights  exceedinc 
h for  j-1  consecutive  waves  after  one  exceedence  and  then  a 
failure  to  exceed  the  level  at  the  j+1  wave.  Let  the  probability 
that  each  wave  height  will  exceed  h be  p,  then  the  probability 
that  the  wave  will  not  exceed  h is 

q = 1 - p.  (VI-1) 

If  the  wave  heights  are  independent  of  each  other,  then  the 
probability  of  having  j i consecutive  waves  with  a height  greater 
than  h after  the  first  up-crossing  of  level  h is 

P (j  ) = P j’-1  q (VI-2) 

1 ^ • 


The  mean  and  variance  of  the  length  of  runs  in  a record  are: 


E(j  ) = ± 
i <3 


(VI-3) 


0 (j  > = It 

i <3 


(VI-4) 
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Next,  he  defined  the  recurrence  interval  for  wave  groups 
as  a total  run  at  level  h , which  is  the  number  of  waves  between 
two  successive  groups.  The  probability  of  a total  run  with  a 
length  of  j is  given  by 

pq  -i  -1  i-1 

P2  ( j 2 ) = (pD*  - q]2  ).  (VI-5) 

p-q 


The  mean  and  variance  of  the  length  of  a total  run  are 
calculated  as 


E ( j ) 
2 


a 2 ( j 2 ) 


(VI-6) 


(VI- 7) 


For  comparison  of  this  theory  with  his  simulated  wave  record, 
three  definitions  of  runs  were  used:  waves  exceeding  the  median 

wave  height  11rnedian  / waves  exceeding  the  significant  wave  height 
H i ^ , and  the  total  run  length  from  the  two  successive  wave  groups 

exceeding  Hi^,  . His  results  are  shown  on  Fiqure  VI-5.  The  spectrum 

M-9  has  almost  uniform  spectral  density  and  exhibits  run  lengths 
that  are  similar  to  those  predicted  by  this  independent  random 
process.  The  other  extreme,  namely  a narrow  banded,  highly  peaked 
spectrum,  is  represented  by  spectrum  M-5.  The  distribution  of  the 
length  of  runs  differs  significantly  from  that  represented  by  this 
process.  On  the  other  hand,  a more  realistic  wave  spectrum  is 
represented  by  spectrum  F-4  which  is  of  Pierson-Moskowitz  or  wind 
wave  form.  The  distribution  characteristics  of  the  length  of  runs 
of  this  spectrum  follows,  in  general,  the  shape  of  that  represented 
by  this  independent  wave  process,  but  the  magnitude  is  so  different 
that  this  process  cannot  describe  it  accurately. 

The  same  analysis  has  been  performed  for  aves  measured  in 
San  Francisco  Bay.  A typical  result  for  a wave  record  at  various 
heiqhts  is  shown  on  Figure  VI-6  and  Figure  VI-7.  These  resemble 
the  results  shown  in  Figures  VI-5  from  spectrum  F-4. 
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Figure  Vi-5  Length  or  Kun  vs.  Wave  Height 


Note  that  the  agreement  between  wave  data  and  the  results 
of  predictions  by  this  method  is  better  at  lower  levels  than  at 
higher  levels.  The  latter,  however,  are  of  more  interest.  Also, 
note  that  for  wind  waves,  as  represented  by  the  San  Francisco 
Bay  wave  data  or  spectrum  F-4 , the  probability  of  only  one  wave 
in  a group  is  always  smaller  than  that  predicted  by  Goda's 
method  while  the  probability  that  there  are  more  waves  in  a group 
is  always  larger.  Therefore  it  can  be  concluded  that: 


1.  larger  waves  tended  to  come  in  groups, 

2.  successive  whve  heights  are  correlated  with  each  other, 
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3.  the  distribution  of  the  length  of  runs  described  by 

this  independent  random  process  can  represent  only  the 
shape  of  the  curve. 

D.  The  Probability  of  the  Elapsed  Time  of  the  Envelope  Being 
Above  a Given  Level. 

This  approach  is  performed  by  Nolte  and  Hsu  (1973),  (to 
be  published) , and  Vanmarcke  (1971) . Again,  Vanmarcke  discussed 
mainly  the  mean  values  of  various  parameters.  Although  useful, 
his  analysis  will  not  be  pursued  further  here. 

Nolte  and  Hsu  (1973)  assumed  that  the  probability  that  the 
envelope  height  will  cross  a level  hc  during  an  interval  of  time 
can  be  described  by  a Poisson  process.  The  assumptions  associated 
with  using  the  Poisson  model  are  stated  by  them  as  follows: 

1.  The  probability  of  the  envelope's  crossing  the  level 
h0  in  any  non  overlapping  time  interval  is  independent 
and  time  invariant. 

2.  The  probability  of  crossing  the  level  hD  in  a small 
time  interval  is  proportional  to  the  size  of  the 
interval . 

3.  The  probability  of  multiple  crossings  in  a small  time 
interval  is  negligible  relative  to  the  probability  of  a 
single  crossing. 

They  then  define  the  time  duration  for  a group  of  level  h0 
as  the  elapsed  time  beginning  when  the  envelope  crosses  above  the 
level,  h0  , and  ending  when  the  envelope  crosses  below  the  level 
h0  • This  time  duration  is  shown  on  Figure  VI-4  as  t l . The 
probability  that  the  time  duration  is  greater  than  t is  expressed 
by  them  as  an  exponential 


P (t  j>t|h0 ) = exp [-t/x (h„ ) ] 


(VI-8) 


P(N) 


L26 


Goda’s  "Random  Process" 

Data  of  "Seaway  14  Aug  73" 
Sequence  1 probe  1 
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where  t (h0)  is  the  average  time  duration  for  groups  of  level  h0. 
They  derived  this  mean  time  duration  as 


t (h0) 


1 


2 


(VI-9) 


in  which  h0 
significant 
defined  by 


is  the  level  of  the  envelope  height  divided  by  the 
wave  height.  f is  the  average  wave  frequency 
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J%  S ( f ) f 2 df 
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(VI-10) 


S(f)  is  the  wave  energy  spectrum,  f is  the  centroid  of  the 
wave  energy  spectrum  defined  as 


f 


S(f)f  df 


S(f)  df 


(VI-11) 


th 

and  m , is  the  i moment  of  the  energy  spectrum. 
i 

With  these  known,  they  further  derive  the  average  frequency 
of  occurrence  for  wave  groups  of  level  h as 


v (h  ) 
0 


exp(-2hp) 

t (hQ ) 


(VI-12) 


If  the  number  of  waves  in  a group  is  denoted  by  N,  then 
the  probability  that  a wave  group  contains  more  than  n waves 
higher  than  the  level  h0  is 


P(N>n|h0)  = exp  1-nT/x  (h0  ) ] 


(VI-13) 


129 


where  T is  the  average  wave  period. 

Since  the  number  of  groups  in  a record  of  time  L is 
assumed  to  be  described  by  a Poisson  process,  the  probability  of 
having  M groups  of  waves  of  level  h0  can  be  expressed  as 

[v(h0)L]M 

P (M | h0 ) = exp [-v (h0 ) L] . (VI-14) 

M! 

If  the  highest  wave  height  in  a group  of  N waves  exceeding 
the  level  h0  is  denoted  by  H1,  then  the  probability  that  H1 
is  less  than  h is  derived  as 

P (H 1 <h | h0, N)  = l-exp[-2 (h2-hl) ]N.  (VI-15) 


In  normalized  form,  the  energy  spectra  for  the  envelope 
and  the  square  of  the  envelope  are  shown  to  be  equal  and  can 
be  expressed  as 


E(f) 


rs  (x)  S (f+x)dx 

- oo  0 

f S2 (x)dx 

J o 


(VI-16) 


where  S(x)  is  the  energy  spectrum  for  the  wave  trace  and  x 
is  a dummy  variable  representing  frequency. 

The  assumption  of  a Poisson  model  is  shown  by  Nolte  and  Hsu  to 
agree  very  well  with  data  collected  during  one  storm  in  the  Gulf 
of  Mexico.  Similar  statistics  have  been  determined  for  the 
waves  measured  in  San  Francisco  Bay.  Some  results  are  shown 
in  Figures  VI-8,  VI-9,  and  VI-10  for  different  levels  of  the 
same  wave  record. 

In  comparing  their  results  with  the  data  they  collected, 

Nolte  and  Hsu  noted  that  because  of  the  finite  wave  population 
of  the  record,  x(h0)  obtained  from  Equation  (VI-9)  cannot 
represent  the  mean  time  duration  properly. 
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Figure  VI-8 
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Figure  VI-9 


DISTRIBUTION  OF  THE  AVERAGE  PERIOD 


SEAWAY-2  8 AUG  73,  SEQ.  3,  PROBE- 1 
HEIGHT  AT  LEVEL  1.2173 


KOLMOGOROV  STATISTICS 

- EXPONENTIAL  [t  (h  )]:  0.4783 

- Exponential  [T]  : 0.1826 

- Gamma:  0.1522 

Critical  Value 

(of  level  0.05) : 0.2750 

(of  level  0.20) : 0.2160 


Figure  VI-10 


(This  fact  was  clearly  shown  by  the  San  Francisco  Bay  Wave  data 
Further  study  showed  that 


uo 

Jh  P (h)  dh  - -1 

T (ho  ) = 5 ( h ) 


(IV-17) 
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— /nOO 

h = / h P (li) dh 
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Since  h is  a constant  for  a given  wave  record,  then 


r(h0)  <* 


P (h)  dh 


P(h„) 


( IV-19 ) 


This  relation  may  then  be  used  to  obtain  the  mean  tine  duration 
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x(h0)  . It  is  not  obvious  that  we  can  ignore  h and  obtain  the 
same  statistics.  On  the  other  hand,  the  main  concern  at  this 
point  is  to  find  out  the  statistical  characteristics  of  x , 
rather  than  the  actual  mean  time  duration.  The  actual  data 
mean  is  therefore  used  in  the  plots.  The  results  seem  to  follow 
the  exponential  curve  reasonably  well.  Kolmogorov  statistics 
have  been  calculated  to  evaluate  the  quality  of  fit  and  have 
been  shown  to  accept  the  hypothesis  that  the  data  is  exponentially 
distributed  at  a level  of  0.2  . 

A Poisson  process  is  a very  handy  assumption  to  make  in 

many  situations.  It  is  formed  by  the  so  called  streams  of 

random  events  with  a continuous  time  parameter  (Cramer  and 

Leadbetter,  1967)  . In  most  applications  we  are  concerned  with 

events,  such  as  calls  on  a telephone  line  or  arrivals  of 

customers  in  a shop,  which  occur  successively  in  time.  A Poisson 

process  is  usually  assumed  for  these  events  because  in  some 

cases  the  Poisson  model  will  give  an  approximation  of  the  true 

situation  and  in  others  it  will  at  least  provide  a starting 

point  for  various  generalizations  (Cramer  & Leadbetter.-  1967). 

The  time  interval  between  events  for  a Poisson  process  is  an 

exponentially  distributed  random  variable  as  was  implied  by 

equation  (VI-8) . In  the  time  series  for  a given  level  hD  » 

if  we  consider  any  crossing  of  the  envelope  as  an  "event" 

then  we  should  count  x ' as  well  as  x in  Figure  VI-4  as 

i i 

the  time  intervals  between  events.  On  the  other  hand  if  we 
consider  only  up-crossings  as  "events",  then  it  would  be  more 
appropriate  to  consider  x^  as  the  time  intervals  between 

events.  In  this  case  the  number  of  events  in  a time  interval  is 
the  same  as  the  number  of  wave  groups  in  that  interval.  The 

distribution  curve  for  x has  been  calculated  and  the  results 

2 

shown  in  Figures  VI-11,  VI-12,  and  VI-13  also  agree  reasonably 
well  with  an  exponential  curve. 

It  should  be  noted  that  Vanmarcke  (1974)  had  also  made  the 
assumption  that  x^is  exponentially  distributed.  A more 
complicated  exponential  type  curve  had  been  suggested  by 
Davidan  *t.  at. 
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Figure  VI -11 


Figure  VI-12 
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E.  DISCUSSIONS 


As  can  be  seen  from  the  comparisons  of  wave  data,  improvements 
on  both  of  the  above  approaches  are  necessary  so  that  more  accurate 
probabilities  can  be  predicted.  For  the  elapsed  time  t of 

the  envelope  above  a given  level,  a gama  curve  can  be  shown  to 
better  fit  the  data  than  an  exponential  curve. 

On  the  other  hand,  as  stated  earlier,  the  purpose  of  this 
section  is  to  obtain  the  probability  of  occurance  for  a deter- 
ministic wave  configuration.  In  practice,  a design  wave  is 
usually  described  by  a height  and  a period.  For  most  engineering 
purposes,  especially  for  inputs  to  a ship  capsize  simulator, 
tne  most  important  parameters  needed  to  describe  a wave  group 
seems  to  be  the  number  of  waves  in  a group  N,  the  level  H,  and 
the  mean  period  of  waves  in  the  group  T.  The  joint  probability 
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of  these  three  parameters  is  therefore  desired.  Due  to  the 
comple..  nature  of  multi-dimensional  joint  probabilities,  it 
would  be  desirable  to  separate  this  joint  probability  so  that 
each  parameter  can  be  studied  individually.  Using  Bayes'  rule, 
this  joint  probability  can  be  written  as 

P(H,N,T)  - P(H)P(N|H)P(T|H,N)  . (VI-20) 

Now  if  a time  series  is  interpreted  as  consecutive  outcomes  of 
individual  waves,  each  having  its  own  height  and  period,  this 
joint  probability  can  then  be  established. 

The  analysis  of  wave  record  "SEAWAY  14  AUG  73",  sequence 
1,  probe  I,  will  be  discussed.  In  this  example,  a wave  in  the 
time  series  is  defined  as  one  cycle  starting  and  ending  with 
a trough.  The  height  and  period  of  the  wave  is  defined  as  the 
distance  from  the  peak  to  the  average  of  two  troughs  and  time 
interval  between  the  two  troughs  respectively.  The  results  of 
the  analysis  are  shown  on  Figures  VI-14  through  IV-18. 

P (H) , as  shown  on  Figure  VI-14,  is  approximately  a Rayleigh 
distribution.  P(N  H) , at  four  different  levels  H,  is  shown  in 
Figure  VI-15.  As  stated  earlier  this  can  be  represented  by  Goda's 
method  reasonably  well  although  improver er ‘ is  necessary. 

P(T|H,N)  is  shown  on  Figures  VI-16  through  VI-18.  This  conditional 
probability  is  likely  to  cause  the  most  difficulty.  When  H 
and/or  N is  large,  there  will  not  be  enough  data  points  to 
estimate  this  probability  properly.  However,  from  these  figures, 
a normal  curve  seems  to  be  a reasonable  suggestion. 

A detailed  study  of  these  conditional  probabilities,  both 
theoretically  and  experimentally  is  in  order.  The  parameters  needed 
to  describe  each  conditional  probability  also  should  be  investigated 
so  that  the  joint  probability,  Equation  (VI-20) , can  be  estimated 
from  the  spectrum  rather  than  the  time  series. 
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Figure  VI-18 


VII  CONCLUSIONS 


Since  this  is  the  final  comprehensive  project  report,  it 
is  felt  worthwhile  to  include  n 1 ^ge  body  of  detailed  information 
which  miqht  otherwise  not  have  been  reported.  In  several 
places  in  the  report,  it  is  noted  that  the  present  results  are 
tentative  and  that  further  comparisons  and  analysis  are  definitely 
indicated, especially  when  new  full-scale  data,  e.g.,  storm 
wiv<  statistics,  become  available. 

The.'  study  was  initiated  with  the  optimistic  hope  of  gaining 
both  insight  into  the  fundamental  nature  of  ship  capsizing  and 
gathering  statistics  upon  which  to  base  stability  criteria.  As 
a secondary  objective  it  was  hoped  to  be  able  to  develop  analytical 
methods  by  means  of  which  capsizing  could  be  predicted  and  para- 
metric studies  could  be  performed.  Some  of  these  objectives  have 
been  fulfilled  only  to  a limited  extent.  While  the  fundamental 
nature  of  > .psizi.ng  of  intact  ships  in  heavy  seas  has  been  fairly 
well  revealed,  the  ability  to  obtain  rel i il  le  statistical  estimates 
of  the  probability  of  capsizing  of  a specific  ship  is  still 
lucking.  As  a specific  example,  it  is  clear  now  that  neither 
long  experimental  runs  nor  long  time-domain  computer  simulation 
provide  a feasible  means  of  amassing  the  statistical  base  needed 
for  such  predictions.  On  the  other  hand,  it  appears  possible  to 
simulate  on  a computer  the  short-term  large  amplitude  motion  of 
a ship  in  steep  quartering  or  following  seas.  At  the  same  time, 
other  analytic  techniques  have  been  developed  which  make  it  possible 
to  study  certain  nonlinear  motions  characteristic  of  capsizing 
situations.  These  techniques  offer  the  possibility  of  associating 
ship  capsizing  with  the  occurrence  of  a relatively  small  family 
of  wave  patterns.  The  probabilistic  nature  of  capsizing  is  then 
divorced  from  the  mechanism  of  ship  motion  and  is  derivable  solely 
from  the  statistics  of  the  waves. 

In  Summary 

1.  Linear  ship  mot  on  can  only  broadly  outline  areas  of  speed, 
heading,  and  ship  characteristics  which  may  load  to  trouble. 
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2.  Classical  ship  motion  theories  which  are  augmented  by 
some  nonlinear  terms  may  reveal  phenomena,  not  apparent  in  results 
of  linear  theory,  which  may  lead  to  severe  motion.  The  ability  to 
predict  capsizing  is  still  not  available  through  such  theories  but 
they  are  useful  in  parametric  studies  of  ship  characteristics  in 
relation  to  stability. 

3.  The  open  water  experiments  mentioned  above  yield  great 
insight  into  the  mechanics  of  capsizing.  The  phenomenon  is  a 
"rare  event"  in  the  statistical  sense  however,  and  very  extensive 
experimental  programs  would  be  required  to  adequately  describe 
the  probability  of  capsize  of  just  one  ship  configuration. 

4.  An  important  side  benefit  of  these  experiments  has  been 
the  development  of  methods  of  deducing  the  directional  properties 
of  wind  generated  waves  from  measurements  made  by  a sparse  gage 
array.  These  measurements  show  that  good  similarity  exists  between 
the  wave  conditions  in  the  model  test  area  and  full  scale  storm 
waves  at  sea. 

5.  The  knowledge  gained  in  the  experiments  has  led  to  the 
development  of  a time-domain  numerical  simulation  of  the  large 
amplitude  motion  of  a ship  in  following  and  quartering  seas.  Records 
of  severe  motions  including  capsizes  have  been  computed  by  this 
technique,  and  the  computed  motions  resemble,  quite  closely,  the 
motions  of  the  model  under  similar  sea  conditions.  This  technique 
appears  tobe  useful  in  defining  the  regions  of  wave  properties 
which  may  cause  capsizing. 

6.  The  characteristics  of  wave  groups  encountered  by  the 
model  have  considerable  influence  on  the  tendency  to  capsize. 

The  prediction  of  wave  group  occurrence  has  not  received  much 
attention  in  the  literature  yet  it  seems  to  play  a fundamental 
role  in  the  prediction  of  the  probability  of  ship  capsizing.  It 
appears  that  the  statistics  of  wave  groups  can  be  derived  from 

the  conventional  representation  of  the  random  seaway  and  therefore, 
the  study  of  wave  groups  is  a fruitful  subject  for  future  work. 
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Figure  A-  3 Righting  Arm  Curves  for  the  'h^t  1 mjcr  in  the  Heavy 

Condition 
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TABLE  A-2 

SL-7  CALIBRATION  FACTORS  FOR  1973 


ROLL 

+ 

83. 750  degrees 

/ 

+ 

1.35 

volts 

YAW 

+ 

66.466  degrees 

/ 

+ 

1.35 

volts 

PITCH 

+ 

22.620  degrees 

/ 

+ 

1.35 

volts 

RUDDER 

+ 

47 . 750  degrees 

/ 

+ 

1.35 

volts 

SWAY 

+ 

20.285  ft/sec' 

/ 

+ 

1.35 

volts 

SURGE 

+ 

7.  333  ft/sec2 

/ 

+ 

1.35 

volts 

HEAVE 

+ 

19.731  ft/sec2 

/ 

+ 

1.35 

volts 

RPM 

+ 

1088.71  RPM 

/ 

+ 

1.35 

volts 

+ 

8.6  ft/sec 

0 

1.35 

volts 

SPEED 

+ 

2.2  ft/sec 

0 

0.00 

volts 

TABLE  A- 3 


19  7 3 WAVE 

PROBE  ARRAY 

GEOMETRY 

PROBE  NO. 

DISTANCE 

BEARING 

(ft) 

(compass  deg.) 

1 

o 

o 

o 

o 

2 

9.208 

o 

o 

3 

9.292 

60.0 

4 

5.  125 

25.8 
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APPENDIX  B 
Roll  Damping 

The  coefficients  and  wave  exciting  forces  necessary  for 
linear  ship  motion  predictions  may  be  computed  using  potential 
flow  theory  for  an  ideal  inviscid  fluid.  Values,  thus  obtained, 
have  been  found  to  be  reasonably  accurate  when  compared  to  experi- 
mental results,  except  the  roll  damping  term,  for  which  the 
viscous  contribution  is  considerable.  Roll  damping  affects  the 
motion  response  only  in  a very  narrow  range  bracketing  roll 
resonance.  An  acceptable  value  for  the  roll  damping  coefficient, 
q , may  therefore  be  determined  from  experiments  in  which  the 
model  is  allowed  to  roll  at  its  natural  frequency.  The  theory 
and  procedure  for  evaluating  q from  roll  extinction  experiments, 
summarized  briefly  below,  are  given  in  (Vugts,1968  ). 

Vugts ' analysis  deals  with  the  motions  of  a cylinder  and 
will  be  assumed  applicable  to  a hull  form.  The  equation  of 
motion  in  roll  about  the  center  of  gravity  for  no  wave  excita- 
tion and  no  induced  sway  is 

£ , ,0  + q 4 ,0  + i*  0 - 0 

* r y Vt 

(B-l) 

where  Ais  the  hydrostatic  restoring  moment  ana  F is  a term 

¥ t ♦'  V 

including  the  effects  of  the  transverse  moment  of  inertia  and  the 
roll  added  mass  about  the  center  of  gravity.  /•  is  easily 

t r 

calculated  by: 


H * 


V GM 


( B — 2 ) 


and  P ^ may  be  approximated  from  the  natural  period  of  rolling 
from  the  following  equation, 

4irV, 


* / 7p  r 71 

\J  ri.  rt.  <i.  4.  ^ , 


44  14 


(B-3) 


with  the  assumption  that  the  motion  is  lightly  damped  or  that 
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to  rtn 


■j  - . . )■  is  much  larger  than  q.  ..  This  gives 

: , : ; ^ 4’  '4 


(B-4) 


The  roll  damping  term  is  determined  from  the  measured 
decrease  in  the  oscillations  or  the  logarithmic  decrement, 


(B-5) 


where  : , , and  the  relationship  between  6 and  q A are 

y V 

defined  in  Figure  H-l,  taken  from  Vugts  (1968). 


Figure  B-l  Roll  Decrement 


Substituting  euuation  (D-3)  in  (B-5)  for  the  value  of  TV 

<4 

anc  solving  for  . . results  in 

t 

? 

(B-6) 

A significant  parameter  which  is  often  overlooked  in 
determining  the  roll  damping  is  forward  speed.  Blagoveshchensky 
proposes  that  as  a rolling  snip  moves  ahead  it  imparts  more 
energy  into  the  fluid  and  that  the  resulting  increase  in  damping 
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is  linearly  proportional  with  forward  speed.  The  significance 
of  the  increase  is  readily  noted  from  a comparison  of  the  roll 
extinction  curves  at  zero  speed  and  at  forward  speed  for  the 
SL-’’  model  in  the  same  condition  shown  in  Figures  B-2  and  B-3. 

An  increase  in  roll  damping  due  to  the  forward  speed  is 
taken  into  account  in  the  theory  presented  by  Salvesen,  Tuck, 
and  Faltinsen  (1970)  in  the  calculations  dealing  with  the  blunt 
stern  effect  on  the  ship.  To  use  this  feature  in  SKAWAY  for  a 
ship  like  the  ZL-7  which  has  a very  small  portion  of  its  transom 
below  the  load  waterline,  it  would  be  necessary  to  approximate 
the  point  along  the  hull  where  separation  occurs  and  to  remove  the 
stations  aft  of  this  point  from  the  calculations.  The  SL-7 
offsets  for  SEAWAY  calculations  were  not  modified  for  flow 
separation.  Rather  the  measured  increased  damping  due  to 

forward  speed  was  added  to  that  which  was  determined  from  the 
integral  of  the  two-dimensional  damping  coefficients  along 
the  hull. 

k Experiments  were  performed  to  determine  the  increase  in 

roll  damping  at  forward  speed  with  the  twofold  purpose  of  as- 
certaining a better  figure  for  linear  theory  computations  and 
of  determining  a reasonable  relationship  between  these  two 
variables . 

The  same  analysis  was  applied  to  the  forward  speed  roll 
decay  tests  as  for  the  still  water,  though  with  some  reservation. 
Imparting  an  initial  roll  angle  or  impulsive  roll  moment  to  the 
free  running  model  was  somewhat  difficult  to  do  without  also 
engendering  some  sway  or  yawing  motion.  This  task  was  actually 
accomplished  by  following  the  model  in  the  rubber  boat  and 
pushing  one  deck  edge  of  the  model  down  at  about  amidship  by 
means  of  a fairly  long  pole.  Trom  an  observation  of  this  pro- 
cedure, it  was  noted  that  any  resulting  yaw  or  sway  motions 
were  damped  out  very  quickly  so  that  any  such  motion  could  be 
considered  negligible  after  only  one  half  of  one  roll  cycle. 

Before  commenting  upon  the  results,  one  point  should  be 
clarified.  It  was  desired  to  determine  a value  for  j . ( which 


L 
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would  give  a good  linear  approximation  of  the  model's  motions 
in  the  extreme  seaway  (scaled)  in  San  Francisco  Bay.  However, 
for  large  roll  amplitudes  the  restoring  moment  is  not  linear 
and  indeed  the  period  between  successive  cycles  is  changing. 
Since  the  results  were  to  be  used  as  input  for  the  linear  ship 
motions  program,  SEAWAY,  in  order  to  be  consistent,  all  cal- 
culations were  made  using  a value  of  T.  determined  from  zero 
speed  roll  extinction  curves  after  the  motion  had  diminished 
significantly  and  the  period  was  fairly  constant.  Generally, 
the  roll  period  for  the  initial  roll  cycles  (roughly  a 15 
degree  amplitude)  was  about  10%  less  than  the  natural  period  of 
roll . 

Values  of  the  roll  double  amplitude  were  plotted  against 
the  number  of  cycles.  A constant  $ . , approximately  equal  to 


30  degrees,  was  chosen  and  appropriate  values  of 


picked 


off  the  curves  for  an  increment  of  one  cycle.  An  example  is 
given  in  Figure  B-4 . 

Roll  decay  experiments  were  run  at  zero  forward  speed  and 
at  all  three  speed  settings  for  the  TL- ' in  the  heavy  condition 
and  in  the  light  condition  with  and  without  bilge  keels.  The 
results  of  these  tests  are  listed  in  Table  B-I  and  presented 
schematically  in  Figures  B-5  and  b-6. 

A nor. dimensional  coefficient,  r is  used  to  compare 
roll  damping  ac  different  speeds  for  the  different  model  test 
conditions.  It  was  computed  from  the  following  relationship 

<7 . , 

a - y i 


V"*-  • 

a v 


( B — 7 ) 


q is  the  total  roll  damping  including  both  viscous  and  hydro- 

V > 

dynamic  contributions,  r is  the  roll  restoring  moment,  and 


xx 


is  the  moment  of  inertia  of  the  model  in  air  about  a 


ait- 


longitudinal  axis  through  the  center  of  gravity.  For  these 


computations 


xx 


is  the  difference  between  the  total  moment  of 


ur 


J 

tv 


<\J  K)  ■'i'  L/l 


o n f o 


o 3j  i inn  - s^tonv  7 7g^  (xvw)  v/jv 
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he  actual  speed  will  differ  depending  upon  the  condition  of  the  model,  light  or  hea 
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Figure  IV-5  Roll  Damping  vs  GM 
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inertia  computed  with  the  GV.  and  the  roll  period  in  salt  water 
and  the  hydrodynamic  added  mass.  Mote  that  the  values  for  the 
hydrodynamic  damping  given  in  Table  B-I,  are  small  compared  to 
the  total,  q These  were  computed  by  the  linear  ship  motions 

t y • 

program  "SEAWAY"  with  an  offset  deck  which  did  not  include  the 
bilge  keels  and  therefore  are  somewhat  misrepresentat ive  of  the 
actual  hydrodynamic  damping.  Foi  the  light  condition  without 
bilge  keels,  though  iJ  is  a useful  measure  of  the  total  viscous 
damping . 

The  increase  in  roll  damping  with  forward  speed  is  believed 
to  be  attributable  to  a transverse  moment  developed  from  sway- 
lift  forces  on  the  submerged  portion  of  the  hull.  The  sway 
involved  here  is  not  that  of  the  center  of  gravity.  Rather  it  is 
the  horizontal  motion  of  a point  on  the  hull  due  to  rolling  and 
therefore  porportional  to  the  vertical  distance  from  the  center 
of  roll.  If  the  lift  force  of  the  bare  hull  is  assumed  to  act 
at  a vertical  height  equal  to  the  center  of  buoyancy,  and  the 
effect  of  the  rudder  is  included,  the  equation  for  the  damping 
moment  would  appear  as  follows: 


1 

1 

“ — f ’ V 

r 

.4.  d.  + 9 A 

!cl) 

d 

2 

. 

n \a  f.  tt  r r 1 

1 

' r 

r 

where 

= hull  angle  of  attack 

A,  = one-half  the  submerged  surface  area  of  the  hull  and 

n 

will  be  assumed  to  be  equal  to  the  length,  i times 
the  draft, 

= the  coefficient  of  lift  of  the  hull  per  unit  angle  of 
attack 

dy  = the  vertical  distance  from  the  center  of  buoyancy  to 
the  center  of  roll 

= rudder  anqle  of  attack 

i 

A = the  surface  area  of  one  side  of  the  rudder 
r* 
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C 


L 


u 


r 


the  coefficient  of  lift  of  the  rudder  per  unit  angle 
of  attack 


J = the  vertical  distance  from  the  center  of  pressure  of 
r 

the  rudder  to  the  center  of  roll. 


v = the  forward  velocity 
p = the  density  of  salt  water 
Since  the  angles  of  attack. 


6.  = — 

h v 


= - py-i 


; 0 = 
r 

(^j  d2 

n \ a lh  , 


[C . , 

A — ' i2 
rl  ilr  r 


(B-9) 


If  the  above  expression  is  equated  to  a first  order  damping 
term  which  is  linear  in  forward  velocity, 


and 


where 


j 


2 


q „ . = ( j + D V ) 0 
11  1 2 

is  the  zero  speed  linear  damping  coefficient,  then, 


(B-10) 


Comparisons  were  made  between  the  measured  increase  in 
damping  with  forward  speed  and  that  calculated  using  equation 
B-9.  In  general,  the  theoretical  increase  in  damping  approximated 
quite  closely  the  ave^  ye  increase  which  was  measured. 

Generally,  the  experimental  results  show  a great  deal  of 
scatter  for  relatively  few  points.  Fortunately,  some  useful 
qualitative  conclusions  may  be  drawn.  The  curves  of  q 4 for 

v V 

zero  forward  speed  showed  little  correlation  except  in  general 
magnitude  between  the  light  and  heavy  conditions.  3 values  were 
fairly  close  and  would  appear  to  be  a more  reasonable  measure 
for  comparison  between  different  ships  or  conditions.  Removal 
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of  the  bilge  keels  resulted  in  a general  decrease  in  roll  damping 
of  about  25  to  30  percent.  Damping  was  significantly  greater  at 
forward  speed  by  a factor  of  3 or  4.  The  large  variations  between 
points  illustrated  in  Figure  B-5  indicate  the  desirability  for 
plotting  the  points  with  a series  of  straight  lines  rather  than 
a faired  curve.  Due  to  the  relatively  good  agreement  between  the 
calculated  and  measured  increases  in  roll  damping  associated  with 
forward  velocity,  it  is  reasonable  to  assume  a linear  relationship 
between  roll  damping  and  forward  sueed.  These  results  are  far 
from  conclusive  and  clearly  demon?  trate  the  need  for  a series  of 
experiments  which  better  covers  the  range  of  important  variables 
and  from  which  a valuable  confidence  interval  may  be  determined. 
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APPENDIX  C 

COMPARISON  AND  DETERMINATION  OF  AUTO  PILOT  PARAMETERS 

The  autopilot  control  for  the  SL-7  model  is  a Proportional 
Differential  (PD)  control  circuit  with  a dead  band  or  "weather 
adjustment".  A block  diagram  of  the  autopilot  for  the  model  is 
shown  in  Figure  C-l. 


4 IKiddar  Maponaa  tafia 


Figure  C-l  Block  Diagram  Autopilot  Model 

The  system  is  setup  such  that  the  rudder  angle  is  proportional 
to  the  yaw  and  the  yaw  rate  as  in  the  following  equations: 

6 = 6 +6  = a 8 + a 3 

12  1 2 — o 

6=o  le|<e0 


where 


6 = rudder  angle 

6 1 = yaw  proportional  rudder  response 

6 2 = yaw  rate  proportional  rudder  response 

3 = yaw  angle 

3q  = dead  band  (weather  adjustment) 

a , a are  constant  coefficients 
1 2 


I 
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The  yaw  angle  is  positive  when  the  bow  of  the  model  turns  to 
port;  and  the  rudder  angle  is  positive  when  the  rudder  turns  to 
starboard. 

The  purpose  of  the  dead  bard  or  "weather  adjustment"  is  to 
stop  excessive  rudder  motions,  particularly  in  rough  seas.  It 
is  generally  recommended  for  ship  autopilots  that 


6o=° 


in  good  weather 

in  rough  seas,  in  order  to  prevent 
unnecessary  or  excessive  rudder  motion 


The  value  of  the  coefficients  a and  a depend  on  the  type 

1 2 

of  ship  and  the  operating  conditions.  They  are  generally  de- 
termined by  experiment.  Typical  values  for  B , ai  and  a2  are: 

B 3/4°  to  4° 

c 

a i 1 to  4 deg  / deg 

a2  4 to  20  deg  • sec  / deg 

a 2 

— 4 to  5 see 


During  the  1973-74  test  season , yaw  and  rudder  motions  were  recorded 
as  functions  of  time.  Based  upon  this  data  we  wish  to  do  two  things: 


Compare  the  static  calibration  constants  with 

the  coefficients  (a  , a ) obtained  from  a linear  equation, 

1 2 

describing  the  experimental  rudder  motion  as  a function 
of  the  yaw  and  yaw  rate. 

Determine  a more  realistic  equation  representing  the 
experimental  rudder  motion  as  a function  of  the  yaw  and 
it's  derivatives. 


The  calibration  of  the  autopilot  was  done  in  air  and  the  complete 
description  of  the  procedure  is  reported  in  Haddara  (1971)  . A 
similar  procedure  was  used  for  the  1972  and  1973  experiments. 
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The  linear  equa.ion  used  to  determine  the  response  of  the  rudat  i 
as  a function  of  yaw  and  yaw  rate  was: 

6 = a t 2 + a t 

The  yaw  rate,  3 , was  determined  by  fitting  a parabolic 

curve  through  every  three  points  of  the  yaw  vs  time  recoid  end 

calculating  the  slope  at  each  of  these  points.  The  method  of 

Least  Squares  was  used  to  calculate  a,  and  a . Up  to  20 u [ wi.t 

(6,2,3,  n = 1,200)  were  used.  Since  ntosl  of  the  model  run 

n n n 

had  more  than  200  samples,  a number  (200)  of  evenly  spaced  pc  . i.t 
(same  time  interval)  were  used.  The  results  and  the  calibi  ate  : 
values  are  shown  in  Table  C-I 


TABLE  0-1 


AUTO 

PILOT  COEFFICIENTS  FOR 

A LINEAR 

RUDDER 

L'Ql  A’l  i ON 

DATE 

RECORD 

# SPEED 

V 

HEAD- 

ING 

a 

l 

SEPT  17 

1 

3 

0.54 

F 

1 . JO 

. O l 

2 

3 

0.59 

F 

1.04 

. b i 

3 

2 

0.69 

Q 

1.34 

. 5 5 

4 

2 

0.75 

F 

1.1) 

• / 0 

5 

6 

2 

0.82 

Q 

1.12 

. j] 

7 

3 

0.93 

F 

.95 

•*  . 1 .) 

8 

2 

0.97 

F 

. 9 8 

.29 

9 

1 

1.02 

r 

1 . i 6 

. 4 i 

CALIBRATION 

1 . it) 

The  17th  of 

September  was  chosen  for  analysis  because  it  is 

a 

very  representative 

day  and  included  many  r 

uns  at 

di  f i 'lent 

S ( Ce 

The  a 

j and  a? 

obtained  from 

the  calibration 

were  i 

"*d  .ii 

larger  than 

the  a j 

and  a2  obtained 

from  the 

lineal 

eqit  c tic 

Possible  reasons  for  this  would  be: 


1.  The  rudder  motion  lags  the  yaw  by  some  time  i h : 'Ml 
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2.  The  dead  band  was  not  accounted  for  in  the  equation. 

This  could  significantly  alter  the  constants  a(  and  a2 

3.  If  the  rudder  is  hard  over,  the  rudder  angle  is  constant 
and  no  longer  dependent  on  6 or  8.  If  such  were  the 

case,  then  again,  the  accuracy  of  aj  and  a?  would  be  limited. 

Examination  of  experimental  values  of  the  rudder,  yaw  and 
yaw  rate  show  that  the  value  obtained  for  the  at  and  a2  by  the 
Least  Squares  method  are  consistent  with  the  data.  It  appears 
that  of  the  three  reasons  given  for  possible  error  in  the  co- 
efficients, the  third  may  be  discarded,  as  this  does  not  happen 
in  the  data  to  any  great  extent. 

Another  equation  was  tried  incorporating  a time  lag,  t. 

Assuming  that  the  response  of  the  rudder  lags  the  yaw,  we  assumed 
that 

6 = a1  fci  (t-t)  + a26(t-'t) 

t = t i me 
t = time  lag 

Expanding  the  above  equation  in  a Taylor’s  series  we  obtain 

2 ...  3 2 

- a j ( 3 - t 8 + Bjy  - * ■ • ■ ) a2  ( r - x 8 + ^2~T  • • • • ^ 

which  can  be  approximated  as: 

5 = a,  8 + (a  -xa^Bj  + (yj-  al  - xa2)8 

* (rr  a2  " t W 

and  8 were  calculated  by  taking  the  derivative  of  the  yaw 
rate  (8)  . The  method  of  Least  Squares  was  used  to  calculate 
the  coefficients  a.,  a , and  x.  Again  a maximum  of  200  sets 

of  [joints  was  allowed. 
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The  results  of  this  equation  are  shown  in  Table  C-1I. 


TABLE  C-II 


AUTOPILOT  COEFFICIENTS 

FOR  A 

PARABOLIC 

: RUDDER 

EQUATION 

DATE  RECORD  # 

SPEED 

H 

HEAD- 

a 

a 

X 

7; 

ING 

i 

2 

SEPT  17  1 

3 

0.54 

F 

1 .10 

.99 

.32 

2 

3 

0.59 

F 

1.02 

.90 

.25 

3 

2 

0.69 

Q 

1.32  1 

. 13 

.32 

4 

2 

0.75 

F 

1.10  1 

.04 

.28 

5 

6 

2 

0.82 

Q 

1.12 

.40 

.07 

7 

3 

0.93 

F 

.94 

. 75 

.70 

8 

2 

0.97 

F 

.97 

.59 

.29 

9 

1 

1.02 

F 

1.14 

. 80 

.29 

In  going  from  the  linear  equation  to  the  parabolic  equation 

the  a 's  remained  essentially  the  same,  while  the  a 's  changed  in 
I 2 

the  neighborhood  of  50%.  It  appears  that  the  calculated  time  lag 
is  smaller  than  that  estimated  from  the  model  runs  by  up  to  factor 
of  3.  The  parabolic  fit  seems  to  indicate  that  t is  almost 

insignificant  for  this  day.  The  increased  a ' s are  probably  due, 

2 

however,  to  the  fact  that  x is  included. 

It  must  be  remembered  that  the  dead  band  was  not  included  in 

either  of  the  fitted  equations.  But  as  pointed  out,  the  computed 

values  of  a , a , and  x are  reasonably  consistent  with  the  data. 

1 2 
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APPENDIX  D 


SUMMARY  OF  MODEL  RUNS  AND  SEASTATE 
INFORMATION  FOR  1972 
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FIGURE  D-l 


SUMMARY  OF  MODEL  RUNS  AND  SEASTATE  INFORMATION 
VS.  TIME  OF  RECORDING  - 1972 


GM  No. 

SPEED  No. 

Seaway 
Re co rds 


See  Table  A-I  for  the  American  Challenger. 

See  Table  A-I  for  the  American  Challenger. 

Sequence  of  recordings  including  bad  starts, 
etc. 

Assumed  Seaway  Energy  or  Wave  Amplitude 
90%  Confidence  Bound 
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TABLE  D-I 

STATISTICAL  VALUES  FOR  ALL  WAVE 
RECORDS  FROM  FFT  SPECTRA  1972 


RMS  Ampli tude  Computed  from  the  Spectra  using  the 
Fast  Fourier  Transform  Method  of  Analysis  for  the 
Four  Wave  Probes. 
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TABLE  D-I 

Statistical  Values  for  All  Wave 
Records  from  FFT  Spectra 


RHS  Amplitude  Ft.  (Deviation  from  Mean) 


DATE 

TIME 

I 

— 

11 

— 

Ill 

IV 

r 

Aug. 

10,  1972 

11.20 

0.2015 

0.2232 

0.2190 

11  .34 

0.2207 

0.2523 

0.2729 

1 .14 

0.3037 

0.3301 

0.3154 

0.3382 

2.57 

0.3948 

0.4099 

0.4040 

0.4256 

4.35 

0.3497 

0.3635 

0 . 381 7 

0.3813 

5.40 

0.3159 

0.3395 

0.3461 

0.3532 

Sept . 

5 . 1972 

11.55 

0.1818 

0.1826 

0.1768 

0.1752 

12.25 

0.2248 

0.2231 

0.2048 

0.2046 

1.59 

0.2908 

0.3126 

0.2793 

0.2698 

3.08 

0.3115 

0.2949 

0.2924 

0.2871 

4.43 

0.3068 

0.2997 

0.2767 

0.2860 

Sept. 

7,  1972 

1.00 

0.2634 

0.2316 

0.2327 

0.2408 

2.37 

0.3141 

0.2978 

0.3062 

0.2949 

3.  55 

0.3526 

0.3201 

0.3206 

0.3289 

4.58 

0.3280 

0.3071 

0.2987 

0.3133 

Sept . 

11,  1972 

11.03 

0.2345 

0.1925 

0. 1729 

0.2080 

12.32 

0.2707 

0.2753 

0.2551 

0.2606 

1 . 52 

0.2840 

0.2644 

0.2261 

0.2629 

3.04 

0.2764 

0.2570 

0.2362 

0.2517 

4.34 

0.3018 

0.2740 

0.2929 

0.2846 

Sept . 

13,  1972 

11  . 34 

0.2622 

0.2370 

0.2334 

0.2399 

1.22 

0.2894 

0.2687 

0.2591 

0.2732 

Sept. 

18,  1972 

12.40 

0.1785 

0.1838 

0.1948 

0.1875 

2.21 

0.1940 

0.1978 

0.2186 

0.2076 

Sept . 

21,  1972 

12.41 

0.1532 

0.1566 

0.1698 

0.1605 

1.27 

0.2015 

0.2095 

0.2168 

0.2188 

2.52 

0.2424 

0.2803 

0.2707 

0.2726 

4. 31 

0. 1712 

0.1846 

0.1970 

0. 1830 

Sept . 

25,  1972 

12.36 

0.2221 

0.2246 

0.2400 

0.2429 

1 . 54 

0.2309 

0.2297 

0.2385 

0.2356 

3.  12 

0.2038 

0.2041 

0.2178 

0.2074 

4.51 

0.1929 

0.1507 

0.  1928 

0.1804 
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TABLF  D-I 

Statistical  Values  for  All  Wave 
Records  front  FFT  Spectra 


RMS  Amplitude  Ft.  (Deviation  from  Mean) 


DATE 

TIME 

I 

1 1 

III 

IV 

Sept.  27,  1972 

i 2 . 5 2 

0.1427 

0.1315 

0.1602 

0.1330 

2.29 

0.1880 

0.1650 

0.1949 

0.1855 

3.  31 

0.1970 

0.1833 

0.2046 

0.1904 

4.38 

0.1822 

0.1947 

0.1970 

0.1887 

Oct.  3,  1972 

10.58 

0.1423 

0.1345 

0.1259 

0.1328 

12.29 

0.2280 

0.1983 

0.2292 

0.2140 

2.06 

0.2237 

0.2063 

0.2255 

0.2162 

3.23 

0.2474 

0.2213 

0.2430 

0.2435 

5.10 

0.2797 

0.2717 

0.2665 

0.2616 

1 
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TABLE  D-II 


STATISTICAL  MEANS  OF  SEAWAY  SIGNIFICANT 
AMPLITUDES , RECORDED  IN  SAN  FRANCISCO 
BAY  1972 


i\\.  Significant  Amplitude,  i.e. 

'»  2*RMS  derived  from  the  time 

series. 

% 

Average  of  the  four  probes. 

S^  Standard  Deviation  of  the 

Mean. 
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TABLE  D-II 


Statistical  Means  of  Seaway  Significant  Amplitudes 


Recorded  in  San  Francisco  Bay  1972 


PROBE  NO. 


SAMPLE 

No. 


8/10/72  1.01 


9/7/72 


1.04 
1 .05 
1.06 


2.01 
2.02 
2 . 0 j 
2.06 
2.05 


3. 

3. 

3.04 

3.05 


9/11/72  4.01 


4.04 

4.05 


^1/3  ’ Ft*  ■ 


♦ampl i tudes 


3746  .3618  .3262  .3722  .3587 


.0312  1. 


.2711  .2400 
.3057  .3655 
.3161  .3490  .3274 
.4474  .4485  .4403 
.3694  .4014  .4165 
.3585  .3816  .3840 


1898 

2304 

2859 

3186 

3142 


2711 

3284 

3878 

3798 


1876 

2253 

2817 

3095 

3367 


2517 

3050 

3416 

3416 


2527  .1937 
3006  . 3031 
3177  .2823 
2954  .2620 
2995  .2691 


2698 

2997 


1838  .1826 
2041  .2117 


.1868 
.2758 
.2389 
.2595 
. 3029 


.1948 
. 2242 


3480 

4748 

4219 

3969 


1875  .1789 
2097  .2031 
2757  .2641 
3126  .3066 
2925  .3133 


2499  .2462 
3060  .2969 
3296  .3487 
3125  . 3379 


6 
4 

3010 


.2567 

.3297 

.3351 

.4527 

.4023 

.3802 


.1859 

.2171 

.2/68 

.3118 

.3142 


.2547 

.3091 

.3519 

.3429 


.2159 

.2917 

.2850 

.2699 

.2889 


74 
79 
81 
76 
0113 
0080 


0126 
01  39 


56 
63 
0170 
0085 
0078 


24 

64 

47 

.0171 

26 

.0082 

90 

.0287 
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TABLE  D-II 


Statistical  Means  of  Seaway  Significant  Amplitudes, 


Recorded  in  San  Francisco  Bay  1972  (cont.) 


^1/3*  Ft* 


PROBE 

NO.-*. 

DATE 

SAMPLE 
No . 

9/21/72 

7.01 

7.02 

7.03 

7.04 

9/25/72 

8.01 

3.02 

8.03 

8.04 

9/27/72 

9.01 

9.02 

9.03 

9.04 

10/3/72 

10.01 

10.02 

10.03 

10.04 

10.05 

2 


ampl i tudes 


.1606 

.1806 

.2403 

.1848 


.2424 

.2413 

.1902 

.1569 


.1902 

.2120 

.2155 

.1687 


.1628 

.2334 

.2187 

.2576 

.3066 


.1642 

.1973 

.2791 

.2004 


.2368 

.2313 

.1916 

.1453 


i 


703 

881 

882 

1862 


.1554 

.2016 

.1896 

.2344 

.2883 


1794 

2007 

2685 

2020 


,2602 

2474 

2034 

1438 


1982 

2215 

2154 

1829 


,1488 

2237 

,2183 

2487 

2956 


1678 

2144 

2681 

1925 


1395 


90 
30 
36 
24 
2764 


fil/3 

sn 

sf\ 

Ft. 

Ft. 

^1/3 

,1680 

,1982 

,2640 

,1949 


0040 


2546  .0096 
2391  .0034 
1947  .0030 
1464  .0037 


.1850 

.2083 

.2057 

.1773 


.1522 
. 2154 
.2100 
.2458 
.2917 


45 
78 
69 
49 
.0063 


0243  1. 

0350  2. 

0314  3. 

0203  4. 


0324 

0342 

0314 

0240 


0293 

0364 

0329 

0200 

0217 


* s = Standard  Deviation  of  the  Mean 
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TABLE  D-III 

STATISTICAL  MEANS  OF  SEAWAY  ENERGIES, 
RECORDED  IN  SAN  FRANCISCO  BAY  1972 


k 


nRMS 


Root  Mean  Square  Wave 
Amplitude 

] 

Average  Mean  Square  Value 


Standard  Deviation  of  the 
Mean 


Mean  Energy  Ship  Scale  assuming  X =30. 


j 
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TABLE  D-III 


Statistical  Means  of  Seaway  Energies 


Recorded  in  San  Francisco  Bay  1972 


PROBE 

NO. 

DATE 

SAMPLE 

No. 

8/8/72 

0.01 

8/10/72 

1.01 
1 .02 
1 .03 
1 .04 
1 .05 
1 .06 

9/5/72 

flHM 

IJfll 

2.04 

2.05 

9/7/72 

HH 

i: 

3.04 

3.05 

9/11/72 

4.01 

4.02 

4.03 

4.04 

4.05 

9/13/72 


9/18/72  6. 

6. 


9/21/72  7.01 

7.02 

7.03 

7.04 


.268  .308  .0821  .0049 


,201  .218 
282  .268 
315  .337 
399  .423 
367  .368 
344  .351 


173  .173 
199  .200 
262  .257 
284  .277 


.238 
.280 
307  .31  5 
291  .307 


0441  .0025 
0649  .0070 
1030  .0049 
1646  .0054 
1271  .0052 
1136  .0051 


0309 
0443 
0721 
0825  .0031 
0848  .0040 


0490  .0054 
0826  .0026 
0884  .0073 
0934  .0039 


0376 

0692 

0583 

0647 

0737 


49 
23 
64 
. 1043 
.0032 


186 
566 
379 
0371 
0*67 


099 
319 
825 
0421 


300 
328 
098 
0570 
0429 


mean 

energy 

ship 

scale 

Ft.2 


Seq . 
No. 


44.66 

75.28 

80.57 

85.12 


67 

17 

49 

94 

65 

89 

29 

71 

38 

10 

20 

60 

32. 

90 

57. 

42 

29. 

89 
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TABLE  D-III 


Statistical  Means  of  Seaway  Energies 


Recorded  in  San  Francisco  Bay  1972  (cont.) 


PROBE 

: no.  - 

DATE 

SAMPLE 
No . 

9/25/72 

8.01 

8.02 

8.03 

8.04 

9/27/72 

m 

9.03 

9.04 

10/3  7 7 2 

10.01 

10.02 

10.03 

10.04 

10.05 

nRMS ’ Ft‘ 


ordinates  -*• 


220  .238  .241 
226  .234  .231 
189  .200  .191 
131  .138  .128 


30  .150  .133 
65  .194  .185 
82  .200  .183 
1 75  . 1 76  .1  67 


.124  .121 
.216  .201 
.205  .196 
.240  .240 
.259  .259 


mean 
ene  rgy 


52  7 
526 
368 
.0191 


.0193 
336 
366 
.0288 


.0160 
423 
391 
554 
.0709 


26  .0500 
08  .0155 
0011  .0295 
0017  .0878 


13  .0666 
23  .0672 
0015  .0420 
0013  .0453 


mean 

energy 

ship 

scale 


48.03 

47.94 

33.54 

17.41 


1 7.59 
30.62 
33.36 
26.25 


14.58 

38.55 

35.63 

50.49 

64.62 
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TABLE  D-IV 


STATISTICAL  MEANS  OF  SEAWAY  ZERO  CROSSING 
PERIODS,  RECORDED  IN  SAN  FRANCISCO  BAY 
1972 


T lezo  crossing  periods  obtained 

from  the  time  series 

T Mean  period 

SE  Standard  Deviation  of  the  Mean 
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TABLE  D-1V 


Statistical  Means  of  Seaway  Zero  Crossing  Periods, 
Recorded  in  San  Francisco  Bay  1972 


PROBE 

NO.  - 

T . 

Sec . 

T 

Sec . 

SE 

Sec. 

Seq . 
No. 

1 

2 

3 

4 

DATE 

SAMPLE 

No. 

«-  periods  -*■ 

mean 
peri od 

* 

8/8/72 

0.01 

2.0000 

1 .7121 

1 .6377 

1.8226 

1 . 7931 

0.0787 

0.0439 

1 . 

8/10/72 

1.01 

1.6359 

1.4902 

1.5701 

1.5464 

T .5606 

0.0302 

0.0193 

1 . 

1.02 

1 . 7393 

1 .6351 

1.6031 

1.6483 

1 .6564 

0.0292 

0.0176 

2. 

1.03 

2.0078 

1.9007 

1.8767 

1.9669 

1.9380 

0.0301 

0.0155 

11  . 

1.04 

2.3779 

2.3141 

2.3473 

2.3473 

2.3466 

0.0130 

0.0056 

12. 

1.05 

2.5052 

2.5721 

2.4601 

2.4787 

2.5040 

0.0245 

0.0098 

13. 

1.06 

2.4545 

2.3956 

2.4453 

2.3956 

2.4227 

0.0158 

0.0065 

14. 

9/5/72 

2.01 

1 .6368 

1.6483 

1 .6016 

1 .6516 

1 .6346 

0.0114 

0.0070 

1 . 

2.02 

1 . 7301 

1 .8004 

1.7965 

1.7926 

1 . 7799 

0.0167 

0.0094 

2. 

2.03 

2.1558 

2.2661 

2.2630 

2.31 74 

2.2506 

0.0340 

0.0151 

3. 

2.04 

2.2536 

2.2915 

2.3207 

2.3011 

2.2917 

0.0141 

0.0061 

4. 

2.05 

2.1473 

2.1874 

2.1962 

2.1  701 

2.1752 

0.0108 

0.0050 

5. 

9/7/72 

3.02 

1 .8471 

1.9029 

1 .8703 

1.8703 

1 .8726 

0.0115 

0.0061 

1 . 

3.03 

2.2883 

2.2322 

2.2201 

2.2567 

2.2493 

0.0151 

0.0067 

2 

3.04 

2.4165 

2.3849 

2.3574 

2.4094 

2.3920 

0.0134 

0.0056 

3. 

3.05 

2.2475 

2.3044 

2.3076 

2.2413 

2.2752 

0.0178 

0.0078 

4. 

9/11/72 

4.01 

1 3650 

1 . 7847 

1 . 7264 

1 .7246 

1 .7252 

0.0244 

0.0142 

1 . 

4.02 

1 .8471 

1 . 8725 

1 .8597 

1 .851 3 

1 .8576 

0.0056 

0.0030 

2. 

4.03 

1 . 7393 

1 .8084 

1.8576 

1 .8144 

1.8049 

0.0245 

0.0136 

3. 

4.04 

2.0054 

1 .9669 

2.0871 

1 .9883 

2.0119 

0.0263 

0.0131 

4. 

4.05 

2.1 586 

2.1  758 

2.1904 

2.1816 

2.1  766 

0.0067 

0.0031 

5. 

9/13/72 

5.01 

1 . 7542 

1 . 71 74 

1 . 7693 

1.8225 

1 . 7658 

0.0218 

0.0123 

1 . 

5.02 

1 . 7809 

1.8388 

1.9366 

1.9692 

1.8814 

0.0434 

0.0231 

2. 

9/18/72 

6.01 

1 .8388 

1 .8534 

1.8327 

1.8044 

1 .8323 

0.0103 

0.0056 

1 . 

6.02 

1.8347 

1 .8388 

1 . 8985 

1 .9029 

1.8687 

0.0185 

0.0099 

2. 

9/21/72 

7.01 

1.6926 

1 .6566 

1 .6583 

1.6368 

1.6611 

0.0116 

0.0070 

1 . 

7.02 

1 .9716 

1 .9140 

1 .9669 

1.9575 

1 .9525 

0.01  32 

0.0067 

2. 

7.03 

2.1501 

2.2171 

2.1  558 

2.1195 

2.1606 

0.0204 

0.0095 

3. 

7.04 

1.9669 

1 .9908 

1 .9321 

1.9622 

1.9630 

0.0121 

0.0061 

4. 
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TABLE  D-IV 


Statistical  Means  of  Seaway  Zero  Crossing  Periods, 


Recorded  in  San  Francisco  Bay  1972  (cont.) 


PROBE 

NO.  * 

T , 

Sec. 

1 

*T 

Sec . 

SE 

Sec . 

SE 

BB 

I! 

1 

2 

3 

4 

T 

DATE 

E 

periods  -*■ 

mean 
peri  od 

* 

9/25/72 

8.01 

1.7067 

1.7598 

1.6770 

1 .635  1 

1.6946 

0.0262 

0.0155 

i . 

8.02 

1.8450 

1.8746 

1.7828 

1.8492 

1.8379 

0.0195 

0.0106 

8.03 

1.9859 

1.9811 

1.9253 

1.9575 

1.9624 

0.0139 

0.0071 

8.04 

1.9932 

1.9908 

1.9389 

1.9208 

1.9609 

0.0183 

0.0093 

9/27/72 

9.01 

1.5142 

1.5545 

1.5004 

1.5341 

1.5258 

0.0118 

0.0077 

i . 

9.02 

1 .6961 

1.7542 

1.6735 

1.6335 

1.6893 

0.0252 

0.0149 

9.03 

1 . 7828 

1.8164 

1.7467 

1.7712 

1.7793 

0.0145 

0.0081 

9.04 

1.9118 

1.9482 

1.9366 

1.9230 

1.9299 

0.0079 

0.0041 

10/3/72 

10.01 

1.8492 

1.6839 

1.6400 

1.6933 

0.0547 

0.0323 

i . 

2.0531 

2.0609 

2.0078 

2.0687 

2.0476 

0.0137 

0.0067 

2. 

■ 'HtI 

2.0153 

2.1141 

1.9787 

2.0353 

2.0358 

0.0286 

0.0140 

3. 

1.9789 

2.0177 

2.0128 

2.0029 

2.0030 

0.0087 

0.0043 

5. 

i 

K AMI 

1.8854 

1.9185 

2.0103 

1.9645 

1.9447 

0.0272 

0.0140 

6. 

s = Standard  Deviation  of  the  Mean 


A 
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TABLE  D-V 

OVERVIEW  OF  MODEL  TESTS  IN  SAN 
FRANCISCO  BAY  1972 


GM  No.  See  Table  A-I  for  the  American 

Challenger. 

Heading  Visual  description  of  the  model 

heading  vs.  the  apparent  wave 
direction. 

Speed  See  Table  A-I  for  the  American 

No . Challenger. 

Duration  Length  of  the  Model  Run. 

Digitized  Rate. 

Seq.  No.  Sequence  Number  assigned  during 

Calibration  of  the  Digitized  data. 


MOTION  DATA  RECORDED  FOR  THE  CHALLENGER  MODEL 
IN  THE  HEAVY  CONDITION  DURING  THE  1972  SEASON. 


YAW 


RUDDER  ROLL  PITCH  RPM 


13  SEP  - 


03  OCT 


* The  number  refers  to  the  channel  on  which  the  data  was 
recorded. 


Overview  for  Model  Tests  in  San  Francisco  Bay  1972 
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TABLE  D-VI 

STATISTICAL  RESULTS  FROM  MEASURED 
MODEL  DATA  1972 


The  highest  amplitude  statistics  were  com- 
puted from  an  analysis  of  the  time  series. 


The  seaway  statistics  were  interpolated  to 
the  time  of  the  model  run  and  are  intended 
only  as  a guide. 
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APPENDIX  E 


SUMMARY  OF  MODEL  RUNS  AND  SEASTATE 
INFORMATION  FOR  1973  AND  1974 
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Figures  E-l 


SUMMARY  OF  MODEL  RUNS  AND  SEASTATE 
INFORMATION  VS.  TIME  OF  RECORDS 
1973  and  1974 


Model  runs  are  indicated  by  symbols: 


GM  No. 
SPEED 


See  Table  A-I  for  the  SL-7  Model. 


See  Table  A-I  for  the  SL-7  Model. 


COURSE  F - 


Folowing  seas  (visual  observation) 
Quartering  seas 


Beam  seas 
Head  seas 


1.01  - Experiment  run  number 

Model  run  without  capsize: 

vertical  position  corresponds  to 
course  and  speed; 

horizontal  position  indicates  time 
of  day; 

horizontal  length  indicates  duration 
of  run. 

Model  run  with  capsize. 

Linear  interpolutions  of  seaway  measurements  are  indicated 
by  lines: 

Significant  wave  height. 


— Direction  of  spectrum  peak. 
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FIGURE  E-1.5 
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SUMMARY  OF  MODEL  RUNS  AND  SEASTATE  INFORMATION 
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TIME  OF  DAY  (HRS) 


SUMMARY  OF  MODEL  RUNS  AND  SEASTATE  INFORMATION 


FIGURE  E-I.8 
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TIME  OF  DAY  (HRS) 
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SUMMARY  OP  MODEL  RUNS  AND  SEASTATE  INFORMATION 
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TIME  OF  DAY  (HRS) 


SUMMARY  OP  MODEL  RUNS  AND  SEASTATE  INFORMATION 


FIGURE  E-l . 16 
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TIME  OF  DAY  (HRS) 


FIGURE  E-1.18 


SUMMARY  OF  MODFL  RUNS  AND  SEASTATE  INFORMATION 
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SUMMARY  OF  MODEL  RUNS  AND  SEASTATE  INFORMATION 


TIME  OF  DAY  (HRS) 


SUMMARY  OF  MODEL  RUNS  AND  SEASTATE  INFORMATION 


228 


ME  OF  DAY  (HRS) 


SUMMARY  OF  MODEL  RUNS  AND  SEASTATE  INFORMATION 


TIME  OF  DAY  (HRS) 
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FIGURES  E-2 

POIMT  SPECTRA  OF  PROBE-IV  FOR  1973  and  1974 


The  plots  contain  the  wave  spectra  of 
probe  IV  measured  on  a single  day  in  order  to  show 
the  spectral  growth  (or  decay)  throughout  the  test 
period. 
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SAN  FRANCISCO  BAY  WAVE  SPECTRA 
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Figure 
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Figure 
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Figure 


Figure  E-2.10 


SAN  FRANCISCO  BAY  WAVE  SPECTRA 
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Figure 


Figure 


SCO  BAY  UAVE  SPECTRA 


Figure  E-2 


Figure 


Figure  E-2.19 


Figure  E-2.20 
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Figure  E-2.21 
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Figure 
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FIGURES  E-3 
DIRECTIONAL  SPECTRA 

The  directional  spectra,  S(^f),  are 
presented  as  contour  polar  plots  with  fre- 
quency in  Hz  in  the  radial  direction.  Note 
that  the  computations  we  carried  out  at  constant 
frequency  and  may  have  resulted  ir  too  much 
spreading  in  0.  The  contouring  routine  was  not 
perfect  and  some  small  spurious  contours  are 
evident.  The  contour  levels  vary  but  not  enough 
to  include  levels  for  the  very  low  seastates. 

These  appear  as  blanks  in  their  respective  locations. 
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DIRECTIONAL  SPECTRA 
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FIGURE  E 3.10 
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DIRECTIONAL  SPECTRA 
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FIGURE  E 3.12 
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DIRECTIONAL  SPECTRA 
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FIGURE  E 3.16 
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DIRECTIONAL  SPECTRA 
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DIRECTIONAL  SPECTRA 
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DIRECTIONAL  SPECTRA 
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DIRECTIONAL  SPECTRA 
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DIRECTIONAL  SPECTRA 
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FIGURE  E 3.35 
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FIGURE  E 3.38 
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DIRECTIONAL  SPECTRA 
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DIRECTIONAL  SPECTRA 
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DIRECTIONAL  SPECTRA 
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DIRECTIONAL  SPECTRA 
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FIGURE  E 3.60 
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DIRECTIONAL  SPECTRA 
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FIGURE  E 3.61 
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FIGURE  E 3.62 


DIRECTIONAL  SPECTRA 
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DIRECTIONAL  SPECTRA 
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FIGURE  E 3.65 
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FIGURE  E 3.66 
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TABLE  E-I 

OVERVIEW  OF  WAVE  RECORDS  1973  and  1974 


Abbreviations 

BD 

8 

Depth  of  water  at  the  buoy  (ft) . 

BDR 

■ 

Buoy  direction  from  probe  1 to  probe  2 
unless  otherwise  noted  (deg) . 

CDR 

8 

Current  direction  (deg) . 

CS 

- 

Current  Speed  (knots) . 

E 

= 

Ebb  current. 

F 

S 

Flood  current. 

HVs 

f 

Significant  wave  height. 

s 

Slack  water. 

WDR 

= 

Wave  direction-visual  estimate  (deg) . 

WH 

8 

Wave  height-visual  estimate  (ft.) 

OVERVIEW  OF  WAVE  RECORDS, .1973 
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OVERVIEW  OP  WAVE  RECORDS,  1974 
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OVERVIEW  OF  WAVE  RECORDS,  1974 
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TABLE  E-II 


WAVE  STATISTICS  FOR  1973  & 1974 


nRMS 

RMS  value  for  each  probe. 

—2 

'rms 

Average  mean  square  value 
four  probes. 

for  the 

SE 

'1 

Standard  deviation  of  the 
square  value. 

mean 

ffv= 

Average  significant  wave  ] 

height. 

sfi 

Standard  deviation  of  the 

RMS. 

NOT  DIGITIZED 

The  data  corresponding  to  the  sequence 
number  was  not  digitized.  It  was  often 
either  non-existent  or  defective  and 
therefore  not  processed. 
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03  AUG 


28  AUG 


05  SEP 


IV 

ft. 

ft. 

ft. 

192 

n 

.0036 

0.665 

.0301 

296 

ESI 

.0093 

1.010 

.0463 

162 

.02  8 

.0017 

0.571 

.0246 

248 

.065 

.0060 

0.399 

.0566 

251 

.075 

.0058 

0.973 

.0565 

302 

. 100 

.0063 

1.120 

.0423 

244 

.059 

.0024 

0.856 

.0133 

331 

.110 

.0056 

1.179 

.0287 

357 

.126 

. 0074 

1.222 

.0422 

134 

.018 

.0046 

0.443 

.0463 

286 

.089 

.0152 

1.002 

.0625 

318 

.108 

.0055 

1.166 

.0318 

362 

.142 

.0082 

1.338 

.0535 

246 

.058 

.0048 

0.804 

.0322 

286 

.080 

.0107 

0.989 

.0842 

384 

.141 

.0118 

1.357 

.0556 

247 

.059 

.0033 

0.811 

.0270 

243 

.053 

.0032 

0.786 

.0269 

316 

. 100 

.0068 

1.149 

.0401 

264 

.070 

.0033 

0.915 

.0191 

32  3 

.101 

.0030 

1.099 

.0177 

206 

.042 

.0039 

0.708 

.0331 

228 

.055 

.0041 

0.623 

.0199 

332 

.106 

. 0067 

1.127 

.0321 

300 

.090 

.0032 

1.046 

.0245 

262 

.070 

.0043 

0.889 

.0286 

361 

. 128 

.0043 

1.256 

.0268 

206 

. 039 

.0028 

0.650 

.0200 

216 

.045 

.0026 

0.742 

.0206 

182 

.033 

.0017 

0.611 

.0084 

302 

.087 

.0033 

1.007 

.0132 

122 

.014 

.0015 

0.356 

.0192 

159 

.017 

.0055 

0.402 

.0519 

199 

.040 

.0033 

0.689 

.0340 

131 

.015 

.0011 

0.379 

.0090 

263 

.066 

.0035 

0.858 

.0216 

328 

. 104 

.0077 

1.119 

.0448 

320 

.107 

.0110 

1.157 

.0687 

210 

.042 

.0060 

0.704 

.0463 

262 

.069 

.0073 

0.906 

.0376 
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TABLE  E-II 


WAVE  STATISTICS  FOR  1973 


01  NOV 


.202  .194 
.252  .235 
.300  .267 


NOT  DIGITIZED 


0.373 

0.670 

1.255 

1.348 

0.621 

0.790 

1.190 

1.298 

0.477 

0.614 

0.860 

1.060 

0.664 

0.869 

0.992 

0.629 

0.765 

0.443 

0.684 

1.046 

0.843 

1.359 

1.252 


.0179 
.0348  2 

.0206  3 

.0231 

.0353 
.0302  2 

.0305  3 

.0554 

.0267 
.0252  2 

.0362  3 

.0277 

.0217 
.1134  2 

.0343  3 

.0208 
.0218  I 2 

.0214  3 


.0279 
.0577  I 2 
.0397  3 


300 


TABLE  E-III 

OVERVIEW  OF  MODEL  TESTS  IN  SAN  FRANCISCO  BAY  FOR  1973  AND  1974 


Abbreviations 


Heading  B 

Beam 

(Relative 

wave  heading 

Q 

Quartering 

with  respect  to  the  model 
heading,  obtained  visually) 

F 

Following 

H 

Head 

GM  No . - See  Table  A- 1 for  the  SL-7  model. 


Model  Heading  - Degrees  magnetic  (approximate) 


Capsize  S 


Capsize  to  starboard 


P - Capsize  to  port 


Movies  X 


Film  records  were  taken  but  not  necessarily 
of  the  entire  run 


4-STBD  - The  autopilot  was  "stepped"  in  two  degree 

increments  four  times  altering  the  course 
to  starboard.  These  corrections  were  re- 
quired to  keep  the  model  in  the  same 
relative  encounter  heading  with  the  waves. 


No.  of  Samples  - Number  of  digitized  samples. 
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TABLE  E-III 

MOTION  DATA  RECORDED  FOR  THE  SEAL AND- 7 MODEL 
DURING  THE  1973  SEASON. 


03  AUG  - 
10  AUG 


I 

X 


1* 


5 

6 


PS 

Q 

s 
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w 
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►4 

Q 

> 

>2 
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EH 

►4 
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< 

< 

PS 

Pi 

H 

O 

D 

w 

D 

(O 

Pi 

Pi 

PS 

s 

W 

to 



3 

4 

5 

6 

7 

_ 

14  AUG  - 
27  AUG 
(Run  11.07) 


(Run  11.08) 
27  Aug  - 
14  SEP 


2 3 4 


6 7 


17  SEP  - 
28  SEP 


2 3 


4 5 fi  7 


04  OCT 
01  NOV 


2 3 4 


NO  MOTIONS  RECORDED 


6 7 


* The  number  refers  to  the  channel  on  which  the  data  was 
recorded. 
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TABLE  F-III 


Significant  changi 


15  AUG  6.06  15:26  94  Q 2 4 P X 8-PORT  1156  26  6 224 

6.27  .’.5:34  358  F 2 4 P X « 27  7 893 


17  AUG | I | | I II  III  Stbd  catamaran  angina  blew 
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OVERVIEW  OF  MODEL  TESTS  IN  SAN  FRANCISCO  BAY 


02  115  STBD  24  X Manual  control,  rudder 

Turn  hard  over  • 44  2 272 

03  13:10  55  Q 3 4 S X ’ 45  3 107 

04  13:22  45  Q 3 4 S X " 46  4 67 

05  13:58  900  Q 1 4 33/5  X 18-STBD  Le«»  wind  at  end  o£  " 47  5 204S 


OVERVIEW  OF  MODEL  TESTS  IN  SAN  FRANCISCO  BAY 
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OVERVIEW  OF  MODEL  TESTS  IN  SAN  FRANCISCO  BAY 


OVERVIEW  OF 


EXPERIMENTS  WITH  MODEL  SL-7,  HEAVY  CONDITION,  1974 
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